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A stack of foils of very pure iron was bombarded with 
deuterons of energy 10 Mev. Excitation curves are thus 
obtained for a radioactive manganese of half-life 21 minutes 
made by a (D,a) reaction and a radioactive cobalt of half- 
life 18 hours, made by a (D,n) reaction. These radioactive 
' jsotopes are both derived from the same parent iron iso- 

. Hence the radioactive yields reveal the relative 
probability of the escape of a neutron and the escape of an 


alpha-particle from the same excited nucleus. As the 
energy of the incident deuterons is increased the intensity 
of the 18-hour (D,n) cobalt attains a maximum value at 
about 7.5 Mev. There is a progressive increase with energy 
in the activity of the 21-minute (D,a) radioactive man- 
ganese. At 10 Mev the escape of a neutron is about three 
times as probable as the escape of an alpha-particle. 


HEN iron is bombarded with deuterons, 
radioactive isotopes of manganese and 
cobalt are formed. Two of these radioactivities 
have half-lives of 21 minutes and 18 hours, 
respectively, and decay by the emission of posi- 
trons. They have been assigned! with a reasonable 
degree of certainty to Mn® and Co, respec- 
tively, both being derived from the same parent 
isotope, Fe**, whose abundance in iron is only 6.5 
percent. These reactions are therefore as follows: 


oF + + + on! 
or 


It was previously suggested that this would 
allow a unique observation of the relative escape 
probability for a neutron and for an alpha- 
particle from the same excited nucleus and a 
preliminary report? of results was made. 

In the present investigation a stack of foils 
of very of very pure iron was exposed to the deuteron 
{010 (1937); J, Livingood and G. Seaberg, 
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eM Cork and B. R. Curtis, Phys. Rev. 53, 946 
1938); 55, 1264 (1939). 


beam of range 64 cm in air. Each foil had a 
stopping power equivalent to about 3 cm of air, 
experimentally determined. After a bombard- 
ment of about 100 microampere minutes, at a 
current of 15 microamperes, the activity in each 
foil was studied by the use of a freon-filled ioniza- 
tion chamber connected to a string electrometer. 
With pure iron, no other radioactivities are pro- 
duced in sufficient amounts to confuse the 
analysis of the decay curves. It was thus possible 
to determine the mean energy for the deuterons 
passing through any foil and the total yields of 
the cobalt isotope of half-life 18 hours and the 
manganese isotope of half-life 21 minutes in 
every foil. 

These yields were first expressed in terms of 
divisions per second of the electrometer at the 
end of the bombardment. It is now possible to 
convert these observed intensities actually to 
the total number of activated nuclei formed. This 
was done by collecting on a small clean metallic 
surface a small quantity of actinon gas. In de- 
caying to Pb?” every actinium B atom emits one 
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life of 36 minutes associated with the beta-decay 
of actinium B is long compared to the half-lives 
of the subsequent emitters. The alpha-particles 
could be completely counted by an appropriate 
ionization chamber and linear amplifier. This 
then became an exact source of beta-particles, 
which could be used to calibrate the electrometer. 
It is assumed that the energies of the beta- 
particles are such that the number that do not 
enter the thin window and completely traverse 
the ionization chamber is negligibly small. It is 
also necessary to cover the actinon with a very 
thin absorbing layer sufficient to prevent the 
escape of the heavy nuclei by recoil upon dis- 
integration, yet not enough to stop the 6.6-Mev 
alpha-particles. 

The data obtained with such a stack of foils, 
with the energy expressed in Mev are shown in 
Fig. 1. Four such sets of foils were studied giving 
in each case results similar to that shown. The 
ordinates, while labeled relative cross sections, 
may be made absolute by the fact that at the 
number 100 the cross section is about 1xX10-” 
cm?. The cover foil was in each case rejected as it 
appeared to have too little activity. This is due 
to the fact that some of the recoil atoms from 
every foil are driven into the next underlying 
foil. Even after rejecting the cover foil there still 
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Fic. 1. Excitation curves for (A)Fe(D,n)Co and 
(B) Fe(D,a)Mn expressed in terms of cross section. At 100, 
equals 1 X10-*? cm’. 
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particle. 

seems to be a very slight increase in the 18-hour 
cobalt activity in going to deeper foils before g 
rapid decrease is observed. No such saturation was 
observed for the 21-minute manganese activity, 

In Fig. 2 is shown the ratio of the yield of the 
18-hour (D,n) cobalt to the corresponding yield 
for the 21-minute (D,a) manganese. In both 
processes the deuteron as a whole has entered the 
nucleus through the potential barrier; that js 
the Oppenheimer-Phillips splitting of the dey. 
teron in the external field would play no part in 
these reactions. If the masses of the incoming 
deuteron and its incident energy be compared 
with (a) the mass of the escaping neutron and 
(b) the mass of the escaping alpha-particle, and 
if we assume the same nuclear defect for original 
and final atoms the limiting energy of the escap. 
ing particles can be determined. For an incident 
deuteron of 10 Mev the neutron could have a 
maximum energy of 15 Mev and the alpha 
particle an energy of 20 Mev. 

It may perhaps be reasonable that the ob 
served yield of 18-hour (D,n) cobalt attain a 
maximum value and then fall off slightly as the 
energy of the incident particles is increased. This 
would be the case, if with increasing energy the 
probability of the competing (D,a) process be 
came greater and the total number of excited 
atoms remained the same. Since the incident 
deuteron energies are much in excess of the 
nuclear barrier it might be expected that the 
number of excited nuclei would be constant 
At 10 Mev the possibility of a neutron being 
emitted is about three times as probable as the 
emission of an alpha-particle. 

This investigation is made possible by a grant 
from the Horace H. Rackham Fund. 
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The cross sections of a number of elements for neutron scattering were measured as a function 
of neutron energy throughout a range of 0.5 Mev. Anomalous cross sections were found in the 


cases of magnesium and aluminum. It was possible to account quantitatively for the anomaly 


in magnesium on the assumption that the effect was due to resonance in Mg®. 


INTRODUCTION 


HE first measurements of nuclear cross sec- 

tions for fast neutron scattering were made 
by Dunning in 1934.' Using neutrons from a 
mixture of radon and beryllium whose energies 
ranged from about 0.1 to nearly 14 Mev, he 
concluded that the cross section of a nucleus of 
atomic weight A was proportional to A!. This 
cross section was, of course, an average with 
respect to neutron energy. More recently, Laden- 


burg and Kanner,* and Kikuchi and Aoki* 


measured the cross sections of a large number of 
nuclei, using the nearly homogeneous neutrons 
produced in the deuteron-deuteron reaction. 
They found that the cross section for these 
neutrons was not proportional to A!, and in fact 
was not a monotonic function of the atomic 
weight. Evidently then, even for fast neutrons, 
nuclear cross sections must have been appreci- 
ably dependent on neutron energy. It therefore 
seemed worth while to measure throughout as 
wide an energy range as possible the cross sec- 
tions of a number of elements in order to verify 
the expected energy dependence, and, more im- 
portant, to see if there were anomalous scatter- 
ing. The experimental results were reported some 
time ago.* Aoki‘ had previously reported a similar 


‘effect in silicon, and recently in more extended 


experiments’ found a number of additional 


PR: er with the Standard Oil Company of Venezuela, 
i; Dunning, Phys. Rev. 45, 586 (1934). 

937) Ladenburg and M. H. Kanner, Phys. Rev. 52, 911 
% S, Kikuchi and H. Aoki, Sci. Papers of the Inst. of 

Phys. and Chem. Research 34, 865 (1938). 

1999) MacPhail and J. Giarratana, Phys. Rev. 56, 207 
‘H. Aoki, Phys. Rev. 55, 795 (1939). 
°H. Aoki, Proc. Phys. Math. Soc. Japan 21, 232 (1939). 


anomalies. Other work has also been reported by 
Zinn, Seely, and Cohen.‘ 


THE APPARATUS 


The neutrons were produced by the action of 
100-kv deuterons on targets of heavy ice. The 
deuteron source and discharge tube were similar 
to the source and discharge tube described by 
Zinn and Seely.”| * The ion beam was unanalyzed 
because it was thought undesirable to have a 
contamination target of variable strength formed 
where the discarded beam would hit the walls of 
the discharge tube. Ion currents used were from 
600 to 1000 microamperes, the upper limit being 
set by the amount of heat which could be dissi- 
pated at the target. This was a cylindrical copper 
box cooled with liquid air, on the front face of 
which heavy water was condensed. The 100-kv 
accelerating voltage was supplied by an x-ray 
transformer-rectifier set. 

For counting neutrons, a system consisting of 
a paraffin-lined ionization chamber, linear ampli- 
fier,* and thyratron counter” was used. The 
ionization chamber and linear amplifier were 
mounted on a light table which could be easily 
rolled about on wheels. A Lauritsen electroscope™ 
lined with paraffin was used to measure the 
neutron intensity. | 

The arrangement of the apparatus is shown 
schematically in Fig. 1. The target, scatterer, 
ionization chamber, and electroscope are repre- 
sented, respectively, by T, S, J, and E. Scatterers 
were cylindrical and were slightly larger in di- 

* Zinn, Seely, and Cohen, Phys. Rev. 56, 260 (1939). 
7™W. H. Zinn, Phys. Rev. 52, 655 (1937). 

® W. H. Zinn and S. Seely, Phys. Rev. 52, 919 (1937). 
* J. R. Dunning, Rev. Sci. Inst. 5, 387 (1934). 


10 J, Giarratana, Rev. Sci. Inst. 8, 390 (1937). 
"C. C. Lauritsen and T. Lauritsen, Rev. Sci. Inst. 8, 


438 (1937). 
669 


APRIL 15, 1940 

18-hour 

before a 
ation was 
Ctivity, 
id of the 
ing yield 

In both 
tered the 

the deu- 
part in 

incoming 
>ompared 
tron and 
ticle, and 
Original 
he escap- 
1 incident 
d havea 
1e alpha- 
t the ob- 
attain a 
tly as the 
ased. This 
nergy the 
rocess be 
of excited 
incident 
ss of the 
| that the 

constant. 
ron being 
ible as the 
by a gramt 


670 


DEUTERONS 


%, 


Fic. 1. Arrangement of neutron source 7, electroscope E, 
scatterer S, and ionization chamber J. 


ameter than the ionization chamber (4.2 cm) and 
between 4 and 6 cm in length. 

The neutron energy was varied by rotating the 
counter and scatterer about the target, that is, 
by varying the angle @ between the paths of the 
neutrons and the deuterons. The following rela- 
tion obtained from the equations of conservation 
of energy and momentum gives the dependence 
of neutron energy E, on the deuteron energy Ep 


and angle 6: 
E,,— (EpE,/2)* cos 6— (En+3Q)/4=0. 


Q, the energy released in the reaction D?+D?* 
—He’*+n'!+(Q, was taken to be 3.32 Mev.” The 
dependence of the neutron energy on @ is shown 
graphically in Fig. 2 for deuteron energies of 
0.10 and 0.05 Mev which are the energies of the 
D+ and D,* ions, respectively. Because of the 
presence of both kinds in the beam and because 
of the slowing down of the deuterons in the ice 
target, the neutrons will have a spread of energy 
corresponding to that of the deuterons from 0 
to 0.10 Mev. 


EXPERIMENTAL PROCEDURE 


The cross section for a particular scatterer was 
calculated from the measured transmission, i.e., 
the ratio of the number of neutrons indicated by 
the counter with the scatterer in front of it to 
that indicated with the scatterer removed, the 
times being such that the electroscope deflections 
in each case were the same. A correction was 
made for the number of neutrons deflected by 
the scatterer into the ionization chamber on the 
assumption that the scattering was isotropic. 


2, G. Bonner, Nature 143, 681 (1939), 
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For the dimensions indicated in Fig. 1 and q 
measured transmission ¢, this was 0.01(1 —t) and 
was to be subtracted from ¢. A second small 
correction due to the finite width of the neutron 
energy spectrum will be given later. From the 
corrected transmission to thus obtained, the crogs 
section o was calculated by means of the formula 
o=(MA/Nm) log (1/to) where A is the area of a 
face of the scatterer, m its mass, M its atomic or 
molecular weight, and N Avogadro's number. 
The exponential law of absorption on which this 
formula rests was tested by measuring the trans. 
mission as a function of length through scatterers 
of aluminum. The results of two tests are given 


Fig. 3. 


The intensity of scattered neutrons coming 
from the walls and floor of the room was esti. 
mated by measuring neutron intensity as a 
function of the distance D of the counter from 
the target. The results of two tests made at 0° 
and 90° to the deuteron beam are shown graphi- 
cally in Fig. 4. All measurements of transmissions 
were done at a distance of 56 cm from the target 
for which D* = 320 X10~* cm~. The intercept at 
D-*=0, i.e., the proportion of neutrons scattered 
by the walls and floor, is less than 3 percent of 
the number at 56 cm coming from the target. 
The effect of these scattered neutrons would be 
to decrease the cross section from 0 to 3 percent 
depending on their intensity. But as the intenisty 
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Fig. 2. Neutron energy as a function of the angle of 
emission for a thin target. 


SCATTERING OF 


of the scattered neutrons is probably about the 
same for all values of 6, a correction, if it could 
be made, would only increase the cross sections 
at each energy by a small amount. 
Because of the finite width of the energy 
trum and the finite angle subtended at the 
target by the neutron counter, the measured 
transmission is an average value, the average 
being taken first with respect to neutron energy 
and then over the face of the ionization chamber. 


A simple calculation gives 


=to(Eo) 
1 dty d*to 
h?+ ? 
+(h?+k?) 


@ 
4 dE» 


where (6) is the transmission measured at angle 
@ (the angle between the direction of the deu- 
terons and the line through the target and the 
center of the ionization chamber), to(Eo) the 
transmission for homogeneous neutrons of energy 
Ey, Eo the average energy of the neutrons 
emitted at angle 0, » the angle subtended by the 
radius of the ionization chamber at the target, 


a 


2 3 
SCATTERER THICKNESS IN CM 


Fic. 3. Experimental test of the exponential law of 
absorption for the case of aluminum. 


h=}n4E,/d0 the root-mean-square variation of 
neutron energy across the ionization chamber, 
and k the root-mean-square variation of the 
neutron energy about the average, Eo, due to the 
lack of homogeneity of the neutron beam at a 
given angle 6. The term in square brackets is to 
a first approximation the effect on the transmis- 
sion of finite resolving power determined by h 
and k. We shall neglect terms beyond it which is 
equivalent to assuming that the transmission- 


FAST NEUTRONS 


NEUTRON INTENSITY 


.0002 .0003 0004 
INVERSE SQUARE OF DISTANCE FROM TARGET IN cm™® 
Fic. 4. Testing of the inverse square law; if scatteri 


from the walls of the room were appreciable, the lines woul 
not pass through the origin. 


energy curves show no fine structure. Replacing 
to by ¢ in the square brackets, we can then solve 
(1) for ft) in terms of the measurable quantities, 
t(@) and its first two derivatives; i.e., 


dt 1 d*t 
| (2) 


o(Eo) — 
de? dE, dE,? 


The first term in the square brackets turns out 
to be 10~* or less, so is negligible. The other 
term represents the combined effects of an ioniza- 
tion chamber of finite size and of inhomogeneous 
neutrons. It is seen to sharpen the transmission 
curves by raising the peaks and lowering the 
hollows where d*t/dE,? is, respectively, negative 
and positive. This is of course what one would 
expect. 

To obtain E» and k as functions of @, it is 
necessary to calculate the neutron spectrum. 
Thus, if the number of neutrons emitted at angle 
6 with energies between E and E+dE is f(E)dE, 
then 


DE,= f 
0 


and 


Dk? = f [>- f | 


The function f(£) was calculated on the assump- 
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Fic. 5. Distribution in energy of the neutrons somes off 
the target in the forward direction calculated for 100 kv 
accelerating potential for the cases (a) pure molecular 
beam, (b) pure atomic beam. 


tion that the ice target was thick (an effort was 
made always to fulfill this condition experi- 
mentally), and that the unanalyzed beam of 
deuterons was 20 percent atomic and 80 percent 
molecular. In addition, the 1 Mev group of 
neutrons reported by Bonner” and others was 
not taken into consideration because of general 
uncertainty regarding it. The excitation function 
used was that given by Baldinger, Huber, and 
Staub." By way of an example, the neutron 


13 Baldinger, Huber, and Staub, Helv. Phys. Acta 11, 
245 (1938). 


TABLE I. Average corrected 


spectrum calculated in this way for @=0° jg 
given in Fig. 5. Values of Ey and k for @=0° and 
certain additional angles are listed in the follow. 
ing table: 

4 0° 76° 92° 109° 4990 
Eo(Mev) 2.80 2.65 249 241 2.34 
k (Mev) 0.065 0.035 0.020 0.009 0.010 0.025 

The other quantities in square brackets in (2) 
namely, h=ndEo/2d0, d*Eo/d6, dt/dEs, and 
d*t/dE,, were found graphically. As mentioned 
before, the correction was entirely due to the 
term —}(h?+k?)d*t/dE,’. In absolute magnitude 
it was in every case between 0 and 0.008; that is, 
it was of the order of 1 percent or less of ¢. 

The cross section for each material and for 
each chosen neutron energy was calculated from 
the corresponding transmission after the correc. 
tions mentioned, above had been applied. This 
value of the transmission was the average of from 
five to twenty measurements made on different 
days. The probable error of the average trans. 
mission was calculated from the deviations of the 
individual values about their means. In Table | 
are given as functions of neutron energy the 
average transmission (corrected), its probable 
error, and the corresponding cross section in 
units of 10-*4 cm?. 

The cross sections given in Table I are plotted 
against neutron energy in Fig. 6. The coordinates 
of the center of each rectangle are the most 
probable values of the cross section and the 


transmission coefficients, their probable error (P.E.) and cross section in units 
of 10-* cm* as functions of neutron energy. 
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SCATTERING OF 


corresponding energy. The height of the rec- 
tangles for a given element is the experimental 
probable error in ¢. This was calculated from the 
average observed probable error of the corre- 
sponding transmissions. Widths of the rectangles 
indicate the amounts by which the neutron 
energy would change due to a change in propor- 
tion of atomic deuterons in the beam from 10 to 
30 percent about a mean of 20 percent. The side 
of the rectangle corresponding to 10 percent 
atoms is the one nearer the energy of 2.50 Mev. 

There are two additional possible sources of 
error in the neutron energies: 


1. An error dQ in the Q value for the D—D 
reaction would cause the cross-section curves to 
be shifted horizontally by approximately 3dQ/4. 
Bonner’s estimate for dQ is +0.04 Mev.” 

2. Fluctuations of the accelerating voltage. 
This voltage was measured frequently and be- 
lieved to be 100+-5 kv. The effect on the neutron 
energy would be similar to but smaller than that 
produced by fluctuations in the proportion of 
atoms in the deuteron beam. 

The results of Aoki® and of Zinn and others*® 
for the elements reported here have been added 
to Fig. 6 after they were recalculated with 
Q=3.32 Mev and allowances made for the 
widths of the neutron spectra. Each of their 
values is represented by the chemical symbol of 
the element, the former’s surrounded by a circle 
and the latters’ by a triangle. 


A THEORY OF ANOMALOUS SCATTERING OF 
Fast NEUTRONS 


We shall now attempt to discuss on the basis 
of present theory the anomalous behavior of the 
magnesium cross section around 2.5 Mev and 
later make a few remarks about the case of 
aluminum. 

What happens when magnesium is bombarded 
by fast neutrons is complicated by the presence 
of three isotopes of masses 24, 25, and 26 and 
respective proportions 0.78, 0.11, and 0.11; but 
because of the greater abundance of Mg", we 
shall attribute the anomaly in cross section to an 
excited level in the intermediate Mg* nucleus 
formed from Mg*‘. The possible reactions to be 
considered are elastic and inelastic scattering and 


FAST NEUTRONS 


the following disintegrations and capture: 


1. Mg**+n'—Na**+H'—4.0 Mev, 
2. Mg**+-n'—Ne*!+ Het —2.3 Mev, 
3. (7.0 Mev). 


With 2.5-Mev neutrons, the first disintegration 
is energetically impossible, and the second highly 
improbable because of the small amount of sur- 
plus energy for the a-particle. The third can be 


Fic. 6. Observed scattering cross sections as a function 
of neutron energy. The rectangles represent the nt 
observations, their heights and widths correspon ing to 
the range of error in cross section and energy determina- 
tions. The circles give the results of Aoki, the triangles 
a of Zinn, Seely, and Cohen (both recalculated; see 
text). 


neglected too, being about 10~‘ times as probable 
as scattering." 

Inelastic scattering also is apparently small. 
The cross section for this process has been given 
by Kikuchi and others” as 0.08X10-** cm? or 
only about 3 percent of the total cross section. 
We therefore conclude that the only important 
process to be considered is elastic scattering by 
the nucleus Mg**. 

The cross section for elastic scattering of par- 
ticles of zero spin by nuclei also of zero spin is"® 


o (2n+1) sin? (3) 


where X=h/(2ME)'. M is the reduced mass of 
the incident particle, E its kinetic energy, and n 
its orbital quantum number in the field of the 
nucleus. The 4’s are phase shifts of the outgoing 
or scattered waves with respect to an incident 
wave unaffected by the scattering nucleus. If, 
xuth f- Bethe, Rev. Mod. Phys. 9, 160 (1937), Table 

Kikuchi, Aoki, and Husimi, Proc. Phys. Math. Soc. 
Japan 18, 115 (1936). 


N. E. Mott and H.S. W.sMassey, The Theory of Atomic 
Collisions, p. 24. 
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Fic. 7. If the elastic scatterin cross section is multiplied by neutron energy, theory predicts in certain 
simple cases that the result is a function_of energy which is composed of a constant part and a part which 
varies with energy as y=sin* 5, where 5 =5-—arc cot (E—E»)/4l'=5—arc cot x. The theoretical curves given 
here for y as a function of x are compared with the observations to find the value of 5 which gives the best 


agreement. 


however, the incident particles have spin 3, the 
nuclei still having spin 0, the cross section will be 


o=4rx? (n+}) sin? (4) 


s=+}, n=0, 1, 


For orbital quantum number n=0 there is only 
one phase shift to consider, 59= do, ; = do, -:. 

The phase shifts occurring in (4) are in prin- 
ciple determinable by solving the wave equation 
for the compound nucleus temporarily formed in 
the scattering process. We are interested in the 
case when this compound nucleus possesses a 
certain semi-stable level which is narrow in 
comparison both with the kinetic energy of the 
neutron and with the distance between levels. 
We denote by 2, s the orbital and spin quantum 
numbers of the neutron which is required to form 
this state from the original nucleus, by E,, the 
energy which the incident neutron should have 
to be in resonance with this state, and by Ins 
the width of the resonance level. The mean life- 
time, 7, of the compound state is connected with 
Ine by the relation T=h/T'ns. According to the 
theory of nuclear scattering,'” the neutron 
wave corresponding to the given orbital and spin 


17G, Breit and F. L. Yost, Phys. Rev. 48, 203 (1935). 
#8 P, L. Kapur and R. Peierls, Proc. Roy. Soc. 166, 277 


(1938). 


quantum numbers will be shifted in phase with 
respect to an unperturbed wave by an amount 
which depends on the neutron energy E and 
which is given by 


= —arc cot (E Ens) ne (5) 


provided the neutron energy does not come into 
the neighborhood of another resonance level. The 
constant 6,,. expresses the influence of all other 
levels of the compound nucleus which belong to 
the given spin and orbital quantum numbers. 
It is evident from (5) what will happen as the 
neutron energy is increased. Considerably below 
the resonance energy we shall have a phase shift 
5ns= dns. As E increases, the phase shift rises more 
and more rapidly until we reach the center of the 
resonance level, where E = En, and Sne=dne+27/2; 
the phase shift then continues to increase, but 
more and more slowly, reaching finally the value 


5ns+7. This is true of course only of the phase - 


shift corresponding to the given orbital and spin 
quantum numbers; the variation with energy of 
the other phase shifts may be neglected if we are 
dealing with the scattering in the neighborhood 
of a single isolated resonance level. The actual 
scattering cross section, according to (4), is the 
product of the slowly varying factor 47%? (in- 
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versely proportional to neutron energy) and a 
sum of terms due to neutrons of various orbital 
and spin quantum numbers. Of these terms only 
one will change appreciably through the reso- 
nance region and be responsible for the variation 
of cross section with energy. This term is (except 
for the constant factor »+}) 


y=sin? dns. 


It may either increase or decrease as we approach 
the resonance level, according as the constant 
j,, happens for the particular nucleus to be less 
than or greater than 7/2. In fact, from the 
manner in which the cross section varies near 
resonance we can even estimate the actual value 
of 5,2 with some accuracy. Curves are drawn in 
Fig. 7 showing the variation with energy of the 
term y=sin? 6,, for a number of values of dns. 
The independent variable in these curves is 
x=(E—Ens)/30 ne; since it varies linearly with 
enerty, the theoretical curves can be compared 
as regards shape (quite without any reference to 
the absolute magnitude of the scale of independent 
and dependent variables) with a plot of the 
observed cross section as a function of energy. It 
was found in this way that 5,,.=—60° gave a 
theoretical variation of cross section with energy 
in closest accord with the observations on mag- 
nesium reported here. 

By comparing on an absolute scale the theo- 
retical and observed scattering, we can determine 
in addition to 65,. (i) the value of the quantum 
number (ii) which appears in the statistical 
factor n+}, the energy E,, (iii) and width I, (iv) 
of the resonance level, and the essentially con- 
stant contribution (v) to the scattering due to 
deflection of neutrons with other spin and orbital 
quantum numbers. Conversely, we may say that 
we have four adjustable constants and our choice 
of a certain integer m with the help of which to 
fit the theoretical formula to the observations. 
To determine , we note that the factor y varies 
as we go through resonance between a minimum 
value of 0 and a maximum value of 1, so that 
the scattering cross section changes by the 
amount 


if 
4nX*(n+}3) if (6) 


(provided we neglect the slight variation of x 


n=0, 
n=1,2,--- 
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itself with energy) its minimum value being 
determined by the contribution (v). 

We shall apply (6) to the experimental results 
for magnesium, but must first multiply 4X? by 
the relative abundance 0.78 of the isotope Mg**. 
With a neutron energy of 2.55 Mev we have 
0.78 X 41x? = 0.84 X 10-*4 cm?. Therefore the total 
variation in cross section through resonance 
will be 


0.84X10-* cm? if n=0. (7) 
0.84 10-24(2n+1)/2 if (8) 


according to theory. The observed value taken 
from Fig. 6 is 0.80 10-*4 cm*. Thus all values of 
n other than n=0 seem to be entirely excluded. 
On the other hand, with n=0, we obtain a real 
check between theory and observation. 

The constant 6,, is to be designated as dy in 
the present case and has already been estimated 
to be —60°. The resonance energy Ep and level 
width Ip are conveniently found from the energies 
E, and E; at which the observed cross section 
reaches, respectively, its maximum and minimum 
values by the help of the relations 


Ey= cos 25o, (9) 
To= (E2 sin (10) 


which follow directly from Eqs. (4) and (5). 
Observing that E;=2.67 Mev and E,=2.50 Mev 
(see Fig. 6) we complete our determination of the 


2, coe, 


CALCULATED MAGNESIUM 
CROSS SECTION G*x 


2200 
NEUTRON ENERGY IN MEV 


_ Fic. 8. The vertical lines represent the present observa- 

tions on magnesium for six different energies and are 
compared with the smooth theoretical curve calculated 
with the best values of the four adjustable parameters 
discussed: in the text. These meters allow no freedom 
as to the vertical scale of the theoretical curve; for the 
increment in cross section between minimum and max- 
imum theory gives 0.84X10-* cm?, tion gives 
0.80 X 10-* cm?. 
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TABLE II. The complexity of the states of the compound 
nucleus Al*. A few of the lowest states are classified according 
to parity, total angular momentum, orbital angular momen- 
tum and multiplicity. 


TOTAL ANGULAR MOMENTUM 
PARITY 0 1 2 


Even | °D 'G 


parameters which enter the theoretical formula 
for the cross section near resonance. They are 
four in number (since »=0 is not really an ad- 
justable constant) : 


(i) —60°, 
(iii) Mev, 
(iv) To=0.15 Mev, 
(v) Cross section’s minimum value 
=1.76X10-*4 (11) 


The value found for I corresponds to a mean life 
of the given level of the compound nucleus of 
4.5X10-*! sec. On the other hand, the time for a 
2.5-Mev neutron to travel a nuclear diameter 
(say 10- cm) is 4.5X10-" sec. That is, the 
neutron spends about 10 times as long inside the 
nucleus as it would take to travel straight acr 

it with no interaction. 

With the constants given in (11), the cross 
section was calculated as a function of energy. 
In this calculation the variation with energy of 
the factor 47x? in Eq. (4) was taken into ac- 
count. The results are plotted in Fig. 8. The six 
experimental values indicated by the vertical 
strokes are added for reference and show satis- 
factory agreement with the theoretical formula 
containing the four adjustable parameters. 

The conclusion, based on the foregoing inter- 
pretation, is that the anomalous scattering found. 
in magnesium is due to an S level in Mg* at a 
height of 9.5 Mev above ground. This height was 


calculated from the masses of the neutron, of 
Mg”, and of Mg” as well as from the neutron 
energy at resonance. 

Because there is but one aluminum isotope 
Al??, the situation is much better in this ; 
than for magnesium. The problem is however 
greatly complicated by the large nuclear spin 
5/2, of Al’. The simplification in the case of 
magnesium was due to the fact that only one 
orbital quantum number 1 could be assigned to 
a state of given parity and given total angular 
momentum. On the other hand, the complexity 
of the states of the compound nucleus, Al, js 
indicated in Table II in which a few of the lowest 
ones are classified according to parity, total 
angular momentum, orbital angular momentum, 
and multiplicity (given by the superscript). The 
scattering cross section for a particular state jg 
now not only determined by the corresponding 
values of but also by the transition probabilities 
between the different components of that state, 
These are unfortunately unknown. 

This method of treating elastic neutron scat- 


tering is accordingly limited in its applicability - 


to nuclei of spin zero. It would however be 
possible to extend it to light nuclei of spin 4 for 
which S states might be expected to be most 
probable. If »=0, the variation of cross sec- 
tion near resonance would then be given by 
3X? sin? 593 or by 7X? sin? 59: according as the 
state was a triplet or singlet. The 6’s would of 
course still be determined by Eq. (5). This is a 
case which it will be useful to test as more ex- 
perimental results are obtained. 

It is a pleasure to express my thanks to Pro- 
fessor Ladenburg, Professor Wheeler, and Dr. 
Giarratana for suggestions and help with this 
problem. 


19S. Millman and P. Kusch, Phys. Rev. 56, 303 (1939). 
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' An Electrical Quadrupole Moment of the Deuteron 
The Radiofrequency Spectra of HD and D; Molecules in a Magnetic Field* 
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Columbia University, New York, New York, 
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The molecular-beam magnetic-resonance method has 
been applied to study the radiofrequency spectra of the 
molecules D: and HD, in the state J=1. All the measure- 
ments were made in applied magnetic fields large enough 
to decouple the angular momentum vectors concerned. 
The spectrum of Dz for the transitions m;= +1 consists 
of six lines as in Hz. The HD spectrum consists of two sets, 
one of nine lines in the region of the Larmor frequency of 
the proton, and another of twelve lines in the region of the 
Larmor frequency of the deuteron. The analyses of these 
spectra into sets of energy levels and a consistent theory 
are presented. As a consequence of the analysis one has to 
assume the existence in the deuteron of an electrical 


(Received February 17, 1940) 


quadrupole moment, Q. The other quantities which are 
evaluated from these measurements are the spin-orbit 
interaction constants for the molecules concerned and the 
magnetic moments of the proton and deuteron, through 
the magnetic dipole-dipole interactions. The values of the 
constants are Q=2.73X10-" cm*, H’(D:)=14.0 gauss, 
H’p(HD) = 20.1 gauss, and H’p(HD) = 20.48 gauss. These 
spin-orbit interaction constants represent the magnetic 
field produced by the molecular rotation at the position 
of a nucleus. The values of the magnetic moments of the 
proton and deuteron agree with the direct measurements 
previously reported. 


[* a previous note! we described experiments 
which indicated the existence of an electrical 
quadrupole moment of the deuteron. In this 
paper we shall present a more complete report of 
these and further experiments and of their 
interpretation. 

It is hardly necessary to emphasize the par- 
ticular importance for nuclear physics of the 
presence of an electrical quadrupole moment in 
the deuteron. Indeed, when these experiments 
were started there was no question of investi- 
gating such a quadrupole moment, because 
current theory predicted, under the assumption 
of central forces between the proton and neutron, 
that the deuteron is in a *S, state which possesses 
no electrical dipole, quadrupole nor higher mo- 
ment. The implication of the quadrupole moment 
is that the forces between proton and neutron in 
the deuteron are not “central.” 

Since the publication of the first note on the 
deuteron quadrupole moment the question has 
received theoretical discussion by a number of 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of aN University. 

ae M. B. Kellogg, I. I. Rabi, N. F. Ramsey, Jr. and 
J. R. Zacharias, 


Phys. Rev. 55, 418 (1939). 
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writers? who have examined the problem chiefly 
from the point of view of the meson theory of 
nuclear forces. 

The possibility that a nuclear charge may de- 
part from spherical symmetry was suggested by 
Racah? and by Pauli.‘ Schuler and Schmidt® were 
the first to present clear experimental evidence 
for such an effect, which they found in the de- 
parture from the interval rule of some lines in 
the hyperfine structure spectrum of europium. 

The way in which a nuclear quadrupole mo- 
ment manifests itself is through the fact that its 
energy in an external electric field depends not 
only on its position but also on its orientation 
with respect to the gradient of the field. In an 
atom the electric field is that of the electronic 


*H. Bethe, Row, 1261 (1939). R 
and S. Kusaka hys 55, 665L 0939). B “Ings, 
Phys. Rev. 55, osgi. (1999): Ph s. Rev. 55 , 329 (1939 
Phys. Rev. 56, 1175 igner and L. 
Phys. Rev. 56, te. (1939). S. M. Dancoff, Phys. Rev. 56, 
(1939). Phys. Rev. 56, 540 (1939). 
M. E. Rose, Phys. Rev. 56, 1064 se). © Kittel and 
G. Breit, Phys. Rev 56, 744 (1939). M . Phillips, Phys. 

ev. 57, 160 (1940). 

3G. Racah, Zeits. f. Ph ik 71, J (1931). 

‘Ww. Pauli, Zurich Conference, 1931. 

a 955). Schuler and Th. Schmidt, Zeta f. Physik 94, 457 
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charge distribution. In a molecule the field arises 
from the other nuclei as wel! as from the electrons. 

The quantum theory of the effect of a nuclear 
quadrupole moment on atomic energy levels was 
given by Casimir.® A similar theory for molecules 
can be constructed along analogous lines. The 
part of the energy operator which concerns the 
quadrupole moment of a nucleus is, according to 


_ Casimir, 


IS +1)}/ 
2J(2J—1)i(2i—1). (1) 


Q is the magnitude of the nuclear quadrupole 
moment and is given by Q= (32*—r*)« and the 
average is taken over the nuclear charges for 
the state which has the largest component 
of the spin 7 in the z direction. The origin is 
at the centroid of the nucleus. The quantity 
q= (2x(3 cos? The sum is taken over 
all of the molecular charges except the nucleus 
under discussion. The electrons contribute nega- 
tively. This average is taken over the state of the 
molecule which has the greatest component of its 
angular momentum J along the z direction. R; is 
the distance of the kth element of charge from 
the nucleus in question. J and i are the molecular 
and nuclear angular momentum operators. 6; is 
the angle which R, makes with the z axis. The 
factor (ge) is the average of d°V/dz2* where V is 
the electrostatic potential at the nucleus. 

The evaluation of our experimental results will 
give us the value of the product gQ. Unlike the 
atomic case there is no empirical information 
like fine-structure splitting from which the mo- 
lecular g can be evaluated. One must? fall back 
on theoretical calculations from molecular wave 
functions to obtain g. Hence the value of Q must 
result partly from a calculation of g. 

Since the effect of artificial external electric 
fields on a nuclear quadrupole moment is not 
measurable at present because of its very small 
size, it is impossible to measure quadrupole 
electric moments as directly as nuclear magnetic 
moments. The existence of a nuclear quadrupole 
moment is therefore at present an inference from 


*H. B. G. Casimir, “On the interaction between atomic 
nuclei and electrons,”’ Prize Essay published by Teyler’s 


Tweede Genootschap (1936). 
7 However, see E. Teller, Phys. Rev. 57, 556(A) (1940). 
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the manner in which certain energy levels are 
shifted with respect to one another. In our ex. 
periments these shifts are of the same nature 
in their manifestation, as those due to the mag. 
netic dipole-dipole interaction between the two 
nuclei, and must therefore be disentangled from 
this type of interaction. However, this js not 
difficult, as we know the magnitudes of the 
magnetic moment of the proton and the deuteron 
and the internuclear distance. 

A quadrupole moment effect is relatively large 
in our experiments whereas in atomic hyperfine 
structure studies the effect is generally small, 
For atomic deuterium in the ground state the 
effect vanishes because it-is an S state for which 
qg is zero. In excited states the effect would be 
smaller than the natural width of the energy level 
due to radiation damping. In general, in the 
atomic case, the quadrupole effect is superposed 
on the large interaction between the nuclear 


magnetic moment and the magnetic field pro.. 


duced by the electrons. These fields may be of 
the order of 10° gauss. In the Dz and HD mole. 
cules which are the subject of our study these 
fields are well under 50 gauss. Since our resolving 
power is about 10° times as great as with optical 
methods we can measure the effects of even the 
small magnetic fields in the molecule. The quad- 
rupole effect is therefore relatively great because 
the gradient, ‘‘g,”’ of the electric field which 
interacts with the quadrupole moment is of the 
same order of magnitude in molecules as in 
atoms. 


METHOD 


These investigations were carried out by the 
molecular beam magnetic resonance method.’ 
The main features of this method as applied to 
gases were described in a previous paper” to 
which we will refer as I. 

All of the spectra in paper I were observed in 
the following way : The radiofrequency oscillator 
was set at a definite constant frequency fo and 
the molecular beam intensity was observed as a 


8 By resolving power we mean the smallest absolute 
interval in wave numbers in which two lines can be re 
solved, that is we mean Ay rather than Av/». ; 

9]. I. Rabi, S. Millman, P. Kusch and J. R. Zacharias, 
Phys. Rev. 55, 529 (1939). 

10]. M. B. Kellogg, I. I. Rabi, N. F. Ramsey, Jr. and 
J. R. Zacharias, Phys. Rev. 56, 728 (1939). 
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function of the value of the homogeneous mag- 
netic field H in which nogadiabatic transitions 
take place. At appropriate values of the magnetic 
field the beam intensity reaches a minimum and 
the various values of the magnetic field for which 
this occurs are called positions of the minima or 
positions of the lines. These minima occur when- 
ever the product (hfo) for the oscillator is equal 
to the difference in energy between two molecular 
states. The whole pattern of minima is referred 
to as the radiofrequency spectrum. In paper I 
we were chiefly concerned with the magnetic 
moments of the proton and of the deuteron. The 
measurements were made on HD and Dz mole- 
cules in the state with rotational angular mo- 
mentum J=0, and He molecules in the state 
J=1. In the state J=0 the HD and the D, 
molecules each give rise to a one-line radio- 
frequency spectrum, for each nucleus. In the 
state J=1 the H2 molecules give a 6-line spec- 
trum in the frequency region in the neighborhood 
of the Larmor frequency of the proton, and 
another six-line spectrum in the region of the 
Larmor frequency of the rotational moment of 
the molecule. Only the first of these spectra was 


considered in I. 


ResuME oF RESULTs 


In this paper we shall report on the nuclear 
radiofrequency spectra of HD and Dz molecules 
in the state J=1, and on our analysis of these 
spectra. Since the treatment of Dz as given here 
is so much like that of He as given in I, and since 
we wish to make use of certain of the He: results 
in the discussion of the analysis of the D2 spec- 
trum, we give here a brief resumé of the impor- 
tant features and results of the He analysis. 

The H; experiments were performed on a beam 
of molecules emerging from a liquid-nitrogen- 
cooled source. Parahydrogen molecules at liquid- 
nitrogen temperatures are practically all in the 
state J=0 and have the total nuclear spin J=0. 
They are inert as far as our experiments are 
concerned, so that our study of He was made on 
orthohydrogen molecules in the state J=1, J=1. 
In this state there are (2/+1)(27+1) =9 energy 
levels. Because of the selection rule Am;=+1 
for transitions with Am,=0, there are 6 transi- 
tions allowed between thése levels giving rise to 
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the six-line nuclear spectrum. We were able to 
analyze this six-line radiofrequency spectrum 
and thereby obtain values for two important 
physical quantities: 1, the value of (up/r*). 
(called H”) which appears because of the dipole- 
dipole magnetic interaction between the two 
protons; and 2, the magnetic field H’ produced 
by the molecular rotation at the position of a 
nucleus. More properly H’ is the constant in the 
expression for the spin-orbit interaction of the 
proton moments with the molecular rotation. In 
high magnetic fields with correspondingly high 
oscillator frequency, the six lines for He are 
symmetrically disposed with respect to Hy (de- 
fined by Hyp=hfo/2up) and with no minimum at 
the position H». Assuming complete Paschen- 
Back effect, the positions of the minima are given 
by the following simple formulae: 


Hot+(3Sp+H’); Hot (3Sp—H’); Hor6Sp, 


where Sp is related to H” by Sp=}H”. The 
values obtained for these quantities were 
H’=27.2 gauss and Sp=6.8 gauss. 

These two constants, H’ and Sp, are sufficient 
to account for the positions of all six lines and 
for the variations in their position with gross 
changes of field. Furthermore the constant H” 
can be calculated from the known value of the 
proton moment and the known internuclear dis- 
tance of H». This value agrees with the experi- 
mental value of H’”’. The consistency and reason- 
ableness of the analysis of the He spectrum shows 
that we have not omitted any important features 
of the molecular constitution in setting down the 
fundamental equation on which our analysis 
rests. 


THE D, SPECTRUM 


Turning now to the discussion of the Dz mole- 
cule we take note of the following facts. The D, 
molecule is homonuclear, the nuclei have spins 
of one and obey the Bose statistics. The lowest 
rotational state of D: is therefore a state in which 
the wave function of the nuclear spins is sym- 
metrical and the total spin angular momentum J, 
is either zero or two. For the next higher rota- 
tional state, J=1, (para-D:) the wave function 
for the two nuclear spins must be antisymmetrical 
and the total spin angular momentum is 1. For 
the state J=2 the wave function is again sym- 


| 

els are 

ur ex. 

lature, | 

> Mag- 

1€ two | 

1from 

is not 

of the 

uteron 

y large 

erfine 

small. 

te the 

which 

uld be 

y level 

in the 

uclear 

d pro-. 

‘be of 

mole- 

these 

olving 

en the 

quad- 

ecause 

which 

of the 

as in 

oy the 

‘thod.* 

to 

ar” to 

ved in 

fo and 

dasa 

bsolute 

| be re- 

charias, 

Jr. and 


680 KELLOGG, RABI, 


metrical and J is either zero or two. At liquid- 
nitrogen temperatures almost all the molecules 
are in these states and they are therefore the 
only ones we need consider. (See Table I of 
paper I). 

This may most easily be done under the classi- 
fication of their total rotational angular mo- 
mentum and spin: 


1. J=0, J=0. 9.3 percent of the beam. These 
molecules are inert in our apparatus since they 


have no magnetic moment. 
2. J=0, I=2. 46.7 percent of the beam. These 


molecules give rise to the deep central minimum 
which is so prominent on any of the Dz nuclear 
curves and which was used in paper J for the 
determination of the deuteron magnetic moment. 

3. J=2, I=2. 8.7 percent of the beam. These 
molecules are divided among (2J+1)(27+1) =25 
energy states in a strong magnetic field. Thus 
each state contains 0.35 percent of the molecules. 
The selection rule Am;= +1, for Am;=0 gives, 
for each my, level, 4 possible transitions. Thus 
there should be a total of 20 transitions possible. 
Since in our experiments transitions from level 5 
to level a are counted as well as those from a 
to b, there would be 0.7 percent of the total beam 
in each of these 20 minima. However the effect 
of the velocity distribution operates to decrease 
this intensity by about 25 percent, and since, 
using a D, beam our apparatus will not detect 
changes of beam intensity of less than } percent 
of the total beam, we did not expect to observe, 
and did not observe any minima from molecules 
in this state. 

4. J=2, I=0. 1.7 percent of the beam. The 
same considerations apply here as far as intensity 
is concerned as to the molecules in state J=2, 
I=2. Furthermore the only transitions possible 
are those for which Am, = +1 and this would give 
rise to lines in the rotational spectrum. 

5. J=1, J=1. 33 percent of the beam. These 
are the molecules with which we are concerned 
in the experiments on D, reported in this paper. 
In a strong magnetic field there are 2/+1=3 
possible orientations of J and for each of these 
2I+1=3 possible orientations of J; altogether 9 
states. Each energy state should therefore con- 
tain 3.7 percent of the molecules of the beam. 
Between these 9 levels there should be 6 possible 
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transitions corresponding to changes of m, o 
+4 while m, remains unchanged. If the effegs 
of the velocity distribution is neglected the ; ia: 
tensity of each of these minima should therefore 
be 7.4 percent of the total beam intensity. If the 
effect of the velocity distribution is included this 
intensity should drop to about 5.4 percent of the 
total beam intensity. ; 

Thus, excepting the intensity consideratj 
the situation for para-D2 molecules in the 
J=1, I=1 is exactly that which obtained fg 
ortho-H: molecules as discussed in paper I. Be. 
cause of our knowledge of the H: constants, we 
expected to be able to predict the positions of the 
six Dz minima at least approximately. This 
should be made clearer by the following dig. 
cussion. 


PRELIMINARY DiscussION OF De SpEctruy 


In« the section called theory we will derive 
expressions for the nine energy levels (given jp 
Table I) of para-Dz on the assumption that there 
are three types of interaction to consider, ie, 
spin-orbit, with constant H’; spin-spin, with 
constant H”; and nuclear electric quadrupole, 
with constant H’”’. The possible transitions be. 
tween these energy levels are given in Table I]. 
Inspection of this table shows that if the electric 
quadrupole interaction constant, H””, is set equal 
to zero, the transitions between the energy levels 
for Dez will be represented by exactly the same 
expressions as for He. The values of the constants 


Taste I. Energy levels correct to the third-order perturba 
tions for the first rotational state of para-D, molecules ing 
pS tic - The terms in H’ and Sp give the first-order per 

ion, the terms in C2 and C,' thé second, and the termsia 
C, the third. The _ of the quantities are as follows: 
=0.515; 
Cy’ = (H’—3Sp)*/(1—a) 
C,=(H 
H" =(up/r*)m; =(—5eqQ/4up). 
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Hl and Sp will of course be different, but these 
should be readily calculable from the known 
yalues of the constants for Hs. The spin-orbit 
constant, H’, can be taken to be one-half of that 
for hydrogen since the angular velocity of rota- 
tion of the Ds molecule is half as great and the 
internuclear distances are so nearly the same for 
the two molecules. We can make this assumption 
with some confidence because it has been shown 
experimentally" that the ratio of the rotational 

etic moments of the two molecules is 
2 to 1. Similarly the value of the spin-spin 
constant H’’ = (up/r*)» can be obtained for Dz 
since (up/uP) and up/r* are both known. These 
constants then become H’=13.6 gauss and 
H" =10.5 gauss. Sp=H’’/5=2.1 gauss. Inserting 
these values into the expressions given above for 
H, we obtain the predicted spectrum with minima 
as indicated by vertical lines in Fig. 1. 

The complete lack of correspondence between 
the experimental and predicted curves is some- 
what startling, especially so when one recalls how 
well the theory fits the case of He. As is evident 
from the figure it was expected that the six-line 
D; spectrum would fall within the wings of the 
deep central minimum arising from the state 
J=0, since the ratio of the depths of the small 
to the large minima is 1 to 7. 

It was suspected that the discrepancy between 
theory and experiment was due to a quadrupole 
moment of the deuteron. 

The apparatus was remodeled to increase the 
resolving power for a more accurate study of the 
phenomenon. It was shown in paper I that the 
resolving power of our apparatus depends upon 
the length of time that the molecule remains 
under the influence of the radiofrequency field. 
In fact, the minimal half-width (width, AH, in 


TaBLe II. Changes of energy associated with the change of 


orientation of the nuclear spin by +1 for -D 
the frst roleiional He) 
my SE/up 

1 0—1 H+H’'—3Sp+0 +C,'+0 

1 -1—0 C3 

0 0—1 H+0 +6Sp+0 +C,'+2C; 

H+0 —6Sp+0 +C,'—2C; 
-1 0—1 C; 
H—H'+3Sp+0 +C,'+0 


"N. F, Ramsey, Jr., Phys. Rev. 55, 595 (1939). 
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Fic. 1. Radiofrequency spectrum of D; at 80°K for 
changes of total spin magnetic quantum number m, of 
+1 with my, unchanged. This curve was taken at a time 
when the length of the radiofrequency field was only 2.7 
cm. The spectrum icted on the assumption of no 
quadrupole moment is indicated by the dotted lines, both 
as to _o- and relative intensity. The theoretical 
minimal half-width for any of these minima, 30 gauss, is 
in good agreement with the observed curve. 


gauss at half-depth) for a transition without any 
fine structure is given by AH=hIv/yL where L 
is the length of the radiofrequency field and 0 is 
the velocity of the molecule traversing the field. 
At the time the experimental curve of Fig. 1 was 
taken the path of the molecules in this field was 
2.7 cm. We calculate for this length of field and 
for Dz at 80°K the minimal half-width of a line 
as predicted by the above equation to be about 
30 gauss. By increasing the length of the radio- 
frequency field to 13.5 cm this minimal half- 
width can be cut down by a factor of 5. 

The apparatus was accordingly lengthened. 
The deflecting magnets were separated enough 
to allow for a homogeneous field 15 cm long. This 
magnet and its radiofrequency field wires were 
also used in the experiments of paper I. 

The curve of Fig. 2 represents data obtained 
with the improved conditions. The subsidiary 
minima are resolved from each other and from 
the central minimum. The positions of the 
secondary minima are in violent disagreement 
with the provisional theory outlined above. In 
fact the interval between the extreme lines is six 
times greater than was expected. 


DetTAILED DiscussION OF D, SPECTRUM 


A set of energy levels which accounts com- 
pletely for these lines will be developed in the 
section called ‘‘Theory of energy levels” of De. 
The only assumption in addition to those made 
for He is that a nuclear electric quadrupole inter- 
acts with the molecular electric field. The nine 
energy levels for Dz in the first rotational state 
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Fie. 2. Rote spectrum of D» at 80°K for Am;= +1, with m, unchanged. The length of the 


radiofrequency fiel 


is 13.5 cm. Dotted lines indicate spectrum predicted on the assumption of no quadru- 


le moment. Arrows indicate positions of lines predicted from the complete theory and the constants 
rom the 2.419-MC data. Theoretical minimal half-width for any of these minima is about 6 gauss. The 
central minimum is not symmetrical about the central dotted line because of the method of producin 
the radiofrequency field, as discussed by S. Millman (Phys. Rev. 55, 628 (1939)). If the directions of all 
of the static magnetic fields are reversed, the asymmetry is reversed. This result was used to show that 


the deuteron magnetic moment is positive. 


are given in Table I which is complete up through 
the third order of the perturbation theory ex- 
panded about the strong field levels. The terms 
in H’ and Sp to the first power are the first-order 
perturbations produced by the three types of 
interaction mentioned above. These terms are 
independent of field and remain even in the pres- 
ence of a large applied field. The terms C, and C,’ 
are second-order terms which are inversely pro- 
portional to H and therefore vanish in high fields. 
They arise because of incomplete Paschen-Back 
effect or incomplete decoupling of the total 
nuclear spin and the molecular rotation. The 
same is true of the third-order terms C; except 
that the decrease with increasing field is even 
more rapid. Throughout the following discussion 
the terms in C3 can and will be neglected. 
Energy differences for transitions of m; of +1 
with my, held constant are listed in Table II. 
Dividing the terms of Table II through by up 
and setting upH»=hf, we obtain expressions for 
the positions of the lines for the frequency fo as 
given in Table III. 
As in paper I we can use the asymmetry of the 


positions of the minima to identify the individual 
lines with the quantum numbers of the transi- 
tions given in Table III. This identification can 
be made unambiguously from the 2.419-MC data 
of Table IV, since, for the fields at which these 
transitions occur, all perturbations higher than 
second order have no appreciable effect, and the 
only terms which render the pattern unsym- 
metrical are C; and C,’. It can be seen from 
Table III that four of the lines should be shifted 
from symmetrical positions about H» by amounts 
— C2’ while two other lines are shifted by amounts 
—(C2—C,'). Thus the midpoint between one 
pair of lines should be shifted from the position 
Hy by —(C2—C2’) while the midpoints of 


two other pairs should be shifted from Hp by 


amounts —C,’. 
This is just the type of shift which is found 


experimentally as can most easily be seen from 
the 2.419-MC data as follows: 


3(Cr+Cz) —Ho= —0.9 gauss, 
3(Br+B 1) —Ho= —6.0 gauss, 
3(Ar+A1)—Hy=—1.1 gauss. 
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We see from this that the B minima correspond 
to the second and fifth rows of Table IV, and 
therefore that (C2’—C2) is about 6 gauss. Simi- 
larly since the A and C minima correspond to the 
other rows we conclude that C,’ is approximately 
1 gauss. Substituting this value and H=3700 
gauss into the expression for C;’ we can arrive at 
a preliminary rough value for (3Sp—H’) of +43 
gauss. Since the A minima are displaced from Ho 
by very nearly that amount and the C minima 
by much more we can take the interval between 
the A minima of 90.3 gauss as 2(3Sp—H’). Thus 
the A minima correspond to the third and fourth 
~ rows of Table III. This leaves the C minima to 
go with the first and sixth rows. We can now 


summarize as follows: 


2(3Sp —H’) =Ar—Axt=90.3 gauss, 
2(3Sp +H’) = Br—Bzi=147.0 gauss, 
2(6Sp) = Ce— Ci = 234.2 gauss. 


In the above summary we have taken Sp and 
H' as positive for the following reason. Since 
(3Sp—H’) is numerically less than (3Sp+H’) we 
conclude that Sp and H’ must have the same sign. 
Now we know that H7’ is positive because it is the 
spin orbit constant for Dz and would have the 
same sign as the corresponding spin-orbit con- 
stant for He, which was shown in paper I to be 
positive. From this consideration we will later 
derive the sign of gQ. This completes the identi- 
fication and the detailed assignment of quantum 
numbers as given in Table III. 

In order to evaluate the constants H/’ and 3Sp 
from the data as precisely as possible the follow- 
ing procedure was adopted. First a preliminary 
evaluation was made from the equations above. 
These values were used to evaluate the constants 

‘Cx, C2’ and C3 as given in Table V. With these 
constants and the 2.419-MC data of Table IV a 
TABLE III, Magnetic fields for which resonances will occur 


for changes of orientation of the total nuclear spin by +1 for 
para-D, molecules in the first rotational state. Ho=hfo/up. 


MAGNETIC FIELD IN GAUSS 


Ho— H'+3Sp—C,'+0 +0 
Ho— Cs 
Hot+0 —6Sp—C2'+0 —2C; 
Hot+0 +6Sp—C.'+0 +2C; 
Cs 
Ho+H'—3Sp—C,'—0 +0 
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set of six equations linear in H’ and 3Sp can be 
set up and solved by the method of least squares. 
The values obtained this way are 


H’=14.00 gauss and 3Sp=58.86 gauss. 


The general consistency of the experiment with 
the theory is best exhibited in Table V. The two 
experimental parameters H’ and 3Sp are used to 
predict the positions of the six minima. The 
greatest deviation between the calculated and 
experimental values is 0.5 gauss. This is better 
than 1 gauss, the estimated probable error in the 
determination of a minimum. The second-order 
perturbations as evaluated experimentally above 
are 1 gauss for the A and C minima and 6 gauss 
for the B minima. These are to be compared with 
the theoretical values of 1.1 gauss, and 5.6 and 
5.9 gauss given in Table V. 

Further agreement of experiment with the 
theory is provided by the 1.300-MC data as ex- 
hibited in Table VI. The constants H’ and 3Sp 
as evaluated above from the 2.419-MC data have 
been used for predicting the positions of the lines 
for the frequency of 1.300 MC. Table VI shows 
that for these low fields (approximately 2000 
gauss) the perturbation theory used does not 
converge. For the minima Cr and Cy, the third 
order is as large as the second. It was therefore 
necessary to calculate the energy levels directly 
from the secular equation and to find the posi- 
tions of the minima from them. For four of the 
lines the agreement between predicted position 
and experimental position is unusually good, 
whereas for the other two lines the agreement is 
just within the sums of the limits of error for the 
ments are considered reliable to +1 gauss for the 2.419-MC 
data and to +1.5 gauss for the 1.300-MC . Note that By 
and Br are much more unsymmetrically placed than the 


others ; that they are all more symmetrical in the higher field, 
Ho= 3710 gauss. 


DISPLACEMENTS IN GAUSS FOR 
OSCILLATOR FREQUENCIES OF 


2.419 MC 1.300 MC 


— 118.0 — 118.0 
— 79.5 — 84.0 
— 46.3 — 47.6 

(3710) (1995) 
44.0 40.7 

67.5 61.3 
116.2 112.0 
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TABLE V. shows of the calculated 
perturbations of the lines in the D in Dz spectrum at 2.419 MC 
and the comparison with the experimentally found perturba- 
tions. The ‘‘Line Designations” are those assigned to the 
minima of Fig. 2. 


PERTURBATIONS LINE DESIG- 
NATION AND 
First SECOND EXPERIMENTAL 
ORDER ORDER ORDER TOTAL |PERTURBATION 
6Sp -—4140/H —7.3X10*/H? Cr 
117.7 —1.1 —0.4 116.2 116.2 
3Sp+H’ —21,300/H  3.6X10*/H? Br 
72.9 —5.6 0.2 67.5 67.5 
3Sp—H’ —4140/H 0 AR 
44.9 —1.1 0 43.8 44.0 
-—(3Sp-H’) -—4140/H 0 AL 
—44.9 -1.1~ 0 —46.0 —46.3 
—(3Sp+H’) -—21,300/H —3.6X10*/H? BL 
: -5.9 —0.2 -79.0 —79.5 
-6SD —4140/H 7.3 X108/H? L 
—117.7 —1.2 0.4 —118.5 —118.0 


constants and the experimental determination of 
the minima. This is most likely due to the fact 
that the 1.300-MC data were taken only once, 
whereas the 2.419-MC data are the result of 
averaging several runs. 

The quantity Sp is defined as }(H’’+H’”) so 
that the experiment with D, is incapable of pro- 


viding a determination of either H” or H’” 


separately. We can, however, take the value of 
H” from the experiment on He: as discussed 
above. That value of H” of 10.5 gauss subtracted 
from 5Sp=98.1 gauss yields 87.6 gauss for H’”’. 
It is to be noted that this value of H’” is positive. 
The significance of this quantity will be discussed 
later. 


EXPERIMENTs ON HD 


Experiments on HD were undertaken to prove 
that the quantity H’” is certainly a property of 
the deuterium nucleus and is not just a peculiar- 
ity of the Dz molecule. This proof takes the 
following simple form. In HD one set of transi- 
tions can occur for changes of the magnetic 
quantum number for the proton and another set 
for the deuteron. In the first case the pattern 
should be completely accounted for without any 
necessity for considering a quadrupole moment 
‘of the deuteron because even if such a quantity 
does exist the changes of energy of the molecule 
accompanying changes of mp would not be 
affected thereby. This follows, since if no re- 
orientation of the deuteron occurs within the 


molecule relative to the extranuclear electricay 
configuration, no changes of energy would Occur, 
On the other hand, if we consider the transitions 
of mp, then such energy changes should occur 
and the spectrum should show the Presence of 
the quadrupole moment. As this turns out to be 
the case these experiments show that the effect 
is a property of the deuteron and not some 
molecular property. 

Furthermore in HD the spin-spin and quadry. 
pole interactions can be evaluated experimentally 
and therefore the value of the deuteron quadrupole 
moment will not depend on subtracting a calculated 
quantity, H’’, as in De. A further check on the 
correctness of the theory will be afforded by a 
comparison of the interaction constants obtained 


with HD with those previously obtained with 


He and Ds. 


NUCLEAR SPECTRA OF HD 


The HD molecule is heteronuclear and there 
is no quantum number of total spin. There are 
therefore three coupled angular momenta J, J, 
and Ip with projections for the strong field case 


of m;, mp, and mp. At the temperature of liquid’ 


Ne, 63 percent of the molecules are in the state 
J=0, and for this state all interactions of the 
type discussed in this paper average out to zero, 
These molecules are therefore ideally suited to 
the determination of the nuclear magnetic mo 
ments of the proton and the deuteron but are of 
no use at all in the present experiment except for 


the location of Ho=hAfo/u. The remaining 37 


percent of the molecules are in the state J=1, 
only 0.3 percent being in the state J=2. There- 
fore any fine structure that is observed is im- 


mediately attributable to the molecules in the 


TABLE VI. Comparison of experimental D in perturba. 
tions at 1.300 MC with the perturbations as calculated with the 
constants obtained from the 2.419-MC data of D in Dy. 


A 

DESIG- First SECOND THIRD MENTAL 

NATION ORDER ORDER ORDER SECULAR VALUE 
Cr 117.7 -—1.96 —1.63 113.0 112.0 
Br 72.9 —10.36 86 63.7 613 
Ap 449 —2.04 0.00 43.0 40.7 
AL —-449 —2.12 0.00 —474 —476 
By -—72.9 —11.11 -—84.6 — 84.0 

L —117.7. —2.20 2.07 —117.7 | —1180 
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Fic. 3. Radiofrequency spectrum of HD at 80°K for Amp= +1 with mp and my unchanged, 


state J=1. In a strong magnetic field there are 
three orientations of J, three orientations of the 
deuteron spin, and two orientations of the proton 
spin. This makes eighteen possible states which 
in strong field gives eighteen different energy 
levels. There are three types of transitions among 
these energy levels, i.e., (1) those for which mp 
changes by +1 giving twelve different transi- 
tions; (2) those for which mp changes by +1 
giving rise to nine different transitions; and (3) 
those for which m, changes by +1 giving twelve 
different transitions. The first group (called D 
in HD) clusters around Hy=hfo/up and will be 
discussed first. The second group (called H in 
HD) clusters around Ho=hfo/2up and will be 
discussed next. The third group which clusters 
around Ho=hfo/ue will be discussed in a separate 
paper devoted to the general question of rota- 
tional magnetic moments. 


D w HD 


_ Figure 3 shows the group of minima associated 
with changes of mp of +1 at the frequency of 
2.419 MC. The deep central minimum is not 
completely shown in the graph so that the small 
minima for J=1 will show to better advantage. 
Besides the central minimum there are eleven 
clearly resolved minima one of which is taken to 
be double since twelve are expected altogether. 


The experimental procedure for taking this 
curve is similar to that for the other data except 
that for this case particular attention has to be 
paid to the magnitude of the deflecting and 
refocusing fields. In the HD molecule the de- 
flections are chiefly due to the large proton 
moment so that if a very large deflecting field is 
used to affect the motion of the molecule by 
action on the smaller deuteron moment the ex- 
cursions produced by the proton moment are so 
large that the refocusing conditions are not good. 
Therefore a compromise is sought before the data 
is taken for one curve. 

The column labeled ‘‘Experimental”’ in Table 
IX represents the average of two sets of data for 
D in HD at the one frequency. 

In the section called ‘‘Theory of energy levels 
of HD” we will derive expressions for the eighteen 
energy levels of HD in the first rotational state 
on the assumption that there are three types of 
interaction to consider, i.e., spin-orbit interaction 
of the deuteron with the magnetic field of the 
molecular rotation at the position of the deuteron 
and the similar type for the proton, with inter- 
action constants upH'’p and 2upH'p; spin-spin 
interaction of the deuteron with the proton with 
interaction constant S=(2upup/Sr*)w; and the 
interaction of the deuteron quadrupole moment 
with the gradient of the electric field with inter- 
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TABLE VII. Energy levels for HD molecules in the first rotational state. The energies 
VR, MP, and up are given in nuclear magnetons. This notation is used because the econ 


nu ly using experimental values for the constants H'p, H'p, S, and E. S=(2upup/5r*)w. E= =2Q/4yp, 


KELLOGG, RABI, RAMSEY, JR. AND ZACHARIAS 


are in gauss nuclear magnetons if 


and third-order expressions are given 


ENERGY IN GAUSS NUCLEAR MAGNETONS 
ms mD ZEROTH ORDER First SECOND Tuirp 
1 4 1 | —uoH —upH'p + S — 0 +0 
1 4 O|} —uR —yupH'p +10,160/H +1.75K 100m 
1 —1| +upH —yrH —ypH'p +unH'n — S — +30,660/H 
1 +upH'p — S — +621/H +0 
i- 0 | +up upH'p +2uvE +34,400/H — 4.75 108) 
1 +ueH +uoH —yrH +urH'p + S — +23,600/H —8.24 x 
0 | 1 | —poH —2S +2upE — 10,780/H — 1.75 Xx 
0 0 | —ueH —4upE —11,110/H +2.14X 
0 —1 +upH +2S +2upE +30,400/H +4.72X106/m 
1 | —uoH +2S +2unE — 30,400/H 
O| +urH —4u +11,110/H +2.14X 
0 +ueH +2 —2S +2yupE +10,780/H —1.75 
1 | —yoH +urH +upH'p + S — — 23,600/H — 
=i 0 | —ueH +ur +upH'p +2upE —34,400/H —4.75 
—1 | +urH +urpH'p —uvH'n — S — —621/H +0 
+prH —peH'p +uoH'p — S — — 30,660/H +6.34X 
-1 0 | —upH'p +2uvE — 10,160/H +1.75 
—1 +ueH +uoH —upH'p + S — +0 


TABLE VIII. Magnetic fields at which resonances will occur for changes of orientation of the deuteron spin in HD 


molecules in the first rotational state. Ho=hfo/up. S’ =(2up/5r*)m 


MAGNETIC FIELDS AT WHICH RESONANCES OCCUR 

my mp Amp FIRST ORDER SECOND ORDER THIRD ORDER 
1 O<>1| Ho —H'p + S'-3E —11,900/H — 2.05 

1 } Ho —H'p + S'+3E — 24,(00/H — 5.38 X 

O<1| Ho —H'p — S'—-3E — 39,500/H +5.56X 

1 Ho —H'p — S'+3E +12,600/H +4.09 
0 O<>1| Ho —25S'+6E +386/H —4.55 

0 —1<0] Ho —2S'—6E —48,500/H — 3.02 

O<>1 Ho +25S'+6E — 48,500/H +3.02 X 

0 Ho +2S’—6E +386/H +4.55 
O<1)| Ho +H'p + +12,600/H — 4.09 
| Ho +H'p + S'+3E —39,500/H —5.56X10/R 
—§ O<>1 | Ho +H'p — S'-3E — 24,000/H +5.38 X10/B 
A +H'p — S'+3E — 11,900/H +2.05 X10/R 


TaBLe IX. Collected data on D in HD. Frequency 2.419 MC. 


PERTURBATIONS 
PREDICTED FROM Hz AND D2 CONSTANTS 
EXpeERI- LINE BATIONS 
First SECOND THIRD MENTAL DESIG- From HD 
my mp mD ORDER ORDER ORDER TOTAL VALUE NATION CONSTANTS 
1 } 0<1 — 59.11 —3.26 —0.14 — 62.5 — 62.2 CL — 62.6 
1 | —1le0 +45.59 — 6.39 —0.39 +38.8 +39.0 Br +390 
1 — 86.39 — 10.89 +0.40 —96.9 —98.7 —9%6.3 
1 —} -1le0 +18.31 +3.38 +0.30 +22.0 +23.2 Ar +228 
0 7 O<l +77.42 +0.1 —0.33 +77.2 75.9 Dr +780 
0 -—1<0 — 131.98 — 13.55 —0.22 — 145.8 — 145.9 Fr 145.4 
0 -—}3 0<1 +131.98 — 12.62 +0.22 +119.6 +118.8 Fer 1193 
0 —77.42 +0.1 +0.33 —77.0 —77.7 Di | 
-1 0<1 — 18.31 +3.41 —0.30 —15.2 —15.0 AL —160 
—-1 | —led +86.39 — 10.40 —0.40 +75.6 +75.9 Er +751 
-1 - —45.59 —6.55 +0.39 —51.8 —51.5 Br -519 
+59.11 —3.15 +0.14 +56.1 +55.8 Cr +562 
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action constant = e*gQ/4. These energy levels 
are given in Table VII which is complete up 
through the third order of the perturbation 
theory expanded about the strong field levels. 
Algebraic expressions for the second- and third- 
order terms are not given in the tables as the 
expressions are too clumsy to print. Since these 
expressions involve the constants H’p, H’p, S, 
and E, the numerical values given in the table 
were calculated after the constants were known. 
For our purposes they are most conveniently 
given in the peculiar energy units gauss-nuclear- 
magnetons. They can be converted to ergs simply 
by multiplying by 5.04X10-** erg/gauss/n. m. 
Just as in Ds, these terms vanish in very high 
field when the Paschen-Back effect is complete. 
However, for HD this field must be unusually 
large to make the magnetic energy of the rota- 
tional moment different from that of the deu- 
teron moment by an amount large compared 
with the other interactions. This is true because 
the rotational moment (0.660 n. m.) is so nearly 
equal to the deuteron moment (0.855 n. m.) 

Table VIII gives expressions for the positions 
of the lines accompanying changes of mp by +1. 
This table is derived from Table VII just as was 
done for the corresponding table (Table III) 
for De. 

The simplest procedure to follow in the analysis 
of this spectrum is to estimate the values of the 
constants in the theory from the experimental 
values of the constants obtained with Hz and D, 
and then to use these constants to predict the 
positions of the spectral lines. If this procedure 
proves successful in predicting the observed lines 
we will automatically have an identification of 
the lines with the quantum numbers of the transi- 
tions between energy levels. Once such an identi- 
fication has been made we can reverse the process 
and calculate the values of the constants from 
the experimentally determined line positions. 
This direct intercomparison of the constants for 
the three different molecules is possible because 
the internuclear distances are almost identical. 
In fact, the values of (r-*), for He and HD differ 
by only 1 part in 400. This enables us to make 
an excellent estimate of S’ = (2up/5r*), = 13.64 
gauss. Since the electron densities of these mole- 
cules depend primarily on the nuclear charge 
and the internuclear distance, the gradients of 
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the electric field with which the quadrupole 
moment of the deuteron interacts will be equal 
in Dz and HD. The interaction constant E can 
be taken over directly from H’” of De. Similar 
considerations for the constants H’p and H’»p are 
not as certain for two reasons. First, there is no 
complete quantitative theory for the correspond- 
ing constants in Hz and Dez which could be 
modified for HD. And second a naive theory 
suffers because of the asymmetry of the HD 
molecule. For instance, it is not immediately 
evident that H’p should equal H’p and this 
question does not arise for Hz or De. However, 
we know that the rotational magnetic moments 
of He, HD, and Dz are 0.88, 0.66 and 0.44 nuclear 
magneton, respectively"' and these quantities 
are in inverse ratio to the reduced masses of these 
molecules. We also know that (H’)He is 27.2 
gauss and (//’)pe is 14.0 gauss again in inverse 
ratio to the reduced masses. We therefore assume 
that H’p =H’ p= 20.4 gauss. 

The values indicated above have been sub- 
stituted into the formulas for the first-, second-, 
and third-order perturbations for the D nuclear 
transitions of HD. The results of these calcula- 
tions are shown in Table LX. The column marked 
“Total” is simply the sum of the entries in the 
three columns to the left of it. The experimental 
values and their designations are placed so that 
they correspond most nearly to the calculated 
values. The agreement in all but two cases is 
good. The two experimental values which do not 
agree with the calculated perturbations are 
necessarily poor because they are merely the 
center of gravity of a pair of unresolved lines. 
No further argument, then, is necessary in 
justification of the assignment of the quantum 
numbers of the transitions to the experimental 
lines. 

Evaluation of the constants H’ p, S’ and E from 
the experimental data given in Table IX was 
made in the following way. The second- and third- 
order perturbations were evaluated using the 
preliminary constants from Hz and De». These 
appear in Table IX. The algebraic expressions 
for the first-order terms were equated to the ex- 
perimentally evaluated perturbations diminished 
by the calculated second- and third-order per- 
turbations. The lines Ag and A, were omitted 
from the group for reasons that are immediately 
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96 
Op 
FREQUENCY 
15.750 MC 
T,*.55 AMP. 
3675 3695 3715 


MAGNETIC FIELD IN GAUSS 


Fic. 4. Radiofrequency spectrum of HD at 80°K for Amp= +1 with mp and my, unchanged. 


evident in Fig. 3. The experimental lines Ez and 
Dr were taken together by averaging the theo- 
retical first-order formulas and equating the re- 
sulting expression to the displacement of the one 
unresolved minimum which represents them 
both. Using this latter equation twice we get ten 
simultaneous equations which were solved by the 
method of least squares for the three constants, 
and we obtain: 


H'p= 20.07 gauss; S’ = 13.348 gauss; 
E=17.495 gauss. 


The perturbations calculated with these con- 


stants are given in Table IX: Agreement with the 
experimental values is naturally better than with 
the preliminary values. In fact, except for those 
minima (Apr, Ax, Dr, Er) which have already 
been discussed, the experiment can be regarded 
as internally consistent to better than 0.5 gauss. 


H HD 


Figure 4 shows the group of minima for changes 
of mp of +1 at a frequency of 15.750 MC. For 
these data particular attention has to be paid to 
the constancy of the homogeneous magnetic 
fields, while the intensity measurements are 
made. This is difficult because the intervals 
between the minima are a very small fraction of 
the total field. In fact, for the data to be of real 
quantitative value comparable with the other 


results of this set of experiments it was 
to hold the field constant to 0.2 gauss in 3709 
gauss. This implies that the field excitation 
current had to be adjusted and read to 0.01 amp, 
in 200 amp. Further the slight hysteresis in the 
electromagnet makes it necessary to vary the 
current in only one direction during a set of 
measurements. Because of the slowness of the 
action of the Pirani gauge detector (15 sec. per 
deflection and nine deflections per point) at least 
five hours is required to take one experimental 
curve. These requirements taxed the rheostat 
control and storage battery bank to their limit. 
It can be seen in Table X why it was not ex- 
pedient to take these measurements at much 
lower field values where these experimental 
difficulties would be reduced proportionately. 
Control and measurement of the radiofrequency 
were not difficult because of the stability of the 
r.f. oscillator and the ease of using the General 
Radio No. 620A wave meter. 

Reference to Fig. 4 shows that besides the 


deep central minimum there are eight small. 


minima for J=1, whereas it is stated above that 
there should be nine. It will be shown presently 


that the theory predicts that one of the small 


ones should be in the position occupied by the 
deep central minimum, and therefore completely 


hidden. 
The procedure followed in the analysis of the 


Hoss = 


2 


OOO 
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H nuclear spectrum of HD is entirely similar to 
the procedure previously outlined for the analysis 
of the spectrum of D in HD. Table X, which is 
obtained from Table VII and the relationship 
hfo= 2urHo, gives the expressions for the posi- 
tions of the H in HD lines. The constants in- 
volved here are H’p and S” =}(up/r*)w. A pre- 
liminary value of H’p is found by taking H’p 
=H'p=20.4 gauss as in the case of D in HD. 
The value of 5S’”’ may be estimated from 
(up/r*) = (uo 10.45 gauss. These con- 
stants substituted in the expressions given in 
Table X yield the predicted perturbations listed 
in Table XI. It is there evident that the “Total” 
predicted perturbation is in excellent agreement 
with the experimental value. The assignment of 
the quantum numbers to the experimental lines 
is therefore accomplished. 

The line expected at — 1.98 gauss is completely 
hidden by the central minimum. The other eight 
minima are clearly resolved and can be used to 
form eight simultaneous equations for evaluating 
the two constants H’p and S” (see Table X). 
The values obtained by the method of least 


squares are H’p=20.48 gauss and S’”=2.12 
gauss. Using these values to recalculate the 
positions of the minima we get the last column 
of Table XI. The internal consistency of this 
set of data is well within the limit of 0.2 gauss. 


THEORY OF ENERGY LEVELS 


D, 


The energy levels of the Dz molecule in the 
state J=1 and situated in an external magnetic 
field H, taken in the direction of the z axis, are 
given by the characteristic values of the operator 


KH = —upH' (is -J 
+ (up?/r*) [is -i2 — 3 (i: (ig 2) 

25(2J—1)i(2i—1) 

(2) 


where i, and i, are the spin operators for the two 

deuteron spins, r is the internuclear distance, and 

the other symbols have their usual meanings. 
The first and second terms represent the effect 


J)? 


X. Magnetic $ for which resonances will occur for changes of orientation of the proton spin in HD molecules 


MAGNETIC FIELDS FOR WHICH RESONANCES OCCUR 

my mp mD FIRST ORDER SECOND ORDER THIRD ORDER 

1 -je 1 Ho —H'p +25" —109/H 

1 0 Ho —H'p —4270/H +1.15 x 10*/H? 

1 -1 Ho —H'p —2S" +1250/H +2.57 X 10°/H? 

0 —te} 1 Ho +3460/H — 1.15 10*/H® 

0 0 Ho — 3880/H 

0 Ho +458" +3460/H +1.15 x 10*/H? 
-1 -te 1 Ho +H'p +1250/H — 2.57 XK 10*/H® 
-1 0 Ho +H'p —4270/H — 1.15 10*/H? 
-1 -1 Ho +H'p —2S" —109/H 

TABLE XI. Collected data on H in HD. Frequency 15.750 MC. 
PERTURBATIONS 
PREDICTED FROM Hz AND Ds CONSTANTS PERTUR- 
EXpeEri- LINE BATIONS 
First SECOND Turd MENTAL DESIG- From HD 

my mp mD ORDER ORDER ORDER TorTaL VALUE NATION CONSTANTS 

1 —je} 1 — 16.22 —0.03 0.00 — 16.25 — 16.09 CL — 16.27 

1 —j<e} 0 — 20.40 —1.15 0.08 — 21.47 — 21.64 Di — 21.55 

1 —}e>4 —-1 — 24.58 0.34 0.18 — 24.06 — 24.11 Ex — 24.20 

0 1 — 8.36 0.93 —0.08 —7.51 —7.57 Bi 7.63 
'0 -be 0 0.00 — 1.98 0.00 — 1.98 A 

0 —he} -1 8.36 0.93 0.08 9.37 9.55 Br 9.49 
-1 -le 1 24.58 0.34 —0.18 24.74 24.95 Er 24.88 
-1 -he 0 20.40 —1.15 —0.08 19.17 19.22 Dr 19.25 
-1 - 16.22 —0.03 0.00 16.19 16.37 Cr 16.21 
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of the external field on the deuteron and rota- 
tional magnetic moments, respectively. The third 
term is the “spin orbit” interaction of the deu- 
teron magnetic moments with the molecular rota- 
tion, and the coupling constant H’ has the signifi- 
cance of the magnetic field produced by the 
rotation at the location of the nuclei. The fourth 
term is the familiar dipole-dipole interaction be- 
tween the two nuclear magnets. 

The last term represents the effect of the 
nuclear quadrupole moment and is summed over 
both nuclei. The quantity g has a meaning 
similar to that of Eq. (1). 

In the appendix it will be shown that for 
identical nuclei the fourth term may be written as 


(27 


(3) 


where J is the total spin for both deuterons, 
which is 1 in our case, and 7, the spin of the 
individual nucleus, is also 1. The factor in curly 
brackets is therefore 2. It will also be shown that 
for the case of total spin J=1, the summation in 
Eq. (2) is equal to the negative of the operator 
in square brackets in formula 3. Thus the total 
quadrupole moment effect is as if we had one 
nucleus with a nuclear spin of 1 and a quadrupole 
moment of sign opposite to that of the quadrupole 
moment of the actual particles.* 

Putting in the values of J and J we may 
therefore write Eq. (2) as 


KH = —pemsH —ppH'l-J 


eqQ 
2 


(4) 


From this form it appears that a quadrupole 
moment effect cannot be distinguished from 
dipole-dipole interaction between the two nuclei. 

For purposes of convenience we will write 
as (2up/5)(H”+H’”’), where 
H" =yp/r*, and H’” = —5eqQ/4up. H” and H’” 


128 We are indebted to Dr. Nordsieck for bringing this to 
our attention. 


may be expressed in gauss and are a measure of 
the magniture of their respective interactions 

In an infinitely strong field (complete Paschen. 
Back) the diagonal elements of 3¢ give the value 
of the energy levels of the molecule. These diag. 
onal elements are given by the formula (see 
appendix) 


= — wpm; H — prmsH —ypH'mym, 


+ +H") 
(5) 


when we set J=1 and J=1. 

If H is not infinitely great these formulae hold 
only approximately and must be corrected -for 
the effect of the off-diagonal matrix elements, 
The energy levels including their perturbations 
are listed in Table I. The necessary calculations 
are indicated in the appendix. 


HD 
For the HD molecule the calculation is some- 
what lengthier. The energy expression is 
KH = —2upmpH 
—2upH' J) J) 


2upup 3(i1- 1) 


(2I 
+4 J) (6) 


where 7; and 72 are the proton and the deuteron 
spins, respectively, and H’p and H’p the spin 
orbit constants for the proton and deuteron 
which need not necessarily be the same; mp, mp, 
and me are the magnetic quantum numbers of 
the proton, deuteron, and rotation, respectively. 
We may write 


[3(is- J)* 


2 
= 
(2J+3)(2J—1)- Gi D 
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TaBLE XII. Summary of results. All quantities in gauss. 


SPIN-ORBIT 


(ub/r*) CONSTANT (—Se*gQ/4upD) 


27.2 
HD 33.441 10.24 20.141. 87.541 


HD 345 10.600.2 20.4840.2 


The diagonal elements of (6) for infinitely strong 
fields are therefore 


3¢(mpmpm s) = — 2upmpH — ppmpH 
_—2upH' pmpms —unH'pmpmy 
+4upH"[3m—J(J+1) jmpmp 
— +1) +1) J. (8) 


when we use the same definition of H’” and H’” 
as in Eq. (5). The values of 72 and J are each 
equal to 1. 

These elements together with their perturba- 
tions are listed in Table VII. 


RESULTS 


Table XII is a summary of the constants ob- 
tained from the experiments on H2, Dz and HD 
in the state /= 1. The quantities have the dimen- 
sions of magnetic field and the results are given 
in gauss. Each row gives the final results obtained 
from the spectrum which arises from the in- 
dicated transition. The first row gives the results 
of paper I. The second row gives the results ob- 
tained from the 6-line spectrum of De, and the 
third and fourth rows the results of the two 
spectra in HD. The italicized quantity 10.6 in 
the second row is the value of yup/r® from the 
fourth row. This quantity is subtracted from the 
sum [yp/r’—(5e’gQ/4up)] which is measured 
experimentally to obtain the quantity in the 
fourth column. Similarly the italicized quantities 
up/r and up/r® in the third and fourth rows are 
obtained by multiplying the directly measured 
values of up/r’ and ywp/r by the accurate ratios 
up/up and up/up, respectively, as obtained from 
the previous paper. These are entered in the table 
for the purpose of comparison. 

From the quantities given in the first two 
columns one may calculate the nuclear magnetic 
moments in question. The results agree with the 


previously published” values of the proton and 
deuteron moments as is evident from the small 
variations of the numbers in the first two 
columns. 

The quadrupole moment of the deuteron is 
calculated from the results of the fourth column. 
We precede this calculation with a discussion of 
the determination of the sign of the moment. 

In the detailed discussion of the nuclear 
spectrum of Ds» the following argument was 
advanced to lead to a positive value of the sign 
of H’’”’. The De spectrum requires that H’ and Sp 
have the same sign, and, since H’ has the value 
predicted from the He spectrum, H’ is most 
likely positive. Therefore Sp and H’” are posi- 
tive. The only seeming weak point in this argu- 
ment is the sign of H’, but it is almost incon- 
ceivable that it be negative. 

However, the D nuclear spectrum of HD per- 
mits us to be certain even with regard to this 
point. Let us make the violent assumption that 
H’ for the Dz molecule is negative. Consequently 
Sp is also negative. The internal consistency of 
the D, nuclear spectrum is not affected by this 
assumption. We would therefore have Sp= 
— 98.10 gauss, and since H”’ is positive and equals 
10.5 gauss, the quantity H’’=—108.6 gauss. 
For HD we now have the quantity 6£ =6H’’’/5 
= —130.3 gauss. This value combined with 
2S’ = +26.7 gauss gives 157.0 gauss for the first- 
order term for the line with quantum numbers 
(0, —}, —1-+0). The second-order perturbation 
is positive and the other perturbations are neg- 
ligible. Therefore a line is predicted at least as 
far from the center as 157 gauss. Experimentally 
the largest observed positive displacement is 
118.8 gauss, and this line, of course, fits with the 
assumption that H’” is positive. Therefore we 
conclude that H’” is positive. 

The measured value of (—5e*gQ/4yp) is 87.5 
gauss from Table XII. In a forthcoming paper, 
one of us (N.F.R., Jr.) will present data obtained 
from the analysis of the spectra of D, and HD 
arising from the transitions my=+1. If we 
include these values in the average the quantity 
becomes 87.2+0.5 gauss. 

To calculate Q the value of g must be inserted. 
Nordsieck” has calculated the quantity q’ =[1/R® 


# A. Nordsieck, Phys. Rev. 57, 556(A) (1940). 
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pdr(3 cos? 6—1) for the deuterium mole- 
cule. In this expression R is the internuclear 
distance and the integral is taken over the elec- 
tronic wave function of the deuterium molecule; 
r is the distance of the electron from the nucleus 
and @ is the angle which r makes with the inter- 
nuclear axis. Thus q’ is $0?V/dz*. Since our q is 
defined with respect to the rotational state of 
the molecule m=J=1, it is related to q’ by 
q= —2Jq' /(2J+3) = —2q'/5. Nordsieck’s value 
of is 1.193X10** cm-', i.e., the effect of the 
positive nuclear charge is greater than the com- 
bined effect of the two electrons. Since H’”’ = 
— 5e’qQ/4up has been shown to be positive and 
since g is negative, Q the nuclear quadrupole 
moment is positive. Thus 


87.2X0.855X5.04X10-* 2 
4.802 10-29 1.193 X 10% 


and therefore Q= 2.73 X cm?. 

This value is more precise than the value pre- 
viously reported! because of the refinement of the 
experiments and the refinement of Nordsieck’s 
calculation of g. The chief source of error is in q. 


DISCUSSION 


The positive sign of the quadrupole moment of 
the deuteron shows that the charge distribution 
is that of a prolate spheroid spinning about its 
major axis. Most of the other known nuclear 
quadrupole moments also show this surprising 
property. 

Although the quadrupole moment of the 
deuteron, 2.73 10-*? is the smallest value 
yet reported for any nucleus, it does not mean 
that the deuteron is exceptional. It must be 
remembered that the optical methods which 
have been used to study other nuclei are not 
sensitive enough to detect quadrupole moments 
very much smaller than 10-* cm?. Furthermore, 
it should be noted from the definition of Q, 
eQ=e/S p(32°—r*)dr, that for a given asymmetry, 
i.e., a given value of 32*/r?, the larger nuclei 
should have larger quadrupole moments. More- 
over, more than one nuclear charge may con- 
tribute to the moment. In fact for the heavier 
nuclei some of the quadrupole moments are so 
large that this must necessarily be the case. It 


would perhaps be more revealing to consider 
Q/Z rather than Q itself for purposes of com, 
parison. 

It may be of interest to calculate the quantity 
3z*/r? which is a rough measure of the departure 
of the deuteron from spherical symmetry, We 
write For we take the a 
squared distance of the proton from the centroid 
of the deuteron. This quantity calculated from 
the deuteron wave function" is 2.38 x 10-8. The 
value of 3z*/r? is therefore 1.115. 

It is of importance to consider whether the 
effect which we ascribe to the presence of g 
quadrupole moment in the deuteron can be due 
to some other form of interaction between the 
deuteron spin and the molecular rotation. The 
experimental result which is observed is ap 
interaction of the dipole-dipole form. As J. H, 
Van Vleck has pointed out to us, such effects can 
arise in the second-order interaction of the 
nuclear spin with the electrons of the molecule. 
However, the excellent agreement of the values 
of the proton moment, as obtained directly, and 
from the analysis of the Hz spectrum in which 
the moment is evaluated from the dipole-dipole 
interaction, precludes this possibility. This point 
of view receives further support from the results 
presented in this paper, which also show that the, 
dipole-dipole interaction in HD is exactly the 
magnitude to be expected from the known values 
of the nuclear moments. To this evidence must 
be added the experimental proof that the elec. 
tronic wave functions of He, De, and HD are not 
greatly changed because of the differences in 
their moments of inertia, as is shown by the fact 
that the rotational magnetic moments" are 
closely proportional to the angular velocities 
of rotation in the state J=1. This is likewise true 
for the spin-orbit coupling which, as has been 
pointed out, may be considered as the magnetic 


field produced by the rotation of the molecule at 


the positions of the nuclei. 

As has been suggested by Schwinger," the 
presence of a quadrupole moment in the deuteron 
indicates that the ground state of the deuteron 
is not a *S, state but may be a mixture of *S; 


3H. Bethe and R. Bacher, Rev. Mod. Phys. 8, 112 
(1936) (Eqs. (44a) and (44c)). 
4 J. Schwinger, Phys. Rev. 55, 235 (1939). 
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3p,. On this view the deuteron magnetic ment of the neutron moment by Alvarez and 
moment will be a resultant of the contributions Bloch. Within the experimental errors the sum 
not only of the spin moments of the proton and of the moments of the proton and the neutron is 
neutron but also of the orbital moment of the equal to that of the deuteron. Whether this 
ton. If the moments of the proton and neutron apparent additivity conceals a real nonadditivity 
were strictly additive, the effect of the quadru-_ will not be known until an adequate theory of 
moment would result in a seeming departure the deuteron is given which accounts quanti- 


. : from this additivity in that the deuteron moment  tatively for all the known facts including the 
oa would not be just the sum of the proton and magnitude of the quadrupole moment. 


’ neutron moments. This investigation has been aided by grants 
Fortunately it is possible to make this interest- from the Research Corporation and the Carnegie 
ing comparison because of the recent measure- Institution of Washington. 


of a 
be due APPENDIX 
a The magnetic dipole-dipole interaction of Eq. (2) of the text contains the operator [(i:-is)r* 


—3(i,:1) (is) J/r*. If we introduce the total nuclear spin I=i,+i: this expression can be written as 


is an 
s J. — — 3 - ] — — 3 (9) 
Cts Can 
bs the By a general theorem"* we can write 

ecule, 
| (x-+ty)? 
ly, and r (2J —1)(2J+3) 
(x-+iy) 
dipole +i) + (10) 
s point (2J—1)(2J +3) at 
results 

=— 
tly the r (2J—1)(2F+3) 
| values 
e must When formula 9 is expanded and Egs. (10) substituted it becomes 
he elec. 
nces in (2J-1)(2J+3) 
locities If i; and ig are equal, then for all matrix elements, m;, m,;, I, J->m,', mz’, I, J, this operator is 
ise true equivalent to the operator 
i 
(12) 
(27 —1)(21+3)(2J —1)(2J+3) 

r,* the for any value of the total spin J. There are matrix elements where AJ = 2 which can be obtained from 
euteron formula 11 but not from formula (12). For Dz in the state J=1 these additional terms do not exist. 
euteron The coefficient multiplying the operator is obtained by using the expressions for the matrix ele- 
» of 4S, ments of the spins 7, and 7, when combined to give a total spin I.” 
5. 8, 112 # L. Alvarez and F. Bloch, Phys. Rev. 57, 111 (1940). 


“H. A. Bethe, Handbuch der Physik, Vol.'24, No. 1, Eq. (65.36). 
”E. U. Condon and G. H. Shortley, Theory of Atomic Spectra (Cambridge University Press, 1935), Chapter III. 
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To show the total effect of the two deuteron quadrupole moments in Eq. (2) of the text we simply 
add the two operators in the sum. From Eqs. (11) and (12) this operator becomes 
I(I+1) +1) 
(2I—1)(2I+3) 
For the special case of De, i: =i2=1, J=1, the first parenthesis of 13 is equal to —1. Thus it is seen 
that the combined effect of the two deuterons is equivalent to that of a single spin of unity with 3 


quadrupole moment of opposite sign. 
The matrix elements of 


(14) 
are easily written down from the expansion 
Utilizing the rules of matrix multiplication and the matrix elements of angular momentum, 
(m|L.|m) =m, 
(16) 


we obtain all the matrix elements of F. The first term on the right-hand side of 15 gives all the 
diagonal elements; the second and third all the elements where m,;, m;—m ;+1, m;#1; and the 
last two terms the elements ms, m;—m s+2,m,;—2 and m;—2,m,+2, respectively. These ele. 


ments are 
(mym;| —I(I+1) [3m ?—J(J+1)], 
(mymy| 
(mym;| F| my+1m;,—1) }}, (17) 
(mymz| F| msy—2m,+2) 

(I —m, 
(mymz| F|ms+2m,—2) =3[(J—mys) (J —ms—1)(J+ms+1)(J+ms+2) 
}, 


These elements are evaluated for D2 by setting J=J=1. 
In the case of HD there is no total spin angular momentum. We multiply up the expression in 


the bracket of the right-hand side of Eq. (7) and write 
G=3(is- J) J) +3 (G2: J) J) — (J +1) 
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SECONDARY EMISSION FROM SILVER 


The matrix elements are obtained with the use of the relations of Eqs. (16). The letter i without 
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subscripts is ./—1. For the particular case of HD we set i;=}, i2=1, and J=1. The spin orbit 
interactions of the form I-J are obtained in the usual manner by setting 


together with Eqs. (16). 


The perturbations of the energy levels of Eqs. (5) and (8) are obtained from the matrix elements 


by the usual second-order perturbation theory. 
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PHYSICAL REVIEW 


Secondary Emission from Films of Silver on Platinum* 
A. E. Hastincst 


VOLUME 57 


Brown University, Providence, Rhode Island 


Total secondary emission and energy distribution of 
secondary electrons have been measured for films of known 
thickness of silver on platinum. The depth of origin of the 
total secondary emission in silver as a function of primary 
energy, and the depth of origin of secondaries having a 
given energy were determined. Practically all the second- 
aries from primaries of 20 ev energy originate at a depth 


INTRODUCTION 


HILE the general characteristics of second- 

ary emission from pure metals have been 
well established, there is lack of agreement 
among investigations as to the depth at which 
the secondary electrons originate and as to the 
relative effects of the work function and the 
characteristic production of secondaries within 
the metal. Secondary emission is generally as- 
sumed to be dependent on the surface work func- 
tion through which the electrons must escape and 
on the production and absorption of electrons 
within the metal, these two effects being inde- 
pendent and characteristic of the particular 
metal. 

Treloar' concludes that a sufficiently thin film 
of a foreign material on a metal acts chiefly to 
alter the work function, the emission from the 
film itself being important only at greater film 


thickness. Plotting secondary emission against 


* Part of the dissertation presented for the Degree of 
Doctor of at Brown University. 

t Now at Naval Research Laboratory, Washington, D.C. 
*L, R. G, Treloar, Proc, Phys. Soc. 49, 392 (1937). 
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of less than 15 atomic layers, and those from primaries of 
50 ev energy originate at a depth of less than 30 atomic 
layers. For higher primary energies, it was found that an 
appreciable amount of emission comes from a depth greater 
than 150 atomic layers. Secondaries with energies close to 
that of the primaries originate at a depth small compared 
to that of the low energy secondaries. 


thickness of barium deposited on tungsten, he 
found an increase in the secondary emission as 
the work function was lowered to that of 
barium, the production of secondaries still oc- 
curring mostly in the tungsten. At greater thick- 
ness of barium film more of the secondaries 
originated in the film, which has a smaller pro- 
duction than that of tungsten, and the emission 
decreased, resulting in a maximum in the curve. 
The thickness at which half the emission came 
from the film and half from the base metal he 
estimated as 1.4 atomic layers, the mean depth 
of origin of secondaries from 300-volt primaries. 

Coomes? has not observed this maximum for 
thorium on tungsten, but Bruining* has observed 
it for barium on molybdenum. He found that 
the photoelectric emission also reached a max- 
imum for the same thicknéss of barium, indicat- 
ing that these changes were due in part to the 
decreasing work function. From experimental 
work with rough surfaces and from an experi- 
mental value of the absorption coefficient, 


2 E. A. Coomes, Phys. Rev. 55, 519 (1939). 
5H. Bruining and J. H. DeBoer, Physica 6, 941 (1939). 
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Fic. 1. Diagram of apparatus. 


Bruining‘ has calculated the mean depth of origin 
for secondaries from primaries of 500 ev energy 
to be about 14 atomic layers. He has also*-* 
investigated the mechanism of secondary emis- 
sion in detail. 

In the present work, measurements were 
carried out to determine the depth of origin of 
the total secondary emission in silver as a func- 
tion of primary energy, and the depth of origin 
of secondaries having a given energy. Total 
emission and energy distribution of secondary 
electrons were measured for films of silver on 
platinum. Energy distributions were obtained 
by the retarding potential method. Silver films 
were obtained from a calibrated source, so that 
the film thickness was determined absolutely 
from the time of evaporation. While the measure- 
ments were affected by changes in work function, 
resulting in the maximum in the curves de- 
scribed, no means for measuring the work func- 
tion directly were included in the apparatus. 


APPARATUS AND PROCEDURE 


The electron gun is based on the one described 
by Farnsworth,° and is shown in Fig. 1. F is an 
oxide-coated filament, surrounded by a guard 
ring. Electrons from F are accelerated through a 


4H. Bruining, Physica 3, 1046 (1936). 
5H. Bruining and J. H. DeBoer, Physica 4, 473 (1937). 
*H. Bruining, Physica 5, 17 (1938). 
7H. Bruining, Physica 5, 901, 913 (1938). 

(1998) Bruining and J. H. DeBoer, Physica 6, 823, 834 
*H. E. Farnsworth, Rev. Sci. Inst. 15, 290 (1927). 


circular disk and cylinder A, and again to g 
where the beam is defined by the 1.5-mm holey, 
C is a disk connected at filament potential t 
repel any electrons leaking around the gun, The 
beam is so defined that all the electrons Pass 
through the hole in the disk D and strike the 
platinum target T. The secondaries coming from 
T are collected by D, and the current to D j, 
measured by a galvanometer. The current to T 
and D is the primary current, and the ratio of 
secondary to primary current is called R. 

The energy analyzer is shown below the gun, 
Part of the secondary emission from T’ pasges 
through a hole in D and through a small shielding 
tube through the hole in E. The beam of second. 
aries is defined by-the target and the hole in E, 
and can pass through the disk H, without striking 
the sides of the hole, and into the shielded ool. 
lector G. A retarding field is applied between the 
cylinder E, and H and G. The current to G js 
measured with a d.c. amplifier, using an FP.54 
in a DuBridge and Brown circuit and a sensitivity 
of about 5X10-'* ampere per mm. Only thoge 
secondaries from T will reach the collector whose 
energy in electron volts is numerically greater 
than the retarding potential applied between the 
collector and target. Curve 1, Fig. 2, shows a 
typical curve of secondary emission obtained 
with this analyzer. The retarding potential must 
be corrected for the contact potential of the 
collector. This analyzer does not give a true 
distribution curve in energy, but actually gives 
the integral of the distribution curve between the 
energy at the point in question and the maximum 
energy which is the primary energy. The slope 
of this curve gives the true energy distribution 
curve, the number of electrons dm having energy 
between E and E+dE. Instead of plotting the 
complete distribution curves and taking slopes 
at certain energies, these slopes were determined 
directly by measuring the collector current at 
two retarding potentials, one either side of the 


point in question, and by dividing the difference - 


in currents by the difference in potentials. In 
this manner, a quantity called R., which is pro 
portional to the true number of secondaries 
having a given energy, was determined. The 
primary current was of course held constant 
throughout these readings. 

The analyzer was tested by placing it in front 


ST TS SS 


of the electron gun, so that their axes coincided. 
Then electrons from the gun could enter the 
analyzer and their energy distribution could be 

ined. Curve 2, Fig. 2, shows this distribu- 
tion. The lack of a very sharp drop at the right 
js due to the inhomogeneity in energy in the 
beam from the gun and to scattering in the 
analyzer. Curve 3 was taken with a small mag- 
netic field along the axis of the analyzer, and the 
resulting sharper drop indicates that there is a 
slight amount of scattering in the analyzer. 

The evaporator is built similar to one used by 
Farnsworth.!° A block of silver is completely 
enclosed in a molybdenum container except for 
a i-mm hole in the upper face through which 
silver can evaporate. The silver is heated by 
radiation from a flat spiral filament beneath it 
and enclosed in a molybdenum radiation shield. 
A Chromel-Alumel thermocouple imbedded in 
the side of the silver, in conjunction with a type 
K potentiometer, indicates the temperature, 
which is held so that the evaporation rate is 
‘constant within one percent. 

The silver source was calibrated in a separate 
tube by depositing on clean microscope cover 
glasses, which could be moved into the path of 
the vapor stream for known lengths of time with 
the source at 900 degrees centigrade. The source 


_was thoroughly outgassed in a high vacuum 


until the pressure was about 210-7 mm of 
mercury with the evaporator running, so that 
the rate of deposition would become constant 
before calibrating. In weighing on a micro- 
balance, sensitive to about 10~ g, the glass with 
the deposit was counterbalanced by a glass simi- 
larly treated. but not exposed to silver. The film 
was removed either by rubbing the glass with 
lens paper, or by washing in 50 percent HNOs, 
or both, all methods giving substantially the 
same result. The glass counterpoise was similarly 
treated. A maximum error of 5 percent was 
estimated for this determination. It was found 
that 4.65 X10-* gram of silver was deposited on 
a square cm in one hour, with the source 12.2 cm 
from the plates. For a (100) face, this corresponds 
to 21.6 atomic layers per hour, and for a (111) 
face, 30.6 atomic layers per hour. All data here 
are expressed in (100) units. 


“H. E. Farnsworth, Phys. Rev. 49, 605 (1936). 
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The electron gun is placed in the main tube, 
and the analyzer in a side tube sealed at 45 
degrees from the axis of the main tube. The 
evaporator is built in a small side tube and 
sealed to the experimental tube. The vacuum 
system consists of two Apiezon vapor pumps 
backed by a Cenco Hyvac pump. No grease or 
wax seals were used in the high vacuum side. 
The condensing traps between the experimental 
tube and pumps and between the fore pump and 
the vapor pumps were cooled with CO: snow in 
acetone. The experimental tube was kept sealed 
to the pumps throughout the experiments and 
the vapor pumps were run continuously. Helm- 
holtz coils were used to compensate the earth’s 
magnetic field. Pressure in the experimental tube 
was measured by an ionization gauge. 

Curves of R vs. film thickness were plotted for 
several primary energies, and curves of R, vs. 
film thickness were plotted for secondary energies 
with primary energy held at 300 ev. These curves 
show the effects of the film of varying thickness 
on the emission characteristics. 


RESULTS AND DISCUSSION 


Thin films 


The curves of Figs. 3 and 4 are typical of the 
results obtained. All the curves of R vs. film 
thickness, except that for 20-volt primaries, 
show a maximum. This is the behavior observed 
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Fic. 2. Curve 1, for secondary electrons from platinum, 
shows the variation of collector current with retarding 
—_ on the analyzer. Curve 2 shows the same variation 
or electrons coming directly from the . Curve 3 is 
similar to Curve 2, but was taken with a small axial 
magnetic field around the analyzer. 
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by Treloar for barium on tungsten, and is as- 
sumed to be due to the decrease in work function 
caused by a thin film of silver, since silver has a 
work function about two volts lower than that 
of platinum. The maxinum occurs at about one 
atomic layer, which is close to the optimum 
thermionic layer observed by Treloar for barium 
on tungsten. The maximum in R is smaller for 
curves taken at lower primary energies, which 
requires some explanation. During deposition of 
the first few atomic layers, two things are 
happening,—R is increasing because of the falling 
work function, and R is decreasing from the 
value characteristic of platinum to that charac- 
teristic of silver because of the emission from the 
silver film. The latter effect is more marked in 
the case of low primary energies, since there is 
smaller penetration of the low energy primaries. 
The former effect is more marked in the case of 
high primary energies, since there are many more 
secondaries with energy near that corresponding 
to the work function, and these secondaries are 
free to escape when the work function decreases. 
These effects tend to produce a greater initial 
increase in the curves of R at high than at low 
primary energies, and for 20-volt primaries, there 
is no maximum at all, but R decreases continually 
with film thickness. 

The maximum does not appear in the curves 
of R, vs. film thickness, since a change in work 
function of the order of two volts affects only 
electrons with very low energies in the metal, 
which are always excluded at the retarding po- 


tentials used. 
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Fic. 3. Secondary emission R (ratio of secondary to 
primary current) is shown as a function of atomic layers of 
silver on platinum for several primary energies. 
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R. increases with film thickness for secondaries 
with energy near that of the primaries, and de. 
creases for secondaries with low energies, because 
of the change in the secondary energy distriby. 
tion from that characteristic of platinum to that 
characteristic of silver. However, the increase at 
high secondary energies is much more sudden 
than the decrease at low secondary energies, 
which would mean that most of the secondaries 
of high energy originate at a smaller depth than 
those of low energy, and that most of the high- 
energy secondaries come from only a few atomic 
layers. This agrees with the results of Farns. 
worth”? for electron diffraction. He used second. 
ary electrons with energy very close to that of 
the primaries, and found that when the condi. 
tions for diffraction were established, the dif. 
fracted beam was entirely eliminated upon 
deposition of a few atomic layers of a foreign 
material. This indicates that secondaries with 
energy near that of the primaries do not originate 
at a very great depth as compared with the depth 
of origin of some of the low energy secondaries, 
This can be explained by saying that those pri- 
maries which penetrate the target to a depth of 
more than a few atomic layers are, because of the 
greater probability of scattering with loss of 
energy, nearly all effective in producing second- 


aries of low energy, while primaries scattered, 


very near the surface have a greater probability 
of being returned with small loss of energy. 


Thick films 

The curves of R in Fig. 3 for thick films give 
the maximum depth of origin of secondaries as a 
function of primary energy, since at the film 
thickness where the curve becomes flat all the 
secondaries are coming from the film. The curves 
indicate that practically all the secondaries from 
primaries of 20 ev energy originate at a depth of 
less than 15 atomic layers, and those from 50-volt 
primaries originate at a depth of less than 30 


atomic layers. For higher primary energies, the - 


curves do not flatten sufficiently to indicate the 
maximum depth from which secondaries come. 
Some data obtained out to 150 atomic layers 
show that the emission is still changing. 

The curves of R, (Fig. 4) at thick films all 
tend to flatten out at about the same rate with 
film thickness regardless of secondary energy. 
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For a primary energy of 300 ev the curves become 
very nearly flat at less than 20 atomic layers. The 
following explanation may account for the fact 
that secondaries with high energies appear to 
come from depths no greater than those at which 
secondaries with low energies originate. The 
number of high energy secondaries, coming partly 
from the base metal even at large film thick- 
nesses, is still changing at these thicknesses. 
These secondaries are effective in producing some 
low energy secondaries in the film, and thus the 
low energy secondaries continue to change as 
long as do the high energy ones. 

Since the total emission R is the integral of R, 
over the whole range of secondary energies, R 
should be affected by changes in R,, and R should 
change with film thickness as long as any curve 
of R, does. Actually R for primary energies above 
20 ev changes much longer than does any of the 
curves of R. for energies measured, for which 
there is no explanation. 

The slope of the R vs. film thickness curve 
would give the number of secondaries originating 
in a layer of thickness dx at depth x. However, 
because the slope changes so little beyond the 
first few atomic layers, and because of the experi- 
mental errors, it was not possible to determine 
this slope graphically with sufficient accuracy to 


be of any value. 


Effects of gas 
" While the pressure in the experimental tube 
was about 5X10-* mm of mercury during read- 
ings, it rose to about 1.5X10~7 mm of mercury 
when the evaporator was running. At this pres- 
sure the values of R for the outgassed target 
dropped as much as 10 percent in a few minutes, 
as determined by running the evaporator without 
exposing the target to the silver vapor. This 
rapid change did not occur with the gun filament, 
but not the evaporator, running, and was not due 
to contamination of the target by evaporation 
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several energies ¢ is shown as a function of atomic layers 
silver on platinum. The primary energy is 300 ev. 


from the filament. The initial maximum in the 
curves of R, Fig. 3, was not observed when silver 
was deposited directly on the outgassed target, 
without first exposing it to gases from the heated 
evaporator. The gas condensed on the target 
lowered R sufficiently to overcome the rise due 
to the silver film. Since it was not possible to 
outgas thoroughly the large block of silver in the 
evaporator, the procedure was to expose the 
clean target to gas coming from the evaporator 
until the emission no longer changed with ex- 
posure, and to take the value of R so obtained as 
characteristic of platinum for this work. The fact 
that with this procedure the maximum in R 
occurred seems to indicate that the surface layer 
of gas remains or is constantly formed during 
deposition, has a constant effect on the work 
function, and does not affect the interpretation 
of the results. Moreover, readings could be 
checked closely over long periods after the initial 
contamination of the target with gas had oc- 
curred. Several sets of curves similar to ‘those 
in Figs. 3 and 4 showed essentially the same 
behavior. 

The author is greatly indebted to Professor 
H. E. Farnsworth for his valuable criticism and 
guidance, and to Dr. John C. Turnbull for 
certain suggestions, during the progress of the 
experimental work. 
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PHYSICAL REVIEW 


Luminescence of Pure Radium and Barium Compounds 
D. H. KaBAKJIAN 


Department of Physics, University of Pennsylvania, Philadelphia, Pennsylvania 


Luminescence, excited in some carefully purified crystal- 
line barium and radium compounds by alpha-, beta-, and 
gamma-rays has been examined. It is found that these 
compounds show a faint luminescence at ordinary room 
temperatures. When the compounds are heated to various 
temperatures and then cooled, they show a continuous 
increase in luminescence with increase of temperature of 
heat treatment until a maximum is reached. Heat treating 
at higher temperatures diminishes the brightness. The 


(Received December 1, 1939) 


temperature of maximum brightness is not directly related 
to the melting point of the individual compounds, The 
maximum brightness in some cases was several hundred 
times the initial brightness of unheated crystals. No a 
changes in intensity of luminescence, such as might be due 
to changes in crystalline form, were detected. The : 
of these experiments on the modern theories of luminescence 
in crystalline solids is discussed. 


N a previous paper,! the writer described some 
experiments to prove that certain zinc borate 
compounds can be made luminescent in a state 
of purity. The compounds were obtained in a 
vitreous form by fusion in an electric furnace at 
about 1100°C. The clear zinc borate glass so 
obtained does not show any appreciable lumi- 
nescence. If, however, it is devitrified by heating 


for a few minutes at 600°C, a strong luminescence . 


of a purple color is developed. The structural 
change from the vitreous to crystalline state was 
verified by Debye-Scherrer photographs of the 
samples, which were kindly obtained for me by 
Professor B. E. Warren of the Massachusetts 
Institute of Technology. 

The close relationship between crystalline 
structure and luminescence in inorganic solids, 
indicated by the above experiment, is already 
well established. A general theory of luminescence 
based on the existence of electronic energy states 
in a crystal lattice, has been found quite suc- 
cessful in explaining a large number of experi- 
mental observations.? Most of the investigations 
are carried on, however, on activated crystals, 
i.e., on crystals containing an impurity atom. 

The possibility of luminescence in very pure 
crystals was doubted by most investigators until 
very recent times. Since the luminescence of 
some pure crystals is now generally admitted,’ 
the nature of the luminescent centers has become 


1D. H. Kabakjian, Phys. Rev. 51, 365 (1937). 

* F. Seitz, Trans. Faraday Soc. 35, 74 (1939). 

3N. Riehl, Ann. d. Physik (5) 29, 637 (1937); S. Roth 
schild, Zeits. f. Physik 108, 24 (1937); W. H. Byler, J. Am. 
Chem. Soc. 60, 1247 (1938). 
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a problem of theoretical interest. The writer 
has found that, with rays from a radioactive 
source, luminescence can be excited in a number 
of pure crystals. The initial brightness in most of 
these is rather faint, but in some cases it can be 
increased several hundred times by a proper heat 


treatment. 
Most impurity-activated phosphors also re. 


quire heat treatment in their preparation, but in . 


these cases the object is either to insure a dif. 
fusion of the impurity into the crystal lattice 
of the base material (solid solution) or to develop 
a crystalline structure in the mixtures when these 
are obtained in the form of fine precipitates. In 
the experiments described below, these considera- 
tions do not exist, since the compounds are’ 
obtained initially in crystalline form and no 
impurity is added. Raising the pure compound 
to a given temperature and then cooling produces 
a modification in its crystalline state which 
makes luminescence possible. It was considered 
desirable, therefore, to investigate the relation 
between the temperature of heat treatment and 
intensity of luminescence. 


MATERIALS AND MODE OF PROCEDURE 


The luminescence of the following compounds" 


has been investigated: RaBre, RaCle, RaSO, 
BaBr2, BaCle, BaSO,y. Pure RaBre was obtained 
from a mixture of radium bromide and barium 
bromide by a large number of successive crystal- 
lizations with distilled water and pure HBr. Its 
weight and radium-element content corresponded 
to pure anhydrous RaBre. Radium sulphate was 
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obtained from pure radium bromide by pre- 
cipitation with C.P. sulphuric acid. The pre- 
cipitate then was dissolved in hot concentrated 
sulphuric acid and was recrystallized by evapora- 
tion. The crystals were well formed and fairly 
large. Barium compounds were similarly purified 
by repeated recrystallizations. 

The pure compounds were placed in hard glass 
or silica tubes of 1.5-2.5 mm outside diameter 
and 0.2 mm wall thickness. They were heated 
to various temperatures in an electric furnace 
until equilibrium was established and _ their 
brightness was measured after cooling to room 
temperature. On account of their small size, the 
samples reach room temperature in a few min- 
utes. Tests showed that the luminescence was 
not affected by the rate of cooling. In radium 
compounds the excitation is produced chiefly by 
alpha-rays. In the case of barium compounds, a 
small quantity of radium bromide or radium 
sulphate was added. to the corresponding barium 
salts for an alpha-ray excitation. Since there is a 
possibility that the added radium might act as 
an impurity activator, these were investigated 
by a different process which eliminates this 
uncertainty. The barium salt was placed between 
the walls of two coaxial glass or silica tubes 
closed at the bottom. The thickness of the 
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Fic. 1. Growth of brightness, in arbitrary units, of pure 
radium compounds with the temperature of heat treat- 
ment. Curve I, anhydrous radium bromide. Curve II, 
radium sulphate. 
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TEMPERATURE OF HEAT TREATMENT ‘C 


Fic. 2. Growth of brightness of pure barium bromide 
with the temperature of heat treatment. (I) Barium 
bromide mixed with 0.24 percent radium bromide, (alpha- 
ray excitation). (II) barium bromide mixed with 1.2 percent 
radium bromide. (III) barium bromide excited with beta- 
and gamma-rays from a sealed radium tube. 


barium salt was about 0.6 mm. A thin-walled 
monel metal tube containing 6 mg of radium 
element in sulphate form was. fitted into the 
inner tube, having a wall thickness of 0.1 mm. 
The radium tube could be removed during the 
process of heating the salt. The excitation in this 
case is produced by beta- and gamma-rays of 
constant intensity since the radium is in equi- 
librium with its decay products. 

The results obtained by this method, which 
were in general similar to those obtained from 
mixtures of radium and barium compounds, 
showed that the production of luminescent 
centers is independent of the presence of radium 
atoms. 

The brightness measurements were made on a 
modified Lummer-Brodhun photometer, which 
has been described in a previous paper.‘ 


RESULTS 


A large number of determinations made on 


the temperature-brightness relationship of these 


*D. H. Kabakjian, Phys. Rev. 44, 618 (1933). 
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compounds gave similar results. Some typical 
curves are indicated below. Fig. 1 gives the vari- 
ation of brightness for pure radium bromide and 
radium sulphate with the temperature of heat 
treatment from 200° to 650°C. It is seen that 
the brightness rises as a continuous function of 
this temperature, reaching a maximum at about 
650°C for both compounds. This maximum ac- 
tivation temperature apparently is not related 
to the melting points of the two compounds. The 
ratio of maximum brightness to the initial 
brightness for radium bromide is 182 and for 
radium sulphate 35.4. This ratio refers to visible 
luminescence only. Since a large portion of the 
radium sulphate emission lies in the ultraviolet 
region, the emission energy of the two com- 
pounds might approach equality. No attempt 
was made to determine this point. 

Figure 2 gives the activation-temperature rela- 
tionship for barium bromide. Curves I and II 
refer to compounds containing 0.24 percent and 
1.2 percent of radium bromide, respectively. 
Excitation in these was produced chiefly by 
alpha-rays. 

Curve III was obtained by beta- and gamma- 
ray excitation from a sealed radium tube. The 
curves show the same activation-temperature re- 
lationship, although there are some differences in 
the two types:of excitation. The ratio of maxi- 
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Fic. 3. Growth of brightness of pure barium chloride 
with the temperature of heat treatment. Beta- and gamma- 


ray excitation. 


85885 


BRIGHTNESS 


wi 


| 


LOSS OF WEIGHT IN MG 


A . 
a i 
° 100 200 300 400 500 600 700 600 
TEMPERATURE °C 


Fic. 4. Curve I. Growth of brightness of pure BaBr;-2H,0 


with the temperature of heat treatment (beta- and ‘ 
ray excitation). Curve II. Loss of water of lization 
during heating. Ordinates to the right. crystal 


mum brightness in samples I and II is roughly 
proportional to their radium content and there. 
fore to the intensity of excitation. In Fig. 3 the 
temperature effect is shown for barium chloride 
also excited by beta-'and gamma-rays. 


CRYSTALLIZATION AT LOWER TEMPERATURES 


The crystals in the above experiment were 
produced by boiling the respective solutions to 
saturation. The bromides contain two molecules 
of water of crystallization. Complete dehydration 
of this compound required, under the conditions 
of these experiments, prolonged heating at 
200°C. The boiling point of the sulphate solution 
in concentrated sulphuric acid is about 190°C, 
therefore the initial brightness readings had to 
be started with a heat treatment temperature of 
200°C. In the case of bromides, it is possible to 
extend this part of the curve to a much lower 
temperature. A saturated solution of barium 
bromide was prepared at room temperature and 
placed in a desiccator. A few large clear crystals 


were produced in the solution by slow evapora- » 


tion extending over several days. It is assumed 
that these crystals approximated a perfect lattice 
better than those produced from hot solutions. 
The crystals were dried in a desiccator in the 
same way without heating. This compound, 
having a formula BaBr:-2H,O was ground, and 
its activation-temperature relationship was ex- 
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. 5. Variation of brightness with gradual removal of 
™ water of crystallization at 200°C. 


amined in the same way by beta- and gamma-ray 
excitation. The loss in weight due to dehydration 
was determined at the same time. 

In Fig. 4, Curve I gives the rise of brightness 
with the heat treatment temperature. Curve II 
shows loss of weight due to dehydration. The 
initial brightness in this case was very low and 
dificult to measure. The ratio of maximum 
brightness to the initial brightness was approxi- 
mately 300. 


EFFECT OF WATER OF CRYSTALLIZATION 


The removal of the second molecule of water 
of crystallization had a very marked effect on the 
intensity of luminescence. As the temperature 
was raised from 150° to 200°C, the luminescence 
began to decrease and continued to do so until 
all of the moisture was removed, when it sud- 
denly rose again. This relation is shown in Fig. 5. 
Asimilar effect might have been produced during 
the removal of the first molecule of water but was 
not detected because of the faintness of lumi- 
nescence at the lower temperature. The points on 
the smooth curve in Fig. 4 represent equilibrium 
values of brightness. 


THE NATURE OF EmIssION BANDS 


The compounds were sealed in clear silica 
tubes and their luminescence was examined by 
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means of a quartz spectrograph. The emission of 
the sulphates consists of a diffuse band extending 
from 5800-3000A with a maximum at about 
4050A. The bromides give a somewhat over- 
lapping double band with maxima at 4700 and 
3000A. The luminescence bands of the same 
compounds, after a small percentage of Mn or Sa 
was added to them were also obtained under 
similar conditions. The added activator in these 
cases produced its own band, quite distinct both 
in structure and wave-length, somewhat at the 
expense of the luminescence of the pure com- 
pound, but did not completely supersede it as 
in the case of zinc borate.' These are shown in 
Fig. 6. 

Although absolute purity cannot be claimed 
for any of these compounds, from the examina- 
tion of the emission bands of the “‘pure’’ and 
“activated” samples, it seems reasonable to 
assume that the wide diffuse bands are those 
due to the compounds themselves and not due 
to an undetected impurity. A more careful 


Fic. 6. Luminescence spectra of radium and barium 
compounds with comparison spectra of Hg arc. Upper 
(a) pure RaSO,; (6) RaSO,+Mn; (c) RaSOy+Sa. 
Middle panel: RaBr2; RaSO,. Lower panel: 
(a) pure BaSO,; (6) BaSO,+Sa. 
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Fic. 7. Decay of phosphorescence of pure barium bromide. 


examination of these bands and their dependence 
on the temperature of heat treatment seems 


desirable. 


RATIO OF PHOSPHORESCENCE TO FLUORESCENCE 


This ratio was determined for pure BaBrz by 
quickly removing the exciting source (the radium 
tube) and determining the decay rate of the 
residual brightness. This is shown in Fig. 7. 
Plotting the reciprocal of the square root of 
brightness against time gives a straight line 
indicating that the phosphorescence is a bi- 
molecular reaction. From the intercept of this 
curve, the ratio of phosphorescence to fluores- 
cence was found to be approximately equal to 5 
percent. 


DECAY OF LUMINESCENCE UNDER 
ContTINuOUS EXCITATION 


The decay of luminescence of pure radium 
sulphate was investigated. The decay was very 
rapid at the start, but soon became asymptotic 
with the time axis reaching a constant value and 
remaining so indefinitely. The decay curve is 
shown in Fig. 8. The curve is similar to those 
obtained for pure radium bromide by Rodman’ 
and pure barium bromide obtained by Smith.° 
It is very complex. Factors such as decomposition 
of the compound, change in the absorption coef- 
ficient of the material due to discoloration as well 
as the destruction of the luminescent centers 


5 J. Rodman, Phys. Rev. 23, 478 (1924). 
6. E. Smith, Phys. Rev. 28, 431 (1926). 
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probably enter into it. It has not been possible 
to find a simple explanation for this decay, The 
compounds will invariably recover their full 
luminescence on being heated to the tempera. 
tures of maximum activation. 


DIscussION 


When a crystal is irradiated by alpha-, beta. 
and gamma-rays from radium, as in the above 
experiments, the electrons can absorb sufficient 
energy to raise them to any unoccupied energy 
bands that may exist in the crystal, including the 
conduction bands (ionization). If these electrons 
could directly return to their normal levels, 
luminescence would result. The fact that only 
very faint luminescence can be obtained jp 
unheated, pure crystals may mean one of two 
things: (a) that permissible and unoccupied 
energy bands do not exist in pure unheated 
crystals, (b) that such bands may exist but the 
conditions necessary for the direct return of the 
electrons to their normal levels, with emission 
of luminous energy, do not exist. The heating of 
the crystals to higher temperatures, and then 
cooling, apparently produces the structural 
changes necessary for luminescence. It is dif- 
ficult to determine from these experiments the 
nature of this change. It cannot be due to a 
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Fic. 8. Decay of luminescence of pure radium sulphate, 
Curve I time in minutes, Curve II time in hours. 
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chemical decomposition, as the compounds 
remain perfectly white after heating. However, 
there is distinct evidence of instability, indicated 
by the fact that heated radium bromide becomes 
almost black in time, while the unheated speci- 
mens show only a slight change in color. 

Seitz’ explains the luminescence of heat- 
treated pure ZnS as due to Zn atoms diffusing 
into the interstices of the crystal lattice and 
acting in the same way as impurity atoms. Asimi- 
lar displacement may be taking place in the case 
of barium compounds. This would explain the 
brightness-temperature curves if it be assumed 
that the number of displaced atoms depends on 


TF. Seitz, J. App. Phys. 10, 479 (1939). 


the temperature of heat treatment. It is still 
necessary, however, to explain why these give 
a much broader emission band than other im- 
purity activators as shown in Fig. 6. A more 
thorough examination of the emission bands and 
their dependence on the temperature of heat 
treatment seems desirable. 

The author wishes to acknowledge his indebt- 
edness to Dr. G. P. Harnwell and Dr. F. Seitz 
for many valuable suggestions and criticisms. 
This research was aided by a grant from the 
Faculty Research Committee of the University 
of Pennsylvania and also by a grant from the 
American Philosophical Society which are hereby 
gratefully acknowledged. 
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A band spectrum attributed to the molecule SbCl, in the region \\4200-5600, has been 
excited by introducing SbCI; into active nitrogen. No other bands appeared in the range 42000 
to 7000. The spectrum consists of bands showing pronounced progressions and degraded 
towards the red. They can be placed in two sub-systems of bands, designated as A; and As». 
Vibrational analysis shows that within experimental error the wave numbers of the heads of 
forty-one bands of SbCI® are represented by the formulas: 


A, : =22395+244.4u’ —2.28u"— 368.0u" +-0.84u', 
Ag : ve=25855+240.2u’ —2.19u’?—370.0u" + 1.00u'". 


The chlorine isotope effect appears in 13 bands and is used to confirm the v’, v’’ designations, 
The analysis of the spectrum given here is somewhat different from that originally presented 
two years ago. Many more bands are here accounted for and identified. An attempt failed to 
produce a spectrum of AsCl by introducing AsCl; vapor into active nitrogen. 


HE electron configurations of the lighter 

diatomic molecules have already been de- 
termined with a fair degree of certainty. Com- 
paratively few measurements have thus far been 
made of the spectra of the heavy diatomic mole- 
cules. Spectra have been produced, however, for 
many of the diatomic halides of the elements in 
group IVb of the periodic table. For group Vb 
there are only the spectra due to SbF ! and the 


four bismuth halides.2 This experiment was 


'G. D. Rochester, Phys. Rev. 51, 486 (1937). 
*F. Morgan, Phys. Rev. 49, 41 (1936). 


undertaken to find out whether the same regu- 
larities exist for the halides of group Vb as have 
been shown to exist for group IVb.* 

The arrangement of the apparatus and the 
methods employed for the excitation of the bands 
of antimony chloride are essentially the same as 
those originated by Strutt and Fowler‘ in their 
work with active nitrogen, and as further de- 


* W. Jevons, L. A. Bashford and H. V. A. Briscoe, Proc. 
Phys. . 49, 532 and 537 (1937); W. Jevons and L. A. 
Bashford, Proc. Phys. Soc. 49, 554 (1937). 

a i J. Strutt and A. Fowler, Proc. Roy. Soc. A86, 105 
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Fic. 1. The spectrum of SbCl; (a) low dispersion, (b) 
higher oy peers 4200 — 5600, with vibrational quantum 
numbers for systems A; and A». 


veloped by Mulliken’ in his early work on the 
band spectra of the copper halides. By introduc- 
ing the vapor of pure antimony trichloride into 
a stream of active nitrogen, the spectrum at- 
tributed to the molecule SbCl was excited. The 
wave-length region 2000 to 47000 was photo- 
graphed and examined. 

A Hilger quartz spectrograph of low dispersion 
(200 A/mm at \5000) was used in the preliminary 
stages with Eastman 40 plates. With exposures 
varying up to 70 minutes only one group of 
bands, between the wave-lengths \4200-A5600, 
was disclosed. Using a glass spectrograph with 
lenses of 492-mm focal length, dispersion 42 
A/mm at 45000, this system of bands was photo- 
graphed on Ilford S. G. panchromatic backed 
plates in ten hours. Intensities were visually 
assigned. The iron arc was used for wave-length 


standards. 
The bands shade towards the red. Enlarge- 


5R.S. Mulliken, Phys. Rev. 26, 1 (1925). 
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TABLE I. Wave numbers of SbCl® band heads, in cm~; A, above A>. 


HUDES 


ments of the photographs, Fig. 1, illustrate their 
main features. Fig. la was taken with the 
quartz spectrograph and shows clearly the exist. 
ence of the group of bands in the region dA4209- 
5600. Taken with the glass instrument of greater 
dispersion, the general characteristics of the 
bands are more easily recognized from Fig. 1b 
This figure contains the assignments of some of 
the vibrational quantum numbers in the two 
systems, A; and Az (see below). 

The heads of forty-one bands which have been 
observed can be arranged in two sub-systems, 
designated A; and Az, where A; identifies the 
less refrangible system (smaller value of »,). The 
wave numbers of the band heads are shown jp 
a —v” array in Table I, in which are indicate 
the intensities in parentheses. A; bands ar 
placed above A: bands with 6», between. Several 
band heads which were impossible of identi. 
fication, with their estimated intensities, fg. 
low: 23970(2); 20141(3); 19789(3); 19490(3). 
18838(3) ; 18114(4) ; 18088(3) ; 17726(5) cm-, 

Within experimental error the heads of thes 
two sub-systems of bands of SbCI* are repre. 
sented by the formulas: 


vy = 22395 + 244.40’ — 2.28u" — 368.0u"’ +0.84y', 
ve= 25855 + 240.2u’ — 2.19’? — 370.0u"’ +1.00u", 


where u=v+}. 
The average of the differences between ob- 


served values of the band heads of SbCI* and 
values calculated by the formulas is +3.1 em+ 
The constants are different from those originally 


Means 


4 
v =3 4 5 6 7 8 9 10 11 12 13 14 act 
0 21248(1) 20880(2) 20516(3) 20157(4) 19802(4) 19446(4) 19093(4) 18743(5) 
3456. 3455. «34543455 


1 21485(2) 21119(2) 20758(3) 20394(4) 20038(1) 19694(0) 
3457 3458 3447 
23851(2) 23496(0) 23141(0) 


21717(1) 21358(3) 


21943(0) 21591(2) 


1 361.0 362.5 361.0 = 


23617(3) 23258(4) 22901(4) 22547(4) 22198(4) 


23375(0) 23021(0) 22664(0) 22319(1) 


21852(4) 21503(2) 21161(2) 20823(0) al 
2363 

18980(2) 

3454 

22434(0) 22087(2) 21743(2) 21397(1) 21054(2) ia 


212870 
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‘ven in a note on the spectrum of SbCI.* Whereas 
only 28 band heads were identified in the earlier 
note, 54 band heads (including isotope heads) 
are taken into account here and identified. 

It was to be expected that the band heads due 
to the individual antimony isotopes would not 
be resolved with the instrument here employed. 
Weaker heads of SbCI*’ are displaced toward the 
system origin from stronger SbCI® heads; the 
displacement coefficient, p—1, is —0.0210. The 
origins of the A, and Az» systems are found by 
extrapolation to zero displacement. The dis- 
placements of 13 bands due to the SbCI*’ isotope 
shift are given in Table II. 

The agreement of observed and calculated 
values of the isotope displacements confirms the 
locations of the system origins (A; and Az). Also, 
the relative intensities of the isotope bands ap- 

satisfactory for chlorine isotopes. 

Evidence for the identity of the emitter of this 
spectrum was weighed. Vapor of SbBrs intro- 
duced into active nitrogen produced only the 


 NBr spectrum.’ Vapor of Sb metal produced no 


spectrum in active nitrogen. The presence of 
SbN seemed therefore to be eliminated. If, 
nevertheless, SbN (or SbO) were the emitter, no 


_ isotope displacement would be observed. So an 


alternative classification of the band heads was 
speculated upon, by considering the above band 
heads of the isotopic molecule SbCI*’ as R heads 
of SbCl* displaced from more intense Q heads of 
SbCl*. This appeared reasonable but was pre- 
cluded, however, because w,’—w,”’ is large. Each 
band would then be characterized by a small 
interval between Q and R head.* This was not 
the case here. 


*I. Hudes, Phys. Rev. 52, 1256 (1937). 
7A. Elliott, Proc. Roy. Soc. A169, 469 (1939). 
Pe. W. . Report on Band Spectra of Diatomic Mole- 
, p. 63. 
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From the similarity of vibrational constants of 
these two band systems it is probable that they 
comprise a doublet system. It is not possible 
from the present analysis to determine in which 
state the multiplicity occurs. The existence of a 
third system, making a triplet, is not impossible ; 
because an odd multiplicity of electronic states 
is expected for a molecule possessing an even 
number of electrons. The above statements have 
likewise been made in reference to SbF. 


TABLE II. Displacements of bands due to SbCl" isotope shift. 
Avcarc. ts 


Ai A: (cm~) 
0,5 41. 40 0,6 45. 47 
0,6 49. 47 0,7 $2. 54 
0,7 50. 54 0,8 64. 62 
0,8 57. 62 0,9 68. 69 
0,9 60. 69 0,10 76. 76 
1,5 36. 34 1,8 56. 57 
1,6 36. 42 


For group Vb halides a comparison between 
the constants of the antimony halides and the 
bismuth halides was made. The constants of the. 
bismuth halides were plotted against the number 
of electrons in the molecule and were then com- 
pared to similar plots for the antimony halides. 
The expected trends were found in both cases, 
namely, with increasing number of electrons and 
mass of a molecule there is a decrease in the 
values of w, and X.w, in both the upper and lower 
states, a decrease in v., and an increase in the 
value of 5y,. Similar trends were found also in the 
values of the constants of any one halide of the 
two metals. 

The introduction of AsCl; into streaming ac- 
tive nitrogen produced no characteristic band 
spectrum in the region \A2000—7000. 

Appreciation is expressed to Thomas J. W. 
O'Neil for his assistance with this experiment. 
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Spectra of Lead Hydride and Tin 


high pressure arc are briefly discussed. 


Revised calculations of energy constants of the PbH molecule are presented as a result of 
analyses of bands in the photographic infra-red. Former v’ values are raised by one unit. An 
isolated weak PbH band is found at 3815A. Additional SnH bands in the red produced in the 
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Hydride in the Near Infra-Red 


Wititi1am W. Watson AND RALPH SIMON 
Sloane Physics Laboratory, Yale University, New Haven, Connecticut 


(Received February 2, 1940) 


N previous reports! *? on the spectra of PbH 
and SnH in the visible, the probable existence 
of spectra of both these molecules in the photo- 
graphic infra-red was mentioned. With further 
development of the high pressure arc in hydro- 
gen, together with long exposure times and wider 
slit widths, PbH bands have been recorded to 
9105A and SnH bands to 7234A. In addition, 
careful search for the PbH analog of the violet 
2A—*II SnH bands reveals a single weak band 
centered at 3815A. This band is not well de- 
veloped and is not accompanied by another at 
any plausible interval for a *II PbH state. 
The Pb and Sn arcs have been run in hydrogen 
at 4 to 5 atmospheres pressure. All spectrograms 


1W. W. Watson, Phys. Rev. 54, 1068 (1938). 
2W. W. Watson and R. Simon, Phys. Rev. 55, 358 


(1939). 


TABLE I. Additional bands in the 2=—>*= system of PbH. cm= units; d denotes fused lines. 


were obtained in the second order of a 21-foot 
grating, average dispersion 2.13A per mm, except 
those of PbH for wave-lengths longer than 72004. 
For the latter the first order with 4.94A per mm 
was used. Because of the wide-open structure of 
these hydride spectra, however, this lower dis. 
persion was quite sufficient despite the cop. 
siderable line widths due to pressure broadening 
and to the use of a wide slit (0.1 mm). 

It was thought that bands involving a new 
electronic state of PbH might be located in the 
infra-red, since strong interaction with other 
states is revealed' by the details of the red 
2y—"> system. We find, however, that the ex- 
tensive PbH spectrum in the infra-red to 91004 
belongs completely to this same system. With 
the inclusion of the data from the analysis of 
these new bands, plus additional lines in previ- 


(0,3) (0,4) (1,3) 
P; Ri Re 
11,831.6 
13,072.8 13,128.7 13,090.6 751.9 768. 819.4 |13,555.7 13,575.7 
042.9 109.6 063.2 13,124.6 723.4 7904 742.7 803.3 527.4 5514 
462. . 0 008.9 086.1 031.9 102.3 691.7 768.2 712.9 783.4 522.3 
344.94 431.5d 373.3 453.9 |12,970.9 056.9 12,997.2 076.5 656.4 742.7 679.7 759.8 460.2 489.4 
9 301.0d 396.7 331.9 420.5 929.2 024.5 957.4 046.3 617.5 712.9 643.1 732.5d| 420.4d 13,520.4 452.7 13,4478 
10 251.7 357.0 285.1 383.6 883.4 12,988.0 913.8 012.3 574.8 679.7 602.7 701.1 376.5 486.9 412.2 = $17.1 
ll 199.5 312.6d 234.8 342.2 833.5 947.3 866.4 12,973.9 528.4 642.3 5584 666.3 329.5 449.1 367.8 
12 264.5 180.3 296.2 779.6 902.4 814.9 931.3 478.5 601.2 510.5 626.7 278.3 407.4 319.6 4427 
13 212.0 121.1 245.5 721.9 853.5 759.6 884.4 424.6 556.2 458.7 583.7 223.5 361.6 267.5 4 
14 013.5 154.5 057 191.5 669.1 800.3 699.7 833.5 367.2 507.6 403.0 537.0 164.9 311.7 211.2 3517 
5 |13,943.8  %090.8d 13,989.4 131.6 594.3 743.2 636.1 778.3 306.0 454.9 343.6d 486.0 102.4d 257.4 150.8 
16 026.5 917.2 066.8 524.4 681.7 568.3 719.0 241.2 398.6 280.6 431.3 035.7 199.2 086.1 2436 
17 790.4 13,956.2d 840.3 13,999.5 450.6 616.1 496.4 655.1 172.94 338.2 213.9 372.6 |12,964.8 136.8 018.0 183.3 
18 707.0 880.5 759.2 926.2d) 3725 546.3 420.1d 587.2 100.6 274.1 143.3 309.9 890.0 068.4 12,945.2 118.0 
19 619.6 800.9 673.6 848.4 290.7 4724 340.3 515.0 024.5 206.1 069.0 243.5 811.4 12,997.2 868.3 
20 527 716.8 583.7 766.3 | 204.9 394.2 256.2 438.7d 134.3 10,990.3 172.9d| 726.8d 921.8 787.2 12,973. 
21 431.2d 628.2 489.4 679.8 114.8 312.0 358.1 058.6 098.9 640.3 841.7 896.2 
22 535.1 391.0 588.5 021.1 6 075.8 273.4 10,979.4 020.9 548.8 757.3 611.2 8137 
23 438.2 6 11,923.4 135.4 11,979.7 184.6 453.2 668.2 517.9 
24 336.3 181.1 821.8  O41.1d 879.6 091.7 4 574.7 420.1d  635.0d 
25 230.3 716.4 11,942.4 775.7 11,994.8 249.4 476.9 317.6 530.0 
26 8$40.0d 667.8 894.4 141.1 374.7 211.0 438.74 
27 733.4d 4 268.0 100.6d 333.7 
28 11,911.9 11,985.4 224.6 
29 791.2 O41.1d 866.3 11084 
30 666.3d 11,921.9 11,9929 
31 7.1d 
32 404.5 7452 
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OF LEAD HYDRIDE AND TIN HYDRIDE 


TABLE II, Corrected PbH band origins. 


2 


3 4 


1 

0 

| 17,9774 1503.1  16,474.3 
448.2 448.1 

2 | 18,425.6 1503.2 16,9224 1444.3 
438.9 438.1 

3 | 18,864.5 1504.0 —17,360.5 
432.6 

4 | 19,297.1 

5 | (19,700) 


14,551.1 1388.5 13,162.6 1324.9 11,837.7 
479.0 478.6 

15,030.1 1388.9 13,641.2 
448.0 


15,478.1 


ously analyzed bands, revision of the calculation 
of the energy constants' for PbH is necessary. 


Rep AND INFRA-RED PbH SpectrRuM 


Assignments of the lines of four additional 
bands at the red end of this system are included 
in Table I. This analysis is based on the usual 
combination relations and indicates that the v’ 
numbering of all the bands described in reference 
1 should be increased one unit. The results of a 
recomputation of all band origins are given in 
Table II. Since we can find no traces of bands 
that would extend this array upwards, we believe 
that these vibrational quantum numbers are now 
correct. New constants of the energy terms are 
collected in Table III. 

These data emphasize the peculiar nature of 
the upper *2 state as discussed in reference 1.* 
The new By’ being markedly lower than the old 
one, the tendency for the B,’ values to increase 
with » to a maximum and then to decrease is 
even more apparent. Since about 500 cm are 
now added to the extent of the stable levels of 
this state, the cut-off of levels comes actually at 
about 0.42 volt above the bottom of the po- 
tential energy curve. These and the other pe- 
culiarities of this *2 state all point to the existence 
of additional near-lying electronic states of PbH. 
The faint band at 3815A is, however, the only 
indication we can find of other PbH spectra. 
Because of competition from the continuous 
background always present in high pressure 
sources and which is particularly strong to the 
violet of the 4019 PbI atomic line, long exposure 
times do not yield better registrations of this 
band. For this reason and also because of its 


isolated character we have not attempted its 


* Also cf. L. Geré, Physica 7, 155 (1940). 


TABLE III. Constants of the energy terms for PbH 
in cm units. 


CONSTANT 
2.478 4.899 
B,* 2.660 4.754 
B,* 2.766 4.609 
B;* 2.770 4.467 
Bf 2.646 4.324 
B® 4.971 
a 0.144 
To 1.848A 
Do — 2.66 10-4 —2.01x10-* 
We ~520 1564.1 
XcWe 29.75 
Ve ~18018 


analysis. We conclude that the only stable states 
of the PbH molecule are those involved in this 
red band system. 


Rep SnH SpectRUM 


SnH bands have been photographed further 
into the red than those already reported.? Promi- 
nent heads of the red-degrading branches occur 
at 6745, 6892 and 7030A with a pile-up of lines 
at 6931A. Between the latter and a long series of 
lines on its violet side, as well as between two 
other branches, we find *II-state combination 
differences.? We have not found differences relat- 
ing these bands to the bands at 6095 and 6214A, 
however, and since the analysis of all these bands 
is far from complete we refrain from presenting it. 
They are definitely of a different type than those 
of the PbH system, probably *2—*II, judging 
from the number of branches. It is altogether 
probable that more SnH bands exist further into 
the infra-red. Their measurement should prove 
helpful in the analysis of the bands already 
photographed, but we find it impossible to obtain 
high dispersion spectrograms in the infra-red 
with the present high pressure Sn arc as a source. 
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The index of refraction of methane has been measured for 49 different wave-lengths in the 
infra-red between 1 and 15y4. These data have been interpreted in terms of the active infra-red 


frequencies of methane, and yield what may be considered estimates of the dipole moment of 
the CH bond and the rate of change of this moment with CH distance. The measurements 
show that the infra-red bands are only about one-sixth as strong as would be necessary to 
remove the discrepancy between ”.—1 as obtained from indices of refraction measured at 
optical wave-lengths and extrapolated to infinite wave-length, and e—1 as determined from 


dielectric constant measurements. 


HE need for quantitative information on 
the intensities of infra-red bands of molec- 
ular gases has frequently been pointed out. Van 
Vleck! has called attention to the fact that such 
data would permit quantitative calculation of 
the contribution of ‘‘atomic polarization’’ to the 
dielectric constant, thus furnishing a check on 
dielectric constant measurements and increasing 
the accuracy of measurements of dipole moments. 
Dennison? and Mulliken* have used intensity 
measurements to calculate ‘‘effective charges” 
of the vibratng molecules (or the rate of change 
of the dipole moment with internuclear distance) 
and have pointed out the need for more quan- 
titative data. This paper contains essentially the 
same information as would be obtained from 
accurate measurements of the intensities of the 
infra-red bands of methane. Since the reasons for 
measuring the refractive index instead of the 
intensity of absorption have already been 
pointed out by Van Vleck,' and have been fully 
explained in a previous communication,‘ they 
will not be repeated here. 

The method consists in measuring the change 
in deviation produced by a change in pressure in 
a hollow prism containing the gas being studied. 
The apparatus used was essentially similar to 
that used for the measurements on HCl,‘ except 
for the thermopile. The thermopile, mounted on 


1J. H. Van Vleck, The T of Electric and Magnetic 
Susceptibilities (Oxford, 1932), Chapter 3. 
?D. M. Dennison, Phil. Mag. 1, 216 (1926); Astrophys. 


2R. S. Mulliken, J. Chem. Phys. 2, 400; 2, 712 (1934). 
a 935) Rollefson and A. H. Rollefson, Phys. Rev. 48, 779 
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a comparator carriage, is used to locate the 
image. A new design of thermopile increased the 
speed of taking observations by a factor of about 
ten and decreased the uncertainty of measure. 
ments by a factor of six. 

The new thermopile has two receivers, each 
10 mm by 1 mm, mounted side by side, with the 
edges slightly overlapping but not in contact. 
Each receiver is mounted on a two-junction pile 
and the two piles are connected in opposition, 
Fig. 1 shows a typical curve of galvanometer 
deflection against position as this thermopile is 
moved across a slit image. In practice it was 
necessary to take only four or five points along 
the straight center portions of the curve, and the 
two zeros, in order to determine the image posi- 
tion. In the high energy regions of the spectrum, 
successive determinations of the image position 
would check to within 0.3 micron. For all the 
wave-lengths less than 11, variations of more 
than 2 microns were rare. The thermopile case 
was so constructed that thermal drift was 
eliminated, and since the thermopile was of the 
compensated type and the two compensating 
piles were almost identical, deflections due to 
electrical surges induced in the galvanometer- 
thermopile circuit by sparks in neighboring 
rooms were reduced to a minimum. A Kipp and 
Zonen type Z C galvanometer was used for al 
the measurements. When necessary, the deflec- 
tions were magnified by a thermal relay. With 
this arrangement, when working at night, it was 
possible to reach the limit set by Brownian 
motion of the primary galvanometer. 
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INDEX OF REFRACTION OF METHANE 


The methane for most of the measurements 
was obtained by heating sodium acetate and soda 
lime. The gas thus obtained was purified by 
repeated distillations in an all-glass system. The 

was generated several times and the agree- 
ment of the refraction measurements on the 
different samples indicated that the gas was 
sufficiently pure. However, since the absolute 
values of the indices of refraction which were 
obtained did not agree with the value obtained 
by extrapolating data available for wave-lengths 
in the visible region of the spectrum, several 
additional checks were made. Since the same 
tus had been used for measurements on 

HCI and CO, gases which are easily obtained in 
high purity and which had given results in close 
agreement (better than } percent) with those 
obtained by extrapolating the data of Cuthbert- 
son and Cuthbertson, it was decided that the 
most likely source of error in the measurements 
on methane was impurity of the gas. The vapor 
pressure of the methane was then measured at 
the boiling point of oxygen. The value obtained 
fitted on the vapor pressure curve as well as 
those tabulated in the International Critical 
Tables. This indicated that no volatile impurities, 
such as hydrogen, carbon monoxide, or air, could 
be present in sufficient quantity to cause a notice- 
able change of index of refraction. Then varying 
fractions of the methane condensed in the liquid- 
air trap were distilled into the hollow prism and 
the index of refraction measured for each frac- 
tion. No change in the index reduced to standard 
conditions was noted. This indicated that if any 
impurity were present it must have a vapor 
pressure curve almost identical with that of 
methane, or an index of refraction nearly equal 
-to that of methane. An additional check on the 
less volatile impurities was furnished by the 
infra-red absorption spectrum of the gas used. 
Absorption measurements were customarily 
made throughout the spectrum in order to deter- 
mine for which wave-lengths there was appreci- 
able absorption, and therefore where measure- 
ments of the index of refraction could not be 
made by the prism method. No bands other than 
those belonging to methane were observed. It 
must be pointed out, however, that this consti- 
tutes merely corroborative evidence and not 
conclusive proof of purity, since the resolving 


711 


power which was used was not high and some of 
the organic impurities which may have been 
present have their strongest bands in the same 
regions as those of methane. 

When preliminary efforts at interpretation 
indicated that the data were not consistent with 
measurements in the visible, new runs were 
made for which the methane was generated by a 
slightly modified Grignard process. Both this 
process and the chemicals used in it were entirely 
different from the sodium acetate-soda lime 
process. The results of the new runs agreed with 
the old to within 0.5 X 10~ for the absolute value 
of m, and so far as shape was concerned the new 
curves were indistinguishable from the old. 
Shortly thereafter new measurements® on the 
index of refraction of methane in the visible and 
ultraviolet appeared. These new values for n 
gave a dispersion curve which when extrapolated 
into the infra-red gave values much closer to 
those we obtained, although the disagreement 
was still 4X10~-*, or about 1 percent in the ab- 
solute value of m—1. The slope of the new extra- 
polated curve allowed a much more satisfactory 


DISPLACEMENT 
0.325 


IN MM 
0.350 


0.375 


Fic. 1. Galvanometer deflection as a function of thermo- 
sey position. The upper diagram is an approximate curve 
or a large total displacement, while the lower curve is one 
actually obtained in the measurements on methane. With 

ints as good as these, the position for zero deflection can 
be determined to Iu or better if the comparator screw is 
reliable. 


*S. Friberg, Diss. Lund. (1933) quoted by T. Larsen, 
Zeits. f. Physik 105, 166 (1937). 
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WAVE-LENGTHS IN MICRONS 


Fic. 2. Dispersion of methane. The solid curve is that 
given by theory; the points are averaged experimental 
results; the dotted curve is (n—1) extrapolated. 


interpretation of our results than the old. Even 
if the 1 percent difference which remains were 
due entirely to systematic error in the present 
experiments or impurity of the gas used, it 
would be of little importance in the interpreta- 
tion of the results, since the changes of index 
which are important probably cannot be meas- 
ured to within one percent. Any impurity which 
had active frequencies in the infra-red would of 
course cause trouble but would probably have 
been detected by its absorption. 

The effect of purity of spectrum was tested by 
trying several different slit widths in the mono- 
chromator and gas spectrometer. The trials 
showed that changing slit widths over the range 
actually used had no effect. The values of the 
index of refraction which were obtained were 
reduced to standard conditions (0°, 760 mm) by 
means of coefficient of expansion and pv data 
given in the International Critical Tables. 
(a=0.003683; the pv correction was at most 


about } percent.) 


RESULTS 


The index of refraction of methane was 
measured at 49 different wave-lengths between 
1.6 and 14.8u. The values of the index for the 
various wave-lengths are given in Table I and 
plotted in Fig. 2. The values given are averages 
for several runs. In taking these averages, runs 
which did not cover the whole range of wave- 
lengths used were adjusted to agree in absolute 
value with the average of the runs which included 
all wave-lengths; that is, the whole curve was 
shifted up or down a definite amount. This was 
obviously necessary in order to avoid spurious 
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kinks in the dispersion curve which would be due 
only to the fact that certain regions of the 
spectrum were more thoroughly investigated 
than others. The largest shift was 1.1 X 10-4 jp ai 
and most of the curves were not shifted at al. 

The contribution of the overtone bands in the 
neighborhood of 2.3u is small but easily notice. 
able. The measurements in this region Cannot, 
however, yield quantitative values of the inten. 
sities of the overtones. Two bands overlap in this 
region, and in order to get estimates of their 
separate contributions to the index of refraction, 
measurements would have to be taken in regions 
very near the band center where there is appre. 
ciable absorption. Such measurements would not 
be reliable if obtained by the prism method, 
since unsymmetrical absorption of the slit image 
causes an apparent shift of the image. Since 
measurements in this wave-length region are not 
quantitative, they are not included in Table | o 
Fig. 2. 

The only other measurements on the dispersion 
of methane in the infra-red which have come tp 
the writers’ attention are those of Koch® who 
used an interferometer and residual rays of 
wave-lengths 6.5574 and 8.678u; that is, wave. 
lengths on each side of the 7.74 band. His values 
for (n—1)X10® were 419.20 and 450.06, respec. 
tively, for the above two wave-lengths, while 
those obtained in the present work were 426.35 
and 458.10. The quantity which is of most im. 
portance for the interpretation of the results is 
the difference between the values obtained for 
the two wave-lengths. For this difference Koch 
obtains 30.86 while we obtain 31.75. This must 
be considered satisfactory agreement, since we 
differ by 1.7 percent in the same direction on the 


absolute values. 


INTERPRETATION OF RESULTS 


The contribution of an infra-red band to the 
square of the index of refraction is given by 
A;(n?—1) = K;/(v2—v*) where »; is the resonane 
frequency, v the frequency of the light used, and 
the K; are constants. A(m +1) due to the infra-red 
bands was obtained by subtracting the exper 
mental values of m from the values obtained by 
extrapolating into the infra-red the measure 


® J. Koch, Nova Acta Soc. Upsala, 2, No. 5 (1909). 


| 


be due 
Of the 
tigated 
inn; 
t all. 
in the 
notice. 
cannot, 
> inten- 
) in this 
of their 
raction, 
regions 
appre- 
uld not 
method, 
it image 
e. Since 
| are not 
ble I or 


spersion 
come to 
ch*® who 
rays of 
S, wave- 
is values 


426.35 
most im- 
results is 
ained for 
ice Koch 
‘his must 
since we 
on on the 


INDEX OF REFRACTION OF METHANE 


ments made by Friberg® in the visible and ultra- 
yiolet. The values of the K; for the two active 
bands were adjusted by successive approxima- 
tions until each K; remained constant throughout 
the spectrum. In this adjustment the extra- 
ted curve was adjusted vertically to give the 
same value of K on both sides of a band, the 
assumption being that the visible and ultraviolet 
data gave the slope of the curve more accurately 
than the absolute value. 
The final results are given in Table I. The 
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average value of K; for the band for which 
v3= 3019.6 cm was 16.50. The average value of 
K, for the band »4=1304 cm was 7.50. For the 
computation of these averages the various values 
are weighted according to the magnitude of 
A(n—1), and values beyond 11 are not used 
since the energy available in this region is very 
small and the points consequently scatter rather 
badly compared with other regions of the 
spectrum. The values for some wave-lengths very 
near the absorption bands (those in parenthesis) 


Tape I. Index of refraction of methane. Here (n—1)ezp is the average experimental value of n—1 where n is the index 
of rdrecion (n—1)ext ts the value obtained by extrapolation of the data of Friberg; 43 and A, are the contributions of v; and 
», to the index of refraction, as determined experimentally. 


WAVE-LENGTHS 


A(n —1) X 10° 


As *(v3? (" —Icarc X 10° 


437.07 
436.74 
* 436.40 
436.35 
436.30 
436.30 
436.17 
436.15 
436.12 
436.11 
436.08 
436.07 
436.06 
436.04 
436.02 
436.02 
436.01 
436.00 
435.99 
435.98 
435.98 
435.97 
435.97 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 
435.9 


436.50 
435.70. 
432.70 
431.50 
427.80 
425.85 
445.25 
442.80 
441.10 
440.10 
437.80 
437.30 
436.80 
435.80 
434.55 
434.00 
432.60 
432.05 
429.05 
426.30 
421.8 
416.8 
413.6 
466.5 
463.6 
462.0 
460.3 
459.0 
457.7 
456.0 
454.7 
453.8 
453.4 
452.2 
450.9 
450.6 
449.6 
448.6 
447.6 
447.7 
446.5 
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431.25 
428.34 
427.22 
445.38 
442.99 
441.07 
439.92 
437.90 
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435.81 
434.58 
433.76 
432.59 
431.77 
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are also left out in computing the weighted 
average, since in these regions the spread of the 
band due to rotation begins to have an appreci- 
able effect. It may be noticed in Fig. 2 that near 
the bands the theoretical curve lies nearer to the 
extrapolated curve than the experimental points. 
Inclusion of all the points would change the 
averages by only about 1 percent. The agreement 
on the high frequency side of vs is not very good, 
but it must be remembered that in this region the 
dispersion of the monochromator is at a mini- 
mum and that there are difficulties due to two 
overtone bands of methane and the absorption 
bands due to CO, and H;0 in the atmosphere. 


T=} +y?+2) + (Sm/8) (Grr +46?) — 2(m/8) (Gidst 
+ Gags t+ +2 (m/ 8) 


In order to get a picture of what might be 
occurring in the methane molecule during vibra. 
tion, a calculation was made to see how much 
charge would have to be placed on each of the 
hydrogen atoms (with, of course, four times ag 
much on the carbon) in order to produce the 
observed effect on the index of refraction. The 
coordinates chosen for the calculation (shown jy 
Fig. 3) were those given by Dennison and 
Johnston’ except that the x, y, and z were re. 
versed in direction. The expressions for the 
kinetic and potential energies were taken from 
Dennison and Johnston’s paper and are 


where 4=4mM/(4m+M) and m and M are the masses of the hydrogen and carbon particles, 


respectively ; 


Ve +92? +98 +96) +2¢(qigst gigs t+ Gigs +9293 +9295 
9396+ + 2e[ x(q: — G2) 9s) }}. 


If we let e be the charge on each hydrogen atom, the potential energy due to an electric field 8 
parallel to the x axis is 4e8x. The Lagrange equations which involve x are then 


—g2) = 408, (1) 
(Sm/8)G1— (m/8) (Gs+Gst+Gs+Ge) +(m/8)G2+bq: +dg2+ex=0, (2) 
+ — ex =0. (3) 


The remaining equations, those involving y and z, will be of the same form as (1), (2), and (3), except 
that &,=6&,=0, and will clearly be satisfied by gs=qs=qs=Qe=¥=2=0, gitg2=0. Subtracting (3) 


from (2) we obtain 


If the g; are of the form A cos wt, @= —w’g. If 
the applied field is given by &=&) cos wt the 
polarization per molecule due to the forced 
vibrations may then be obtained by solving 

(— pw? +a)x —e(gi— ge) = 4886, (5) 

=0. 

This gives for the polarization per molecule p 
(4e)?6(—m/2)w?+b—d 


= . (6) 
(—pw*+a)[(—m/2)w?+b —d]—2e* 


p=4ex 


Since it is more convenient to have the polariza- 
tion expressed in terms of the resonance fre- 
quencies, Eq. (6) can be transformed as follows: 


(5m/8) (G1 — G2) — (m/8) (G1 — G2) — 92) —d (qi — + 2ex =0. 


(4) 


Let As, A4 be the two values of w* corresponding 
to the natural periods of (5). One has then 


(mp/2) (w? — — As) 
= 


so that 
+(m/2)a] 


Asha= (2/mu)[a(b—d) — 2e*]. 
From the first of these: 
(7) 


7 M. Johnston and D. M. Dennison, Phys. Rev. 48, 879 
(1935). 
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. Equilibrium coordinates for methane. The hydro- 

teat oy ea a, b, c, d. The carbon is at the center of 

the cube and x, y, and z give the displacements of the 
carbon with respect to the center of gravity of the hy- 
ns. The gi--*gs give the increases in length of 


According to (6): 


162°& —w*+(2/m)(b—d) 
Da) 


4ex = 


Substituting from (7) one obtains 
16226 
(w?—Xs)(w? — a) 
Aa) 


wl 


4ex= 


(8) 


ole 


Then if N is the number of molecules per cc 
16e°f (As /(As— Xa) 


wl 


4nNp 
A(n?—1)= =4rN 


+ | 


where the first term gives the contribution of vs, 
the resonance frequency corresponding to Xs, 
and the second term gives the contribution of v,. 
Since the oscillator is isotropic the result is un- 
changed when averaged over all directions of 
polarization of the incident light and all orien- 
tations of the molecule. The above calculation 
has been made according to classical theory on 
the assumption that the carbon and the hydro- 
gens carry fixed “effective charges.’’ Comparison 
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with experiment indicates that this hypothesis 
is too simple, and that the two active normal 
modes of vibration must have different e's. Eq. 
(9) may then be rewritten 


Da) 


— ps? 


A;(n? = 1) = 


and 


A,(n? 1) = 


—K, 


ve 
where (27v)? has been substituted for w*, (2xv3)? 
for As, (2av4)? for X44, and K3, A3(m?—1) and 
A,(n?—1) have the same significance as when 
previously introduced in the interpretation of 
the data. It will be noted that a is the only po- 
tential constant which appears in the above 
expressions. This is reasonable since it is the 
constant connected with the displacement of the 
carbon away from the center of gravity of the 
hydrogens, which displacement can be expected 
to produce electrical dissymmetry in the mole- 
cule. It will be noticed that if e, the constant 
giving the coupling between x and q; and qz, is 
set equal to zero, A; will reduce to the result for 
a simple oscillator of mass u and charge 4e and 
force constant a, while A, will be zero. 

The above expressions were used to calculate 
the values of e for the two active modes of vibra- 
tion. In this calculation the following numerical 
constants were used: yw=5.03X10-** gram; 
a=7.670X10° dynes/cm (Dennison and John- 
ston); m=1.67X10-** gram; A3=32.4X10*8;. 
= 6.04 X 1078; = 9.05 10"* vib. /sec.; = 3.915 
X10" vib./sec.; c=velocity of light=3X10'° 
cm/sec.; N=2.69X10'* molecules/cc under 
standard conditions; and e=4.80X10-" e.s.u. 
=electronic charge. From this calculation we 
obtained for e3, the charge required on each H 
atom for the vibration v3, 0.562X10-° e.s.u. or 
0.117 e where e is the electronic charge. For e, 
we obtain 0.277 X10~-"° e.s.u. or 0.0578 e. These 
results should not, of course, be interpreted as 
meaning that the H atoms in methane actually 
carry a charge of 0.562 X 10-" e.s.u. for one mode 
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Fic. 4. Change of the dipole moment of the CH group 
as a function of CH distance. Curves a and } give two 
ibilities; the experiments yield only the value of » and 

its derivative at the equilibrium position. 


of vibration and 0.277 X 10-"° for the other. They 
merely mean that if that were the case the 
observed values of the index of refraction would 
be obtained. It is highly probable that the charge 
on the H atoms is not constant during a vibra- 
tion. The refraction measurements which we 
have made cannot distinguish between the 
various ways in which the dipole moment of the 
molecule can vary during a vibration. 

In order to gain further insight into what may 
actually be happening in the methane molecule, 
it is interesting to examine in closer detail the 
actual motions of the nuclei for each of the two 
active frequencies v3 and v4. For this calculation 
it is desirable to know the coefficients involved in 
the transformation to normal coordinates, and 
these are given in the appendix. A semi-classical 
calculation, which involves setting the potential 
energy at maximum displacement equal to hy 
gives the following results: For v3, if only the 
one of the 3 equal frequency modes which in- 
volves displacement of the C in the X direction 
is excited, the carbon moves 0.0402 X 10-8 cm in 
the X direction with respect to the center of 
gravity of the hydrogens, g; decreases by 
0.0906 X10-* cm, and ge increases by 0.0906A. 
Expressed in another way, the CH distance 
varies between the limits 1.045A and 1.17A while 
the angle between the CH bonds changes only 
from 109°0’ to 110°4’ or each CH bond vibrates 
only about 3° on each side of its equilibrium 
position. For v4, the C moves 0.751A in the +X 
direction while q; increases by 0.962A and q: 
decreases by 0.962A. In more pictorial language, 
the CH bond vibrates between 100°14’ and 
118°42’ or 9° on each side of its equilibrium 
position, while the CH distance changes from 
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1.107A to 1.118A or only about 0.005A on each 
side of equilibrium. 

One interesting possible interpretation then 
presents itself. In the vibration corresponding to 
vs the charges on the hydrogens may remaig 
nearly constant since the CH distance remains gp 
nearly constant and the hydrogens are go far 


_apart. (This, of course, neglects any directiongy 


effect of the carbon on the distribution o 
charge.) The e, above then may be interpreteg 
as actually the charge on the hydrogen and the 
electric moment of the CH bond may be com. 
puted. Since the equilibrium distance between the 
C and the H is 1.11A, the moment of the bong 
is or 0.307X10-8 
This will of course be equal to the dipole moment 
of the CH; group in methane. For 1s, since the 
angles remain so nearly constant, e3 may be 
interpreted as the rate of change of the dipole 
moment of the CH bond with respect to the CH 
distance. If this interpretation is correct, the 
dipole moment of the CH bond may vary with 
CH distance somewhat as shown in Fig. 4. Ip 
this connection it is interesting to note that from 
measurements of the intensity of absorption jn 
overtone bands of several halogen derivatives of 
hydrocarbons, Timm and Mecke® estimate the 
moment of the CH bond to be between 0.3 x 10-# 
and e.s.u. 

One further point which must be mentioned 
is the extrapolation to infinite wave-length, where 
N»*—1 should be equal to e—1. If we take 
435.910 as n,,—1 from measurements in the 
visible and ultraviolet regions of the spectrum 
and add 1.81 X10~* and 4.41 X10~-* as the con- 
tributions of vs and v4, we obtain n,—1=44211 
X10-* and m,?—1=884.4X10-*. Sanger ob 


TABLE A. Coefficients for transformation to normal 
coordinates for methane. 


&: & 


qi (2/3m)} a dD, 
q2 (2/3m)* a —D, -E, 
qs (2/3m)* b—a E; 
(2/3m)t b—a -D, -E; 
(2/3m)+ —b D; E; 
(2/3m)} —b -E; 
x 0 0 2Dye/(urs—a) 
y 0 0 2Dye/(uds—a) 2Exe/ (wa) 
0 0 2Dye/(uds—a) 


* B. Timm and R. Mecke, Zeits. f. Physik 98, 363 (1936). 
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tained 960X10~* for «—1. The contribution of 
the infra-red bands thus reduces the discrpeancy 
between ”.?—1 optical and «—1 by only about 
12.4X10-* or 16 percent of the difference. The 
agreement between various determinations of 
fte-1 makes it extremely improbable that 
changes in its value can reduce the discrepancy 
by more than 10 percent. The remainder must 
then be due to error in measurements of «. 

It is a pleasure to acknowledge the help re- 
ceived from Professor G. Breit in the interpreta- 
tion of the results, and to thank the Wisconsin 
Alumni Research Foundation for a grant in aid. 


APPENDIX 


Table A gives the coefficients for the trans- 
formation to normal coordinates for methane. 
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The D,23, E:,23, @ and } are constants which 
enter because of the degeneracy of the corre- 
sponding modes of vibration. 

If T is to have the form 


8? 


there are the following restrictions on the 
constants: 


a?+b?—ab= im, 


2e 2 
D?+D?+D3= utm 
wA3—a 


2e 2 
( ) utm 
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The Diffraction of Radio Ranges by Hills 


WituiaM R. HASELTINE 
Department of Physics, University of California, Berkeley, California 
(Received January 2, 1940) 


The scattering of electromagnetic waves by hemispherical 
obstacles on a conducting plane has been investigated. The 
results have been applied to the problem of the splitting 
of radio range beacons when the range passes over hilly 
country. The calculations refer particularly to long wave- 
lengths (about 1000 meters) such as used in the Oakland 
range. It is found that even hills of moderate size can cause 


N the standard system of radio ranges for the 

guidance of aircraft the radiating antenna 
array sends out two separate signals on the same 
frequency. One signal, modulated with the code 
letter N (—-) is strong in two opposite quad- 
rants; the other, modulated with A (-—) is 
strong in the remaining quadrants. Along four 
lines passing through the station, the A and N 
signals are of equal intensity, and blend into a 
continuous sound. These lines form the courses 
or ‘‘beams”’ defined by the system. The radiation 
fields are vertically polarized to eliminate some 
of the troubles that occurred when earlier types 
of antenna systems were used. 

In level country these range courses are reason- 
ably well behaved, but over mountainous or even 


wide spatial fluctuations in the received signal strength. 
The results of numerical calculation for one simple model 
of the terrain are reported. The beam is found to split into 
several parallel courses, giving a pattern of the same type 
as is observed in the field. It is concluded that this phe- 
nomenon is a typical diffraction effect. 


moderately hilly terrain curious effects often 
arise. In particular, instead of a single course of 
equal signal strength, there may appear a com- 
plicated pattern of several parallel courses spread 
over an angular distance of up to twenty degrees. 
This effect is very troublesome and dangerous 
to air traffic. 

Early in 1939, field investigations of these 
phenomena were carried out on the range at” 
Oakland, California, by Captain E. C. Dyer, 
U.S.M.C., and Lieutenant D. E. Wait, U.S.N. 
It was found that the seaward course was badly 
split and that this seemed to be related to the 
hills of the San Francisco peninsula and es- 
pecially Mt. Tamalpais in Marin County north 
of the Gate. They requested the author to make 
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Fic. 1. Polar plot of un- me 
perturbed field strength. 
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a theoretical investigation of this problem to 
ascertain the probable cause of the multiple 
beams observed. 

Because the great complexity of the terrain 
prohibits study of the actual conditions, a simpli- 
fied model of one or two hills on an infinite plane 
is considered. For simplicity, the field strength 
of the A signal is taken to vary as the cosine of 
the azimuth measured from an arbitrary zero, 
and the N signal as the sine of the same angle. 

The electric and magnetic field strengths 
radiated must satisfy Maxwell’s equations. If 
E=Ee-*', H=He-‘*', these reduce to 


VE+FRE=0, VH+kH=0, w=2nc/d (1) 


in which c is the velocity of light, and \ the 
wave-length in free space. In free space (or the 
atmosphere), 
Inside matter, 
Rk? = (w?/c*) (eu 


where ¢ is the dielectric constant, » the permea- 
bility, and o the conductivity in electrostatic 
units. This constant will be written k’ to differ- 
entiate it from the value for free space. The 
permeability u is supposed to be equal to one. 
Consider first the solution over a plane earth 
with no hills present. The effect of the earth on 
the wave field is considerably simplified by the 
long wave-length involved (for the Oakland 
range, \= 1267 meters). Sommerfeld! has shown 
that the earth can be considered perfectly con- 
ducting so long as (k?/k”)kR<1 (R is the dis- 
tance from the source). With 108 for average 
land, \=1267 meters, and R=15 miles, this 
quantity has the value 1/10. For sea water, 
with its much greater conductivity, it is still 
smaller. Treating the earth as perfectly conduct- 
ing, as is thus allowed here, is equivalent to 


1A. Sommerfeld, Ann. d. Physik 81, 1135 (1926). 
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requiring that the electric field at the surface be 
vertical. 

Thus the solutions sought for the primary 
undisturbed wave are 


eikR cos 
E,=E sin? x . 


x is the polar distance; y the azimuth; the factor 
cos ¥ applies to the N signal; and the sin y ty 
the A signal. E, is the electric field strength in the 
direction of increasing x; the other components 
are zero. Fig. 1 shows a polar plot of the field jn 
a horizontal plane. 

When a hill is present, this wave will bp 
scattered by it, and the amplitudes of the pri. 
mary and the scattered waves will add with 
differing phase relationships at different points 
An interference pattern is set up. This problem 
is very similar to that found in the quantum. 
mechanical treatment of atomic collision and 
scattering processes. It is almost identical with 
that of the scattering of electron waves by q 
central field of force, except that the electronic 
wave function is a scalar, while the wave functiog 
here is a vector. There are two general methods 
of attack. The first is that of small perturbations 
and is equivalent to the Born approximation, 
The second is that of partial waves. This tech. 
nique was also used by Lord Rayleigh to study 
the scattering of sound waves. The first method 
can be used if the scattered field is small com- 
pared with the primary field even at the hill (as 
it will be when the conductivity is sufficiently 
small). In this case the following approximation 
is valid: 

From the primary field strength calculate the 
currents set up in the hill, I=cE. Then calculate 
the field radiated by these currents(see Hansen’). 
This is the first approximation to the scattered 
wave. This method was tried, but in the case of 
interest the conductivity was so high that the 
approximation broke down. 

In the second method one attempts to finda 
rigorous solution of the wave equation having 
two characteristics. It must satisfy the boundary 
conditions at the surface of the hill (tangential 
components of E and H and normal components 
of D and B continuous across surface), and it 


2 W. W. Hansen, Physics 7, 460 (1936). 
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must represent the sum of the primary signal and 
an outgoing scattered wave. That is, it must 


the form: 
E=E,+E,, (3) 


where Ep is the primary wave of form (2), and 
E, is one which behaves at large distances from 
the hill as e**/r, 7 being the distance from the 
center of the hill. This will be practicable only 
in very simple cases. Here we will investigate the 
case of a hemispherical boss of radius ro. A 
solution of the desired type can then be found in 
series form. 

It has been shown* * that the functions 
6, ¢) and A’,, in 6, (s=2, 3; l=0, 
1,2, -::) defined by Hansen form a complete 
set for the description of any radiative solution 
of Eq. (1). These functions are defined as 


follows: let 
| 
h = 
As X(ré), 
These A and A’ functions are the convenient 
vector solutions of the wave equation for cases 
of spherical symmetry. The A’ functions are 
generated in the same way as A except that the 
Bessel function Ji,,; in is replaced by HZ! 


the Hankel function. of the first kind. 
Any solution can then be written as 


E= as, 1, wAs, 1, nb, 1, 


sin 


(5) 


_ Now the A’ functions have a singularity at the 
origin and behave at large distances like e“*/r; 
J 
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they represent outgoing waves. The A functions 
behave at large distances like (sin kr) /r. Further- 
more, any wave field not having a singularity at 
the origin can be expanded in terms of the A 
functions alone. Thus in (5) the A terms can be 
considered as an incident wave and the A’ terms 
the corresponding scattered wave. Then if the 
origin is chosen at the center of the scattering 
hill, this is of form (3) where 


For any form of primary wave Ep, i.e., any 
given values of a,,:,,, the boundary conditions at 
the surface of the hill determine the 6,,;,. In 
fact, the requirement that the appropriate fields 
be continuous across this boundary gives a rela- 
tionship between any a,,;,, and the 6,., with 
the same indices alone. This separability of the 
components is the reason for using this system of 
functions. If, then, the expansion of Eo is known 
in terms of the A functions, the whole field is 
known. The amplitude of the scattered wave 
made up of the A’ functions may then be added 
to that of the primary signal at any desired point. 
The square of the absolute magnitude of this 
sum is then proportional to the signal intensity 
received. This procedure can be carried out for 
both A and N signals, and the relative intensities 
compared. 

The solutions inside the hill have the form 
where 


n=k’V al 1/ke')V I, ny 
Gs, 


t’ is the same as £ with k replaced by k’ =ke’!. 
The equations for 6 are found to be 


8,1, 


(kro)! (kr’)® dro (k'ry)3 


d (kr)! ro) 


ere)! 


If the hill is assumed to be perfectly conduct- 
ing, the corresponding equations are 


*W. W. Hansen, Phys. Rev. 47, 129 (1935). 
*W. W. Hansen, Proc. I. R. E. 24, 1594 (1936). 


drL (k’r)! Ino 


d 
id 
d (kr)! ro 


Qe, 1, 1, 144(Rro) =0, 


d (kr)! ro 


dfrll™ 
+s, | 


d (kr)! 
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_ Fic. 2, Relative po- 
sitions of station, scat- 
tering hills, and ranges. 


Numerical calculation shows that with reason- 
able values for ¢’, the magnitude of the reflection 
coefficients b/a are not very different from those 
for a perfectly conducting sphere, though the 
phases may be different. In the following calcula- 
tions, the hill was taken to be a perfect reflector. 

If the radius (height) of the hill is small com- 
pared to its distance from the station, the pri- 
mary wave front is only very slightly curved, and 
for purposes of computation can be assumed to 
be plane over the extent of the hill 


cos e**®o Cos 
R sin Ro sin 


=Eyeitr 


where Ro, xo(=7/2), Yo are the coordinates of 
the center of the hill, and 6’ is measured as polar 
distance from an axis through station and hill. 
Furthermore, the electric field is entirely vertical. 

This plane wave can be expanded in a well- 
known form 


eikr cos => 
XJ m(cos 6’). (10) 


Each term of this series can be in turn expanded 
in terms of the A.,;,. functions, bearing in mind 
that these are vectors and that the primary field 
is vertical. The result gives the a,,;,, quantities. 


ka, Ln = (0)h?,, ny 
kas, 1, n= — [2 — ](x/2)! 
The series expansion of the scattered wave can 
always be found in principle; but if the radius of 
the scattering body is larger than a wave-length, 
an inordinate number of terms must be used. 


Here again the long wave-length used is a help. 
For kro=2rro/d= 1.6, about eight’ terms must be 
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(11) 


used. For kro=4, at least 24 must be carried 
This latter figure corresponds roughly to the las 
of Mt. Tamalpais (ro>=810 meters). 

The necessary numerical calculations are 
rather laborious, and have been carried out with 
considerable assistance from Mr. Chaim Rich. 
man working on N.Y.A. time. 

The case of a small hill, kro=1.6, and of a 
single large one, kro=4, showed pronounced 
fluctuations of the signal intensity from the 
normal. However, there was no evidence of course 
splitting found. This is undoubtedly due to the 
simplicity and symmetry of arrangement of 
scattering material. 

Next the scattering by two hills was investi. 
gated. The results here reported are for two hills, 
each of radius kro=4 placed on a circle distant 
about 12 miles (kR»=100) from the transmitting 
station. One hill was at an azimuth of 30°, the 
other at 60°. The field was calculated for both 4 
and WN signals along a circle for which kR=110. 
The geometrical relations are shown in Fig, 2, 
The work was carried out only at ground level 
(x=0= 7/2), because at any other elevation the 
labor would have been multiplied at least three. 
fold. The effect of one hill on the field of the other 
was neglected as a higher order perturbation. 

By choosing properly the relative magnitudes 
of the A and WN signals, the beam may be made 
to fall anywhere in the quadrant. When it lies 
at 45°, directly between the hills, there can be 
no splitting or parallel courses, because of sym- 
metry, even though the fields may be violently 
distorted. When, on the other hand, the range 


Fic. 3. Plot of signal intensities when desired range lies at 
42 degrees. 


he multiple ranges are clearly seen. 
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lies only slightly off this position, the situation 
arises shown in Fig. 3. 

The signal intensities are plotted as calculated 
for points along the circle of Fig. 2 from 30° to 
60°. The intensities that would be observed with 
no scatterers present are represented by the 
smooth, almost straight dotted curves, the NV 
signal decreasing to the right and the A to the 
left. Their crossing point defines the desired 
range here at about 42°. This is indicated by the 
tall arrow below. The magnitude of the dis- 
turbance caused by the hills can be seen from the 
solid and dashed curves, which represent the 
intensities to be expected with the hills interfer- 
ing. The solid curve is for signal N and the dashed 
for signal A. It is seen that they cross in three 
places, marked below 0, 1, 2, thus giving rise to 
three more or less parallel courses. The cross- 
hatched rectangles along the lower axis indicate 
the regions in which the A signal is predominant. 

Figure 4 is a similar plot for the case in which 
the desired course is about 35°. Here there are 
three crossing points and one point of tangency, 
or near-tangency, giving four courses. Again the 
large arrow indicates the desired course, and the 
actual beams are numbered 0, 1, 2, 3, the cross- 
hatched rectangles indicating predominance of 
the A signal. The small arrow at 1 indicates the 
course set by the point of tangency. 

In Fig. 2 the pattern is represented by small 
solid lines attached to the large circle. These lines 
represent the positions along the circle at which 
an on-course signal would be heard. Those that 
extend outward refer to a beam set for 42°, and 
those extending inward refer to one set at 35°. 
In each case the longer line indicates the position 
of the main range. The shorter lines mark sub- 
sidiary courses. 

On examination of these curves there are 
present all the types of beam splitting or parallel 
courses actually observed. The predominant 
signal may be the same on both sides of the 
course, as at 1 in Fig. 4; it may be different in 
either the normal order (as at 0, Fig. 3) or the 
reversed order (1, Fig. 4), and none of the courses 
need lie along that expected. 

The pattern observed at a higher altitude 
would undoubtedly be different. But it would 
show the same type of phenomena. 


40 50 60° 


Fic. 4. Plot of signal intensities when desired ran 
fies at 35 degrees. - 


The number of peaks observed in the signal 
strength between 30° and 60° corresponds to the 
fact that along this circle the wave-length repre- 
sents about 4°. The relative phase with which 
the scattered waves add to the primary thus 
shifts with this period. At a greater distance from 
the station, the wave-length would correspond to 
a smaller angle, and more peaks would be found. 
The pattern of parallel courses would thus 
change entirely in going from kR=110 to, say, 
kR=150. These changes of pattern with both 
altitude and distance have been observed in 
the field. 

In any actual geographical location where these 
phenomena occur, the scattering material will 
be present in much greater quantity, and will be 
distributed much more irregularly than that of 
the model used in this paper. It is therefore to be 
expected that the splitting observed will be more 
pronounced and more erratic. Further, the extra 
courses must vary with ground and weather con- 
ditions. A change in soil conductivity will intro- 
duce phase shifts of different amounts in all 


scattered waves and definitely affect the pattern. 


The author is grateful to Captain E. C. Dyer 
for details on his unpublished field work, and to 
Professor J. R. Oppenheimer for several very 
helpful talks. Without the aid of Mr. Richman 
the work would have taken much longer, and 
the author expresses thanks to him and to 
the N.Y.A. 
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A Theory of Spark Discharge 
J. M. Meek 
Department of Physics, University of California, Berkeley, California 
(Received February 13, 1940) 


The breakdown of a uniform field is considered to occur 
by the transition of an electron avalanche proceeding from 
cathode to anode into a self-propagating streamer, which 
develops from anode to cathode to form a conducting 
filament between the electrodes. A criterion is put forward 
for such a transition, viz., a streamer will develop when the 
radial field about the positive space charge in an electron 
avalanche attains a value of the order of the external 
applied field. For then photoelectrons in the immediate 
vicinity of the avalanche will be drawn into the stem of the 
avalanche and will give rise to a conducting filament of 
plasma, and a self-propagating streamer proceeds towards 
the cathode. The theory thus depends on ionization by 
electrons and photo-ionization in the gasand dispenses with 
the classical assumption of ionization by positive ions in 
the gas or secondary actions at the cathode. On this basis 


an equation for breakdown is developed, and reference to 
a/p—X/p curves enables the potential required for 
down to be determined. Satisfactory agreement between 
calculation and experiment is found in air for values of 
pressure times gap length down to p5~100 mm Hg Xem, 
The theory does not conform absolutely to Paschen’s lay 
but the deviations are within the present day margins of 
experimental error. For lower values of p35 the deviation 
between calculation and experiment is explained by the 
fact that the density of photo-ionization becomes 

and the secondary mechanism, 7, is observed in this reg; 
so that classical theory applies. The theory is indicated ty 
be consistent with all the requirements so far established 
in connection with sparks for large pé up to the lightning 
discharge. 


INTRODUCTION 


N the early 1900’s J. S. Townsend! evolved his 
now classical theory of the mechanism of 
spark discharge on the basis of measurements 
made at high values of X/p, the ratio of field 
strength to pressure, and low values of pé, the 
product of pressure times gap-length. His theory 
has been extrapolated to explain the mechanism 
of sparking under ordinary conditions for longer 
gap-lengths at atmospheric pressure and above. 
It has formed the guiding principle of all con- 
siderations of the sparking mechanism for nearly 
40 years. Many researches have indicated that 
for values of 5 near and somewhat above the 
minimum sparking potential the theory was 
adequate and relatively accurate. The theory in 
this region enables calculation of the sparking 
potentials to be made with adequate accuracy, 
and recently it has been shown by Schade? to be 
consistent with time-lag studies. 


1J. S. Townsend, Nature 62, 340 (1900); Phil. Mag. 1, 


198 (1901). 

2P. O. Pedersen, Ann. d. Physik 71, 371 (1923); J. W. 
Beams, J. Frank. Inst. 206, 809 (1928); J. J. Torok, Trans. 
A. 1. E. E. 47, 177 (1928); 48, 46 (1930); Buraway, Arch. f. 
Elek. 16, 14 (1926); I. Tamm, Arch. f. Elek. 19, 235 (1928); 
W. “7 Arch. f. Elek. 20, 99 (1928); F. G. Dunning- 
ton, Phys. Rev. 38, 1535 (1931); H. J. White, Phys. Rev. 
46, 99 (1934); A. Tilles, Phys. Rev. 46, 1015 (1934); R. R. 
Wilson, Phys. Rev. 49, 1082 (1936); M. Newman, Phys. 
Rey. 52, 652 (1937); R. Strigel, Wiss. Veréffentl. Siemens 
Werken 15, 1 (1936). 


In recent years, however, it has become jp. 
creasingly apparent that the theory is entirely 
inapplicable when extrapolated to sparks for pj 
greater than about 200 mm Hg Xcm in air, and 
analogous values in other gases. The reasons for 
this conclusion may be stated briefly as follows: 

1. The formative time lags of sparks at gj. 
mospheric pressure, i.e., the time required for the 
breakdown mechanism to materialize, have been 
measured by several diverse methods, and all 
lead to time intervals of the order of 107 sec, 
or even less, for a one-cm gap at small over. 
voltages.? These time intervals are far below the 
microseconds required by positive ion movement 
as shown by Schade’ and others.‘ 

2. At atmospheric pressure the sparking poten- 
tial has been found to be independent of cathode 


material to a very high degree,® while the theory’ 


above, even though somewhat insensitive to the 
value of y at higher pressures, requires a definite 


‘ dependence of the order of 10 percent on the 


value of y.° The values of pé at which this de 


3R. Schade, Zeits. f. Physik 104, 487 (1937). 

41. Tamm, Arch. f. Elek. 19, 235 (1928); W. Rom 
Arch. f. Elek. 20, 99 (1928); W. Rogowski and I. Tama, 
Arch. f. Elek. 20, 625 (1928); A. Tilles, Phys. Rev. 


1015 (1934). 
5L. B. Loeb, Fundamental Processes of Electrical Die 


charge in Gases (John Wiley, 1939), p. 415. 
® Reference 5, p. 416 et seg. 
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ce begins rapidly to vanish occur again 

3. In longer sparks, in lightning discharge, and 
in positive point corona at atmospheric pressure, 
it is virtually certain that the character and 
action of the cathode is without any influence 
whatsoever on the process. 

4. The observed appearance of sparks in 
measurements of Townsend’s coefficient, a, at 
higher pressures and adequate x with no measur- 
able value of y indicates in another way that the 
theory based on a ¥ is no longer applicable at 
higher values of pressure and gap length.’ 

5. Observations of actual sparks in short time 
intervals made visible by Kerr cell-shutter 
studies® and by cloud-track pictures of spark 
paths after short time intervals,® as well as the 
usual visual observations of sparks, indicates that 
breakdown at higher pressures occurs along very 
narrow channels or filaments. This is incom- 
patible with the transient rearrangement of con- 
duction currents in a gas caused by accumula- 
tions and movements of positive ions or suc- 
cessive waves of electron ionization when 
ye*>1. This discrepancy is made still more 
glaring by the zigzag nature of the spark filament 
and the frequent branching observed in longer 
sparks. 

6. While only one electron is required to 
initiate a spark on the Townsend theory as inter- 
preted by Holst and Oosterhuis” it is clear that 
the greater the number of photons incident on 
the cathode by external irradiation the more 
rapidly should the space charge accumulations 
occur, so that formative time lags should on that 
theory decrease as the illumination intensity of 
the cathode increases. This has been shown to be 
true by Schade for low values of p5. However, at 


_high pé observation shows no appreciable change 


in the formative time lag and of the sparking 
potential with illumination over a considerable 
range of intensity."' When the illumination 


'F.H. Sanders, Phys. Rev. 44, 1020 (1933); D. Q. Posin, 
Phys. Rev. 50, 650 (1936); W. E. Bowls, Phys. Rev. 53, 
293 (1938). 

*L. von Hamos, Ann. d. Physik 7, 857 (1930); F. G. 
Dunnington, Phys. Rev. 38, 1535 (1931); H. J. White, 
Phys. Rev. 46, 99 (1934). 

*E. Flegler and H. Raether, Zeits. f. Physik 99, 635 
(1936); 103, 315 (1936); H. Raether, Zeits. f. Physik 96, 
567 (1936); 110, 611 (1938); 112, 464 (1939). 

”G. Holst and E. Oosterhuis, Phil. Mag. 46, 1117 (1923). 

"H. J. White, Phys. Rev. 48, 113 (1935); R. R. Wilson, 
Phys. Rev. 50, 1082 (1936). 


intensity is increased by a factor of 105 beyond 
that giving the 10’ electrons/sec./cm* generally 
used, both formative time lag and sparking po- 
tential are lowered.” In this case the potential is 
lowered at a maximum by no more than 10 
percent, and the formative time lag appears to 
be doubled when the intensity of illumination is 
decreased by a factor of 500. 


THEORY 


The above-mentioned objections to the clas- 
sical theory make it imperative to formulate a 
new mechanism which does not so conflict with 
experimental observations. The importance of 
streamer formation in the mechanism of spark 
breakdown has been discussed in a recent article 
by Loeb and Kip'*'** and more extensively by 
Loeb in his book Fundamental Processes of Elec- 
trical Discharge in Gases. The breakdown process 
is there pictured in qualitative fashion, but no 
quantitative criterion is given for streamer 
formation and breakdown. However, it is shown 
that under sparking conditions the positive space 
charge in an electron avalanche which is origi- 
nated at the cathode produces an electric field of 
the same order as the external applied field when 
the avalanche reaches the anode. This fact was 
also indicated previously by Slepian,’* who then 
postulated the development of thermal ionization 
in the avalanche. It is further stated by Loeb that 
the axial field distortion produced by the positive 
space charge inhibits the advance of the electron 
avalanche towards the anode, but favors the 
development of a positive streamer towards the 
cathode. An explanation of the mechanism of 
streamer formation and mid-gap streamers is 
also given, though the fundamental criterion for 
their appearance is lacking. However, it is indi- 
cated that breakdown will not occur until a 
conducting filament exists across the gap. The 
single electron avalanches observed by Raether'® 
just below breakdown do not form such a fila- 


"H. J. White, Phys. Rev. 48, 113 (1935). 

3 L. B, Loeband A. F. Kip, J. App. Phys. 10, 142 (1939). 

3e At the request of essor b, attention is called 
to page 157 of this article where a careless blunder is made 
in the calculation of the field strength produced by space 
charge. However, the order of anaaied of the calculated 


values is not affected by correction, and the argument 
given in the article is by no means invalidated thereby. 

Slepian, Electrical World 91, 761 (1928). 

* H. Raether, Zeits. f. Physik 107, 91 (1937). 
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ment but produce a charge concentration at the 
anode so that positive streamer propagation to 
the cathode is needed. 

In the above-mentioned theoretical discus- 
sions, consideration appears to have been given 
principally to the field distortion produced by the 
positive space charge in the direction of the axis 
of the avalanche, i.e., in the direction of the 
external applied field. Clearly, however, the 
positive space charge will also produce a field 
distortion in the direction radial to the axis of 
the avalanche, i.e., perpendicular to the external 
field. Consideration of this fact led the writer to 
the following criterion for the transition of an 
electron avalanche to a conducting streamer, a 
criterion which makes quantitative calculation 
possible. A streamer will develop when the radial 
field about the positive space charge in an electron 
avalanche attains a value of the order of the external 
applied field. Then photoelectrons produced in 
the gas in the immediate vicinity of the avalanche 
will not only be accelerated in the direction of the 
external applied field but will also be drawn into 
the stem of the avalanche. (It has been shown 
by Cravath,!* Dechene,” Loeb and Kip" that 
streamer propagation depends on photo-ioniza- 
tion in the gas.) In this manner the positive ion 
space charge attracts photoelectrons from the 
surrounding space in sufficient numbers to 
originate a self-propagating streamer. For lower 
values of external applied field, where the multi- 
plication of ions in the avalanche is insufficient 
to produce a radial field of the same order of 
magnitude as the external field, photoelectrons 
produced in the gas proceed to the anode and are 
not deflected to the main avalanche. 

It is clear that the exact equality of the radial 
field and the external applied field need not be 
insisted upon, and it is possible that the requisite 
criterion could be satisfied when the radial field 
is only 50 percent of the applied field. However, 
in the calculations the two fields have been set 
equal, and the application of this condition at 
once leads to equations for field strength for spark 
breakdown that closely approximate experi- 
mental observations and clarifies difficulties in 
the qualitative picture given by Loeb. The 
criterion not only enables quantitative calcula- 


16 A, M. Cravath, Phys. Rev. 47, 254 (1935). 
17C, Dechene, J. de phys. et rad. 7, 533 (1935). 


tions of breakdown to be made, but it also ase; 
in the interpretation of a number of other spark. 
ing phenomena. It thus seems of interest to 
present the results of the analysis and its appli- 
cations, together with its experimental justi. 
fication. 


DERIVATION OF EQUATION FOR BREAKDowy 


Because of the cumulative character of the 
ionization in an electron avalanche and the 
relative immobility of positive ions as compared 
with electrons, the bulk of the positive ions left 
behind by the advancing electron swarm are 
concentrated in a region of several ionizing free 
paths behind the tip of the avalanche. The region 
is also bounded by the periphery of the avalanche, 
which is conditioned by electron diffusion. Since 
the estimation of the field produced by such a 
distribution is difficult it will be assumed for 
convenience of calculation that the positive space 
charge is concentrated in a spherical volume of 
radius, r, equal to that of the avalanche. While 
the volume may be more nearly in the form of an 
oblate or prolate spheroid, the assumption of a 
perfect sphere will not affect the calculated field 
in order of magnitude. 

The field, X:, at the surface of a spherical 
volume of uniformly distributed space charge is 
given by 

Q 

1=— €, 

where N is the number of ions per cm? and ¢ is 
the charge on an electron. To find N, one recalls 
that owing to cumulative ionization the number, 
n, of positive ions formed when the avalanche has 
progressed a distance, x, is m=e**, where a is the 
Townsend coefficient for ionization by electrons. 
In a distance dx the number of ions created is 
dn=ae**dx. These ions are assumed to be con- 
tained in a cylindrical volume of radius, r, and of 
length, dx, so that 

N=ae%/nxr’. 

According to Raether"® the radius, 7, is related to 
the distance, x, of avalanche travel by means of 
the diffusion equation r?=2Dt=2Dx/v, where 9 


is the velocity of advance of the avalanche and 
has been measured by Raether'® and White.” 


18H. J. White, Phys. Rev. 46, 99 (1934). 


a2 & = 


| 

T 

a 
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ce v= KX, where K is the mobility and X is 
the applied field, we have 


4 ae**e 4 


3 r ~ 3 (2Dx/KX)* 


Xi= 


Both D and K are functions of X/p and are 
approximately related by the expression 


D_ P of 


K Ne 


where P is atmospheric pressure, Ne is the 
Faraday constant, c is the speed of thermal agita- 
tion of electrons (1.2 X10" cm per sec. at 20°C), 
and ¢; is the square root of the mean-squared 
velocity of the electron in the applied field. 
Townsend” has shown that c; is considerably 
larger than c. According to Compton” the value 


of ¢; is given by 
e XX 


c?=1. 
m 


where \ is the mean free path of the electron, X 
is the applied electric field, m is the mass of an 
electron, p is the gas pressure, and f is a factor 

which depends on the fractional loss of energy 
per collision. If A» is the mean free path of an 
electron at a pressure of 760 mm Hg, then the 
value of c, for a pressure, p, is given by 


e X Xo 
— —. 
mp /f 


The values of \» and f obviously depend on X/p 
and are not known for the high values of X/p 
which are necessary for spark breakdown. But 
inasmuch as the theory will be found insensitive 


- to slight errors in the value of \o/4/f it is con- 


TABLE I. Breakdown potentials corresponding to different 
values of pandé. 


BREAKDOWN POTENTIAL (KV) 


MM XCM $=0.1 5=1.0 =10.0 
7600 245 248 249 

760 31.5 32.2 32.9 

380 18.0 18.6 18.9 


wy. S. Townsend, Electricity in Gases (Oxford Uni- 
"oy Press, 1914), p. 122ff. 
= 0 (1930). Compton and I. Langmuir, Rev. Mod. Phys. 2, 
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sidered reasonable to calculate its value from 
known data based on experiments with fields 
close to those which produce sparking. Raether"® 
found that the speed of advance of an electron 
avalanche was 1.2510" cm per sec. in a field 
X/p=41 volts per cm per mm. Substitution of 
such values in the well-known Langevin equation 
for electron mobility, viz.; v=0.815 AX/me,, 
enables one to calculate that do/f=5.7 
The mean free path Xo is taken as 3.6 X 10-5 from 
curves given by Brose and Saayman™ and 
corresponds to an electron energy of three 
volts. Whence we have f=0.025 and c,=2.01 
X107(X/p)!. Then we have 
4 ae**¢ 
xX 


3 (2Pc;2x/Nec?X)! 


(1) 
=5.28X 10-7 


volts/cm. 
x/p 


Now for x= 6, the gap length for a plane parallel 
gap, application of the criterion for streamer 
formation, that X,=X, gives 


which may be written in logarithmic form as 


a a xX 
—pi+log.—=14.46+log. ——} log. pi+log, 6. 
p p (2) 


Reference to curves of the Townsend type which 
relate a/p with X/p now enables the value of 
X/p required for sparking to be calculated for 
given 6 and p. In the case of a one-cm gap in 
air at atmospheric pressure the breakdown 
equation is satisfied for X/p=42.4 (a=18.6), 
i.e., X =32.2 kilovolts per cm,* which is in close 
agreement with the observed value of 31.5 kilo- 
volts per cm.” 


COMPARISON BETWEEN CALCULATED AND 
MEASURED VALUES OF BREAKDOWN 
POTENTIAL 


The form of the breakdown Eq. (2) shows that 
it is in agreement with Paschen’s law other than 


1H. L. Brose and E. H. Saayman, Ann. d. Physik 5, 
797 (1930). 

*If we set the criterion for streamer formation as 
X,=0.2X instead of X,=X, breakdown field strength 
evaluted is 31.8 kilovolts per 

2 W. Spath, Arch. f. Elek. 12. 331 (1923). 
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Fic. 1. Experimental V x 
and calculated curves ‘+ 
to show relationship 4 
between breakdown { 
voltage and pé. 
on 


for the term, log, 5. This term has little effect, 
however, on the calculations, as shown by Table I 
for breakdown potential corresponding to dif- 
ferent values of » and 6. The deviations, for the 
same value of »6, are within the present margins 
of experimental error in measurement, and it 
may be that for more accurate measurements 
such differences will be observed. For it should 
be pointed out that the derivation of Paschen’s 
law is based on theoretical considerations 
analogous to those of Townsend. Such consider- 
ations make the mechanism of spark breakdown 
dependent on the total number of ions formed in 
the gap and, since the multiplicative factors in 
ion production depend on X/p, Paschen’s law 
results. However, on the new theory, it is seen 
that the breakdown mechanism is not dependent 
on the total number of ions formed but depends 
on the ion density in an electron avalanche and 
the resultant space-charge field. In consequence, 
we would not expect that Paschen’s law would 
be strictly obeyed. In fact, it has previously been 
shown”? that with space charge distortion 
Paschen’s law is only approximately correct, and 
that the deviation in breakdown potential is of 
the same order of magnitude as that derived 
here. So that unless the precision of experimental 
determinations of sparking potentials in the 
future indicates Paschen’s law to hold beyond the 
prediction of the present theory, the deviations 
of the latter from Paschen’s law are not thought 
to invalidate the theory in any way. 

Both experimental and calculated curves to 
show the relationship between the voltage re- 
quired to cause breakdown and the product, pé6, 
of pressure times gap-length are plotted on 
logarithmic scale in Fig. 1. The calculated curve 
is based on values for the primary ionization 


%R, N. Varney, H. J. White, L. B. Loeb and D. Q. 
Posin, Phys. Rev. 48, 818 (1935). 


coefficient, a, in air as given by Sanders% Sor 
X/p up to 100; for higher X/p the values of is 
are taken from Posin’s observations in nitrogen, 
The experimental curve, shown dotted, is that 
given by Whitehead” as the mean of the results 
given by a number of independent workers, [t 
will be observed that the calculated curve ig jn 
agreement with that experimentally determined 
for values of 5 from 10,000 down to about 100 
mm Xcm. For values of p5<100, the calculated 
values are higher than those observed, and the 
deviation increases steadily with decreasing p§ 
and thus with increasing X/p. The deviation 
occurs in the region of X/p=60 which is about 
the minimum value of X/p for which the second 
Townsend coefficient is observed.” Thus for 
values of p5 <100 the deviation can be accounted 
for by the presence of other sources of ionization 
in the gap. This is in agreement with calculations 
made in the past on a knowledge of the Townsend 
coefficients, and which have yielded approx. 
imately correct values of breakdown potential 
corresponding to low values of p6. The deviation 
between the two curves then indicates the region 
in which breakdown by streamer formation gives 
way to classical theory and the Townsend type 
of mechanism. 

Preliminary calculations of the breakdown 
strength of other gases, such as hydrogen, have 
been made and give values in agreement with 
those observed experimentally. 


CONCLUSIONS 


It is seen that the above theory postulates an 
entirely different mechanism of breakdown for 
longer gaps and higher pressures than that given 
by classical theory. 

1. The proposed mechanism depends on 
photo-ionization and electron multiplication in 
the gas (Townsend’s a). This photo-ionization is 
justified by recent observations on the effective- 


ness of such ionization under the existing condi- 


tions.'* "7, 
2. It requires that the density of positive ion 
space charge created in an electron avalanche is 


* F. H. Sanders, Phys. Rev. 44, 1020 (1933). 

* D. Q. Posin, Phys. Rev. 50, 650 (1936). 

26S. Whitehead, Dielectric Phenomena (D. Van Nostrand 
& Company, New York, 1927), p. 42. 

27 W. E. Bowls, Phys. Rev. 53, 293 (1936). 
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sufficient to draw in photoelectrons and their 

ny, produced by ionization by collision, in 

4 measure which renders the positive space 

charge a conducting plasma, and permits the 

further development of the streamer. This is only 

ible when the field strength and gas pressure 

are such as to allow a sufficient density of photo- 

jonization in the neighborhood of the streamer 
head to permit self-propagation. 

3. It shows that the development of long 
sparks is completely devoid of any dependence 
on cathode material, in conformity with ob- 
servation. 

4. It is in full accord with time-lag studies in 
that it involves formative time lags of ~ 10~’ sec. 
or less and explains the fact that when a poten- 
tial is applied to a gap in excess of the minimum 
required to cause breakdown midgap streamers 
develop and the formative time lag is reduced.® 
For in this case an electron avalanche which 
originates at the cathode will have developed the 
requisite density of positive ions to cause the 
initiation of a streamer before the anode is 
reached. The reduced distance which the ava- 
lanche has to travel before it produces a streamer 
makes it probable that several streamers will be 
initiated simultaneously. A conducting filament 
across the gap will be formed by the junction of 
such streamers, and results in a corresponding 
reduction in time of formation of the spark. 

5. At the sparking threshold the spark is 
initiated by a single electron liberated from the 
cathode. The chance of a spark then depends on 
(a) the chance of a single electron avalanche 
which produces at the anode or in midgap the 
requisite space charge, and (b) the success which 
this space charge has in photo-ionizing the am- 
bient gas and in attracting sufficient electrons to 
ensure self-propagation of a streamer. This 
accounts for the statistical time lags, and for the 
observed indefinite character of the sparking 
potential. That the sparking threshold is at all 
definite under these conditions is explained by 
the cumulative nature of electron ionization, 
Le., i=ipe**, and also that a/p=f(X/p) where 

J(X/p) is of the form Ae®*!? or C(X/p—D)?. 
-6. The theory is completely consistent with 
the observations of H. Raether on the develop- 
ment of sparks in a cloud chamber. Jn fact, 


Raether’s photographs clearly indicate the 
passage of an electron avalanche from cathode to 
anode followed by a streamer which proceeds 
from the anode to the cathode to form a con- 
ducting filament bridging the gap. Simultaneously 
and quite independently, Loeb, on the basis of 
corona studies, and Raether,” on the basis of 
cloud-chamber observations, postulated a quali- 
tative picture of the sparking mechanism which 
is here rendered quantitative. 

7. It can be shown that the field distortion 
produced by cumulative ionization with high 
densities of photoelectric current gives a lowering 
of sparking potential in rough quantitative agree- 
ment with the lowering observed by White," 
Rogowski and Wallraff. 

8. It clearly accounts for the filamentary 
character of spark breakdown at longer gaps and 
higher pressures and explains both the zigzag 
character of the spark and the branching ob- 
served. This differentiates the sparks under this 
mechanism from those in the Townsend regime 
(p6 <100 in air). 

9. It explains the fact that in the study of the 
Townsend coefficients at lower X/p and higher 
pé sparking occurs before values of the second 
coefficient, y, can be observed. 

10. The slight deviation of the theory from 
Paschen’s law is not inconsistent with experi- 
mental data in this field of investigation. 

11. While at present the inexact character of 
the quantitative criterion for streamer formation 
leaves considerable latitude owing to the ex- 
ponential nature of electron multiplication by 
ionization by collision, the equation gives 
numerical results in as good agreement with 
experiment as those based on Townsend’s theory, 
where a hypothetical value of y has to be as- 
sumed, since it cannot be observed. Numerical 
agreement of this character between theory and 
experiment must not be taken too seriously as a 
confirmation of any particular theory in view of 
the latitude given by the exponential term and 
the lack of adequate experimental data. 

12. The theory has been extended to the 
explanation of sparks of all lengths at higher pé 


28H. Raether, Zeits. f. Physik 112, 464 (1939). 
2° W. Rogowski 
183 (1937). 


ki and H. Wallraff, Zeits. f. Physik 102, 
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' in air and the necessary conditions for sparking 
have been established by Loeb. A more exhaustive 
discussion of the application of the theory will 
be given in an article to appear at some future 
date. 

13. It can be shown that the theory leads to 
an explanation of the breakdown of unsym- 
metrical gaps where the field distribution is 
known, and that quantitative calculations can be 
made. The theory may be applied to either surge- 
impulse or static breakdown. 

14. It is found that the self-propagating 
streamer mechanism, as modified for longer 
sparks, will adequately account for the pilot 


streamer in Schonland’s mechanism® for the 
lightning discharge, and it clarifies the 
proposed by the writer® for the “‘stepping” op. 
served in lightning and spark discharges, 

The writer wishes to express his appreciation 
to Professor L. B. Loeb, whose original qualita. 
tive theory stimulated this work and whose sug. 
gestions and criticisms have contributed much to 
the development of the present theory. He also 
thanks the Commonwealth Fund for the grant 
of a Fellowship, during the tenure of which this 
work was carried out. 


8° B. F. J. Schonland, Proc. Roy. Soc. A164, 132 
31 J. M. Meek, Phys. Rev. 55, 972 (1939). (189%) 
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The order-disorder transformation involving Ni;Fe affects both the ferromagnetic anisotropy 
and the saturation magnetization. Iron-nickel alloys in the range 65 to 80 percent nickel 


have been investigated, using spheroidal specimens. With superlattice formation the anisotropy 
becomes more like that of pure nickel with (111) as the direction of easiest magnetization, 
The change is largest near Ni;Fe. The saturation magnetization increases with ordering; the 
greatest observed increase was 5.8 percent for an alloy very near Ni;Fe. Different rates of 
cooling from above the critical temperature (about 490°C) affect the saturation magnetization. 
These changes are attributed to changes in degree of local ordering effects. Long distance order, 
induced by baking for long periods, greatly influences both anisotropy and magnetization at 


saturation. 


INTRODUCTION 


RON-NICKEL alloys in the composition 

range 65 to 80 percent nickel have unusual 
properties which are sensitive to heat treatment. 
They develop extremely high permeabilities 
when cooled rapidly from above 600°C. Baking 
below this temperature or cooling slowly does 
not produce uncommonly high permeabilities. 
Attempts to explain these characteristics in 
terms of an order-disorder transformation in- 
volving NisFe have long been hampered by the 


* Part of a dissertation eet to the Faculty of the 
Graduate School of Yale University in candidacy for the 
of Doctor of 
Now at the General Laboratories, United States Rub- 


ber Company, Passaic, New Jersey. 


lack of conclusive evidence for the existence of 
this superlattice.! 

Changes in resistance, which occur in these 
alloys with proper baking after quenching, are 
not as marked as in many order-disorder trans- 
formations. This is to be expected, however, 
since the differences in size and chemical proper- 


ties between iron and nickel atoms are not great. 


Specific heat measurements have not been con- 
clusive because the anomaly, due to the mag- 
netic transformation, occurs near the critical 
ordering temperature. X-ray analyses have, until 
after the experiments here reported were com- 

10. Dahl, Zeits. f. Metallkunde 28, 133-138 (1936); 


S. Kaya, J. Fac. Sci., Hokkaido Imp. Univ. 2, 29-53 (1938); 
F. E. Haworth, Phys. Rev. 54, 693-698 (1938). 


a 


pleted, given negative results. Small difference 
in x-ray scattering power between iron and 
nickel and the sluggishness of the transformation 
have probably been responsible for ambiguities. 
The ferromagnetic anisotropy of single crystals 
affords an excellent means of examining this 
order-disorder transformation, even though the 
anisotropy of alloys of interesting compositions 
issmall. Its high temperature coefficient indicates 
that it is sensitive to changes in lattice parameter, 
and lattice spacing usually decreases with order- 
ing. Furthermore, anisotropy may be uniquely 
affected, since a regular arrangement of two 
unlike atoms having different magnetic moments 
may well modify the energies required to mag- 
netize a single crystal in the principal directions. 
Calculations by McKeehan,? for purely magnetic 
interactions, predict such characteristic differ- 
ences between ordered and disordered states. 
The first reported measurements of changes 
with baking in the anisotropy and saturation 
magnetization of iron-nickel alloys near Ni;Fe* 
have since been extended to other alloys in the 
range 65 to 80 percent nickel by weight. These 
data would be difficult to explain by other than 
an order-disorder transformation, since they 
indicate changes throughout a single crystal 
specimen primarily related to the crystal struc- 
ture.* X-ray evidence presented since the experi- 
ments here reported®* has proved that this 
superstructure occurs. Exposures as short as 
eight hours have now revealed superlattice lines 
from alloys baked for 500 hours in the tempera- 
ture range 490°C to 370°C.5 These lines were 
absent from similar exposures with quenched 


specimens. 
SPECIMENS 


Large single crystals, of volume greater than 
1 cm’, were grown in a vacuum furnace by the 
slow freezing of melts in high purity magnesia 
crucibles. The stock materials used were homo- 
genized castings containing about 0.15 percent 
of the usual impurities: principally sulphur, 


*L. W. McKeehan, Phys. Rev. 52, 18-30 (1937). 

*L. W. McKeehan and E. M. Grabbe, Phys. Rev. 55, 
505 (1939); E. M. Grabbe and L. W. McKeehan, Phys. 
Rev. 55, 1142(A) (1939). 

especially simple to analyze. . W. G. Burgers and J. L. 
Snock, Physica 2, 1064-1074 (1935). 

*P. Leech and C. Sykes, Phil. Mag. 27, 742-753 (1939). 
*F. E. Haworth, Phys. Rev. 56, 289 (1939). 
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TaBLeE I. Description of specimens. 


Ni (weight percent) 80* 748 70.5 69.8 68.0(a)  68.0(6) 65° 
{west percent) 20° 25.2 29.5 30.2 32.0. 


F 32.0 35* 
Mass, 10% g 412 402 238 262 29.1 48.7 37.4 
Diameter, 10~ cm 5001 5106 4076 4495 4470 5364 5082 
10% cm 360 «3805 802-274 323 356 
Volume, 10-* cm* 476 470 280 309 343 578 44 
Axial ratio 168 135 164 15.0 16.6 143 
Observed 0.626 0.574 0.667 0.552 0.620 0.596 0.561 
) Calculated 0.651 0.548 0.669 0.560 0.606 0.552 0.635 


* Not analyzed. 


silicon, and cobalt, with traces of phosphorus, 
manganese, aluminum, and carbon. 

Specimens were prepared in the form of oblate 
spheroids having an equatorial plane of the form 
{110}. This plane contains all three important 
magnetic axes of forms (100), (111), and (110), 
so that the anisotropy of a given alloy may be 
completely determined from one specimen. A 
slab of correct orientation is sawed out, and a 
disk, cut from this slab, is cemented to the shaft 
of a high speed motor. A highly magnified sil- 
houette of the disk is projected on a vertical 
screen on which an ellipse of proper size and of 
the desired axial ratio has been carefully plotted. 
The disk, while rotating, is cut with fine files to 
match this contour.’ Exceptional care was taken 
in preparing the samples, and their x-ray Laue 
patterns show that distortion by cold work was 
trivial. Departures of the equatorial plane from 
the {110} plane aimed at were not allowed to 
exceed a degree or two. The effect introduced by 
such deviations of the axes chosen for measure- 
ment from symmetry axes in the crystal should 
be negligible. In each spheroid, the axial ratio 
(diameter to thickness) was made large enough 
(greater than 12) to ensure a low demagnetizing 
factor. The volumes were calculated from den- 
sities based upon x-ray data. Table I gives the 
compositions, in weight percent, of the spheroids 
used in this investigation, their dimensions, 
volumes, axial ratios, and demagnetizing factors, 
as calculated from dimensions and as experi- 
mentally determined. 


MEASUREMENTS 


Magnetic measurements were taken with an 
improved quinquefilar pendulum magnetometer. 
The weight was redistributed to increase the 
stability, and the sensitivity was greater than 
that of similar pendulums previously used.’ The 


7L. W. McKeehan, Rev. Sci. Inst. 5, 265-268 (1934); 
L. W. McKeehan, R. G. Piety and J. D. Kleis, Rev. Sci. 
Inst. 7, 494-497 (1936). 
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Fic. 1. J vs. H curves for 74.8 percent nickel. 


most important modification was in the method 
of measuring deflections. A small square of a 
diffraction grating replica fastened to a thin 
glass plate was mounted on the pendulum beam 
with its rulings vertical and in the plane of per- 
mitted motion. This was intensely illuminated 
from behind by a platinum ribbon lamp, using 
circulating water in a cell to absorb heat radi- 
ation. Highly magnified images of several 
adjacent rulings were thrown on a translucent 
glass scale by means of a compound microscope. 
At magnifications of 1200 to 1500, the line 
images were quite sharp. No fixed fiducial line 
is used; any line image may be chosen as the 
zero point for the subsequent deflection. By 
adjustment of the field gradient the deflection 
can always be made an integral number of line 
spacings, so that the displacement of the pen- 
dulum is simply this number multiplied by the 
grating constant. There is then no correction for 
lens curvature, and the magnification need not 
even be known. Furthermore, the deflection 
range is large, and the character of the images is 
such that slight transverse displacements of the 
pendulum, vertical or horizontal, do not diminish 
the precision of measurements. Displacements of 


the order of 0.5 mm were measured to about 1 
part in 1000 by this method. 

The magnetizing field was supplied by a 
solenoid, previously described,* with compen. 
sating coils to level the field near the center 
Two types of gradient compensating arrange. 
ments proved satisfactory; one with two series 
coils placed at the ends of the solenoid, the other 
utilized a single coil of 75 turns with reversed 


OERSTEDS + 411 


or 150 
om 150 300 
Fic. 2. ]—H curves for 68.0 percent nickel. 


current placed around the middle of the sole. 
noid.® The field was made uniform to 1 part in 
2000 for a distance of 8 mm along the solenoid 
axis. All the specimens saturated at low fields, 
and the scatter of points on the magnetization 
curves is small, indicating that the precision 
expected was attained. 


HEAT TREATMENT 


After shaping, the spheroids were annealed in 
hydrogen for 1 hour at 750°C to relieve strains. 
The cooling rate, attained by permitting the 
specimen to cool with the furnace after annealing, 
was about 300 degrees per hour. This cooling 
rate was slow enough to develop local order to 
a small degree. (See Results II.) 


8 L. W. McKeehan, R. G. Piety and J. D. Kleis, Rev. Sa. 


Inst. 7, 494-497 (1936). 
*L. W. McKeehan, Rev. Sci. Inst. 10, 371-373 (1939). 


In order to retain the disordered state at room 
temperature, a method for extremely rapid 
cooling was devised. The specimens were placed 
in the bend of a fused silica U-tube through which 
hydrogen at about atmospheric pressure was 

. This bend of the tube was lowered into 
a furnace with the ends projecting. When the 
temperature had reached 750°C (well above the 
critical temperature), the tube was quickly 
withdrawn and inverted. The specimens fell 
from the red-hot bend of the tube into a cool 
part of the tube. Due to its small volume and 
large surface, a spheroid ceased to glow while 
falling. The cooling rate was estimated to be 
200 to 400°C per second, Maximum disorder was 
obtained by this drastic treatment. . 

The most favorable ordering temperature, in 
agreement with Kaya,’ was found to be near 
490°C. Specimens were baked at this temperature 
for at least a week in evacuated Pyrex tubes. 
Magnetic measurements were then made, and 
the spheroids were returned to the furnace and 
baked for another long period at a slightly lower 
temperature. If no changes occurred during the 
second baking period, it was assumed that 
equilibrium had been attained. It is probable 
that equilibrium was reached in all cases. In 
three instances there is additional ground for 
this opinion, since the same final condition was 
reached in consecutive programs of baking 
separated by a quenching process. 


RESULTS 


I. Magnetic anisotropy 


Baking for long periods at or below 490°C 
caused large changes in the ferromagnetic ani- 
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sotropy of the crystals. For an alloy of 74.8 
percent nickel and the remainder iron (very 
nearly the correct stoichiometric ratio for Ni;Fe), 
the disordered state exhibited nearly zero ani- 
sotropy. In the ordered state the anisotropy ob- 
served, measured by the first anisotropy coef- 
ficient (K,), was almost two-thirds that reported 
for pure nickel. The direction of easiest mag- 
netization was [111], with [110] and [100] fol- 
lowing in that order. Fig. 1 shows the magnet- 
ization curves for this spheroid. 

For 80 percent nickel, the anisotropy after 


TABLE II, Magnetization energy differences. 


ENERGIES IN 10 ERG-CM~* FOR FOLLOWING WEIGHT PERCENTS OF NI 
80 74.8 70.5 69.8 68.0 65 
Wiio— Wioo Ordered —0.31 —0.58 0.08 —0.03 —0.05 0.21 
Wie— Wioo Disordered —0.01 0.00 0.40 0.27 0.26 0.29 
Wioo Ordered —0.48 — 0.83 —0.25 —0.34 —0.19 0.17 
Win— Disordered — 0.03 0.00 0.25 0.21 0.22 0.26 
K, Ordered — 1.26 — 2.30 0.31 —0.11 —0.18 0.83 
K, Disordered —0.04 0.00 1.60 1.07 1.04 1.16 
K: Ordered —1.75 — 1.65 —9.52 — 8.24 —3.55 — 2.98 
K; Disordered —0.33 0.00 —7.76 — 3.82 —3.40 —3.41 
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_ Fig, 4. (A) Saturation magnetization vs. composition for 
iron-nickel alloys. (B) Percent increase in Jsar. with order- 
ing vs. composition. 


quenching is very small, with [111] as the direc- 
tion of easiest magnetization. Baking again 
develops a large anisotropy similar to that for 
the 74.8 percent alloy, but the change in K; is 
less. 

For crystals containing less nickel than Ni3Fe, 
the irregular order” of decreasing ease of mag- 
netization [100], [111], [110] was observed at 
low fields for the disordered state. At higher 
fields, however, the [110] curve jumps to satu- 
ration, crossing the [111] curve. This crossover, 
which has not before been reported, was observed 
in the 65, 68, and 70 percent nickel alloys. In 
going over to the ordered state, a complete 
changeover occurs for the 68 and 70 percent 
nickel alloys. The [111] direction becomes the 
direction of easiest magnetization at all values 
of H. The [110] and [100] curves cross and 
recross, with [100] reaching saturation last. 
This crisscrossing is for almost coincident curves 
and may be due to errors in cutting the spheroids. 
That it is an inherent effect is, however, rendered 
more probable by the fact that it was observed to 
occur in three different specimens. Fig. 2 shows 
the upper part of the magnetization curves for 


10 J. D. Kleis, Phys. Rev. 50, 1178-1181 (1936). 


a 68 percent nickel, 32 percent iron alloy for the 
ordered state, the disordered state, and an inter. 
mediate condition. The changes in aniso 

are small in comparison with those for the 74. 
percent alloy. 

For 65 percent nickel, this complete change. 
over of magnetization curves does not occur, 
(See Fig. 3.) The irregular order [100], [111], 
[110] is observed for the disordered state, and 
the anisotropy merely decreases when the super. 
structure is formed. It follows that for some 
composition between 65 and 68 percent nickel, 
the anisotropy must be vanishingly small for the 
ordered state. 

The changes in anisotropy between various 
degrees of disorder, induced in a specimen by the 
different cooling rates given above, were negli- 
gible. The large changes took place only after 
prolonged baking. Table II tabulates the anisot- 
ropy constants and the energy differences 
Wioo— Wiro and Wii—Wioo for ordered and 
disordered crystals. The general effect of super- 
structure formation is to make the anisotropy 
more like that of pure nickel. 


II. Saturation magnetization 


The large changes observed in the saturation 
magnetization were unexpected. Considering the 
amount of work that has been done on the alloys 
in this composition range, it is surprising that 
such large effects have not been detected before." 
The difference in saturation magnetization 
between ordered and disordered states was a 
maximum for the 74.8 percent nickel alloy. For 


TABLE III. Saturation magnetization values. 


SATURATION VALUES 
wt 
Comme Cole PERCENT INCREASE 
PERCENT (a) LocaL (Dis- 

NICKEL | ORDERED RDER) ORDERED) TO (a) T0 
80.0 908 895 882 1.4 3.0 
74.8 1007 963 952 4.6 5.8 
70.5 1029 986 44 
69.8 1071 1044 1028 2.6 4.2 
68.0 1102 1080 1066 2.0 34 
68.0 1108 1082 1072 2.4 3.4 
65.0 1129 1117 1106 1.2 2.2 


* 10°C to 100°C per minute (normally expected values). 
** About 200°C per second. 


1 Confirmatory results have since been reported by P. 
Leech and C. Sykes, Phil. Mag. [7] 27, 742-753 (1939) 
and by S. T. Pan, Phys. Rev. 56, 933-936 (1939). 
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disorder attained by rapid or normal cooling, the 
change amounted to 4.6 percent, while for drastic 
cooling, the difference between ordered and dis- 
ordered states was 5.8 percent of the saturation 
value for the disordered crystal. These values 
could be reproduced very accurately, regardless 
of previous heat treatment, and, therefore, the 
difference of 1.2 percent for different cooling 
rates is regarded as due to changes in the degree 
of short distance order. 

For 65, 68, and 70 percent nickel, the changes 
with ordering, though marked, were smaller and 
decreased with increasing departure from Ni;Fe. 
One of the specimens of composition near 70 
percent nickel was difficult to magnetize to 
saturation and exhibited an abnormally low 
saturation value. The change with ordering, 
however, fits in with the other values. (See Fig. 
4.) The increase in saturation magnetization with 
ordering falls off more rapidly on the nickel-rich 
side. At 80 percent nickel the maximum change 
observed was 3.0 percent, while near 70 percent 
nickel the shift amounted to 4.3 percent. 

In line with this, it may be noted that ab- 
normal magnetic properties extend over a wider 
range below 75 percent nickel than above this 
composition. This supports the view that the 
order-disotder transformation is responsible for 
permeability changes as well as for saturation 
changes. 

The difference in the saturation magnetiza- 
tions for different rates of cooling was of the 
same order of magnitude for all the alloys tested. 
Table III lists the saturation values obtained 
for different heat treatments. 

The dissymmetry of the ordering forces, which 
is suggested by these facts, implies that the 
ordering forces which act below the critical tem- 
perature for order are greater if, at this tem- 
perature, the degree of ferromagnetic saturation 
is higher. The Curie point has its maximum near 
70 percent nickel, and the magnetization at any 
selected temperature below the maximum Curie 
point has its maximum value at a still lower con- 
centration of nickel, because the magnetic mo- 
ment per atom rises with addition of iron. It is 
therefore suggested that ferromagnetic interac- 
tion promotes order in this series of alloys. 

Curves plotted for the two (111) axes lying in 
the equatorial plane of the 74.8 percent nickel 
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alloy proved to be identical. This test, suggested 
by Dr. R. M. Bozorth of the Bell Telephone 
Laboratories, shows that the effects are not 
determined by the accidental direction of a stray 
magnetic field present during baking, but are 
immediately related to the structure. 


DISCUSSION 


The superstructure transformation involving 
NisFe is a sluggish one, with a critical ordering 


© ORDERED STATE 
(A) © DISORDERED STATE 
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/ 
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Fic. 5. (A) Ki vs. composition for iron-nickel alloys. (B) 
Change in K, with ordering vs. composition. 


temperature near 490°C. Baking just below this 
temperature was most effective in ordering. 
Both the anisotropy and saturation magnetiza- 
tion are changed by baking for long periods 
below the critical temperature. Cooling at dif- 
ferent rates produces appreciable changes only 
in the saturation magnetization. Hence, the 
anisotropy gives a measure of the long distance 
order, while the saturation magnetization change 
is related to local ordering, or the average mag- 
netic moment for any one atom is determined by 
the nature of its nearest neighbors. Fig. 4 gives 
plots of magnetization at saturation versus com- 
position and change in magnetization with 
ordering versus composition. The maximum is 
near the critical composition. Fig. 5 shows the 
change in the anisotropy constant K, by super- 
structure formation for different compositions. 
The maximum at Ni;Fe is quite broad. The values 
for K, (disordered) are slightly lower than values 
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obtained by Kleis. The second anisotropy con- 
stant K, becomes more positive with ordering 
except for the 65 percent nickel alloy. 

The high permeability developed in these 
alloys by quenching is characteristic of the 
highly disordered state. A plausible explanation 
for the normal low permeability is that slow 
cooling forms small ordered regions which are 
frozen into an otherwise disordered lattice and 
give rise to strains. These strains increase the 
coercive force, and the permeability is low, while 
the highly disordered state is relatively strain- 
free and exhibits an extremely high permeability. 

Cooling alloys in this composition range from 
above 500°C in the presence of a magnetic field 
also produces high permeabilities.” The effects 
of such heat treatment in a magnetic field may 
also be related to ordering. The critical tem- 
perature may be altered by the application of a 
field during cooling, if the ordering forces are in 
part of ferromagnetic origin. It does not appear, 
however, that the saturation value can be much 
affected by such treatment, as the following 
additional test indicates. Two specimens of 74.8 
percent nickel and 68 percent nickel were each 
cooled slowly six times from 750°C with mag- 
netic fields of 25 and 800 oersteds directed in 
turn along the three principal magnetic axes. The 
saturation magnetizations were the same as for 
cooling in the absence of a field, indicating, in 
accordance with the views expressed above, that 
these fields had no effect on the local order. 

The addition of small amounts of the elements 
Mo, Mn, and Cu to these alloys increases the 
permeability to high values without the necessity 
of quenching or cooling in a magnetic field.’* If 
the existence of local order is related to the per- 
meability in these cases, the ordering forces are 
likely of a ferromagnetic nature, since ferro- 
magnetic properties of pure alloys are easily 
disturbed by impurities in small quantities. 

If the superstructure FeCo exists, indications 
are that it has an ordering temperature near 
800°C, and it would be almost impossible to 
obtain the disordered state at room temperature." 


12 R. M. Bozorth, J. App. Phys. 8, 575-588 (1937). 

%S. Kaya and M. Nakayama, Zeits. f. Physik 112, 
420-429 (1939). 

4 A, Kussmann, B. Sharnow and A. Schultze, Zeits. f. 
tech. Physik 13, 449-460 (1932). 


In a preliminary test for such order, if it ex; 
one iron-cobalt specimen containing 56.3 percent 
iron (near FeCo) was repeatedly drastically 
quenched from temperatures ranging from 800°C 
and 900°C. The anisotropy and saturation valye 
were not affected. The iron-cobalt spheroid had 
an axial ratio of about 20 to 1 and, Consequently 
saturated easily. Its magnetization at saturation 
was 1937. 

Kaya’ found that polycrystalline samples of 


compositions near Ni3Fe were more difficult to 


saturate after they had been baked for long 
intervals. Much earlier than this, Elmen® had 
shown that such treatment greatly lowered 
initial permeability. This is in agreement with 
the results here presented, which show the ap- 
pearance in the ordered state of one direction 
of very difficult magnetization. 

As mentioned in the introduction, the reported 
changes in resistance with ordering are small, and 
so too are the possible changes in lattice param- 
eter. It may be worth noting that the change in 
resistance is consistent in sign with increase in the 
number of free electrons and corresponding 
increase in the vacancies in the 3d electron group. 
The latter change corresponds to the observed 
increase in magnetic moment per atom. 

McKeehan’ has calculated the anisotropy 
constants K, and Kz for Ni;Fe as a disordered 
solid solution and as a superstructure, from mag- 
netic multipole interaction. The quadrupole 
terms predict the correct sign for the change in 
K, accompanying the order-disorder transforma- 
tion and are of the correct order of magnitude. 


For Ni;Fe in the disordered state, K, (calculated): 


= —0.1810* erg-cm-*, while for the ordered 
array, K, (calculated) = —1.62X10* erg-cm~. 
These are to be compared with the experi- 
mentally observed values, zero and —2.30X10' 
erg:cm~, respectively. The computed sextupole 
part of K,, and the computed Ke, which has no 


quadrupole part, do not change in the observed 


directions for superstructure formation. 

The author wishes to express his sincere thanks 
to Professor L. W. McKeehan, who suggested 
this investigation, for his ready advice and for 
stimulating discussions during the course of the 
work. 


%G, W. Elmen, J. Frank. Inst. 207, 583-617 (1929), 
esp. Fig. 13. 
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It is shown for Poisson distributions that both the probable error and the probability that a 
deviation not exceed the standard deviation approach as limits the corresponding values for 
normal distributions. Curves are given showing the manner of approach. 


F occurrences of a certain event are subject to 
I a normal distribution, probable error is 
defined as that deviation from the average for 
which there is a half-chance that a single trial 
should not yield a greater deviation. For a 
normal distribution 


probable error = 0.674490 (1) 


where o is the standard deviation. For an un- 
symmetrical distribution, this expression does 
not give the probable error; or if it is taken as 
the defining expression of probable error, it does 
not have the same meaning as for normal 
distributions. 

One type of unsymmetrical distribution is that 
which follows a Poisson, or Bateman, distribu- 
tion formula. This formula has received con- 
siderable attention recently because of its appli- 
cation to counts from Geiger-Miiller counters. 
It may be written as 


W,=x"e-*/n!, (2) 


and gives the probability of occurrence of n 
events in a single period of observation for which 
the true average is x events. The purpose of the 
present paper is to discuss probable error for 
distributions following the Poisson formula. 

We shall first take up the question as to the 
probability that a single observation should yield 
a result not differing from the average by more 
than the standard deviation. For a Poisson dis- 
tribution, the standard deviation is simply /x. 
The desired probability is then 


n=z+vz x" 
P.= dS (3) 


n! 


This expression is easily evaluated when x is a 
perfect square and not too large. When x is not 
a perfect square, it can be evaluated only ap- 


proximately. It is found that its value slowly de- 

creases as larger and larger values of x are used. 

This suggests seeking its limit as x approaches 

infinity. To do this it simplifies things somewhat 

if we set x=y’, and set n=y*+t. We then have 
(y?)u*+# 


= Ly 


It can be shown that the difference between 
any two consecutive terms in this sum is an 
infinitesimal of higher order than either of the 
terms as y becomes large. If for any chosen value 
of x in Eq. (2), the factorial be replaced by a 
gamma-function and m be considered to vary 
continuously, a smooth curve may be drawn with 
W as ordinate and m as abscissa. When this is 
done it is seen that the portion of the area under 
the curve bounded by the ordinates correspond- 
ing to any two consecutive integral values of n 
is numerically greater than one of these ordinates 
and less than the other. These two facts enable 
us to replace the sum in (4) by an integral with 
respect to ¢ between the limits —y and +y. 
Making the substitution t=ky, we get 


+1 (y?)u*+ky 
P,=L ———e"ydk. 5 
yon (y?+ky) ! yd ( ) 


If the factorial be considered as a gamma- 
function, the integrand is a bounded continuous 
function of y for y>0 and any chosen value of k 
in the range of k. Also its simultaneous limit 
exists as k-ky and y>~, where kp is any finite 
value of k. To show this it is only necessary to 
set k= f(y), where f(y) is any continuous function 
of y that has the limit ko as y>~, and than take 
the limit of the result as y>. If this is done, it 
is found that the limit is independent of the form 
of the function f(y). Thus the integrand con- 
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Fic. 1. The upper curve 


verges uniformly in the interval —1<k< +1, 
and we may interchange the limit and integral 
sign.! This gives 
+1 
P.= L ———e~ydk. (6) 
(y?-+ky)! 


To evaluate the limit, we replace the factorial 
by the Stirling approximation and proceed as 
follows: 

2)y*+kye—v? 
L F(k,y)= L 


1 
L log Fi(k, AL k 


+(y?+ky) log —— 


ytk 
+ 
( ytk 2(y-+k)? 


1 See E. yy Functions of Real Variables (Henry 
., 1928), pp. 122, 350, 360. 


shows the manner in which 


— — — the probability that a de. 


viation not to exceed ¢ 
as x 
the 

in which the ratio of prob 
able error to « approaches 


700 800 900 1000 


1 
2 
or 1 
L Fk, y) =——e*™”. 
(k, y) (an)! 
Therefore 


2 1 
=—— f =0.68269, (1) 
(2x )! 0 


which is just the probability that a deviation not 


exceed o when the distribution is normal. This 
was to be expected because it is known that the 
Poisson distribution becomes more and more 
similar to the normal distribution as x becomes 
large. 

We can now evaluate the limiting probability 
that a result not deviate from the average by 
more than mo, where m is arbitrary. The same 
line of mathematical reasoning as above leads to 


the result a 
P.=(-)  @ 


which is just expression (7) with different limits 
of integration. If the integral in (8) be set equal 


0.67449 as x becomes large, ° 
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to 0.50, it is known that m must equal 0.67449. 
This means that in the limit, expression (1) gives 
the probable error when the distribution follows 
a Poisson formula. 

In Fig. 1 the upper curve shows graphically 
how the probability that a deviation not exceed 
« approaches 0.68269, and the lower curve shows 
how P.E./o approaches 0.67449 as x becomes 
large. These curves were obtained by direct 
computation for small values of x, and by ap- 
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proximation methods for values of x greater 
than 1000. In both cases it is seen that the ap- 
proach is rather rapid for small values of x. Also 
it is to be noted that, for finite values of x, P, 
is always greater than its limit while P.E./c¢ is 
always smaller than its limit. This means that the 
error committed in saying that P.E. =0.67449c, 
or that committed in saying the probability that 
a single deviation not exceed ¢ is 0.68269 are both 
theoretically on the conservative side. 
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Osmotic Diffusion in Gases 


E. J. HELLUND 


VOLUME 57 


University of Washington, Seattle, Washington 


(Received January 13, 1940) 


The theory of diffusion induced in a ‘‘solvent”’ gas, by the pressure diffusion of two other 


gases, is developed on the basis of formulae derived in a previous paper by the author. The 
analysis is applied to the Loschmidt experiment and the pressures generated, are found 
analytically as a function of the time elapsed, from the onset of diffusion. All formulae are given 


INTRODUCTION 


——- and gases, because of the mobility 
of the molecules of each exhibit a series of 


closely allied phenomena. Particularly is this 
shown in the phenomenon of diffusion, and es- 
pecially with respect to that class of diffusion 
associated with pressure gradients or external 
force fields. The latter seems intuitively to be 
the most understandable and has as a conse- 
quence yielded to the most elementary treat- 
ments, results which are astonishingly accurate. 
In other cases, however, notably in the case of 
thermal diffusion in gases and osmotic diffusion 
in liquids, elementary theories have been power- 
less to explain the experimental results. In the 
case of thermal diffusion in gases it was found 
necessary to develop the theory from a continuity 
equation, the foundations of which seem fairly 
secure. The kinetic theory of liquids, unfortu- 
nately, lacks any such rigorous attack and has, 
as a consequence, given no satisfactory interpre- 
tation of osmotic diffusion. It has therefore 
seemed of interest to examine the corresponding 


in terms of the binary diffusion coefficients of the component gases. 


diffusion in gases even though qualitatively the 
phenomena in gases and liquids do possess essen- 
tial differences which vitiates somewhat the at- 
tempt to establish a close correlation between 
the two. 

Osmotic diffusion in liquids is said to occur 
when motion of a solvent through a _ semi- 
permeable membrane takes place due to the 
action of a solute which cannot pass through the 
membrane. That is, three essential elements are 
present for the process to occur; a solvent, per- 
vading the entire vessel; a substance dissolved 
in a portion of the solvent ; and a substance which 
may be said to occlude the first solute from that 
portion of the solvent which the second solute 
occupies. The membrane here has been desig- 
nated as a solute occupying a certain portion of 
the solvent, and having no effect on the solvent 
in which it exists. Physically this difference in 
point of view is immaterial, but formally it 
expresses more aptly the conditions to be ex- 
pected for such diffusion to occur in gases. 

One may then define osmotic diffusion in gases 
as that diffusion which would occur in a tube 
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pervaded by a gas of constant number density, 
(the solvent), due to the existence in the tube of 
two or more gases of varying concentration 
(solute), which would themselves diffuse because 
of their respective pressure gradients. More 
specifically one may expect in an experiment of 
the Loschmidt type, where a long tube is sepa- 
rated in the middle by a diaphragm, that if, 
initially, one gas pervades both sections with the 
_ same pressure, and if in addition two different 
gases are placed in the two sections, one on each 
side of the diaphragm, that diffusion of the first 
gas will occur when the diaphragm is removed. 
That is, a surge of the first gas will take place in 
a direction of the tube depending on the charac- 
teristics of the other two gases and their position 
in the vessel. The osmotic diffusion current may 
then be expressed as proportional to the number 
density gradient of either of the two gases of 
varying composition, and the constant of pro- 


portionality may be taken as the osmotic diffy. 
sion coefficient. Theoretically this coefficient Was 
shown in a previous paper by the author,! to be 
expressible in terms of the binary diffusion oo. 
efficients, and the densities of all the gases in the 
mixture. 

The Loschmidt experiment does not determine 
directly the osmotic diffusion coefficient since the 
results are necessarily complicated by Pressure 
diffusion. Nevertheless, the observation of the 
maximum pressures generated, provides a crogg 
check on the binary diffusion coefficients, as 
these pressures can be calculated in terms of these 
coefficients. 

In the theory to follow, the results, as in the 
theory of the determination of the binary diffy. 
sion coefficients by the Loschmidt experiment, 
will be carried to that order of approximation 
which is deemed sufficient for the comparison 
with the experimental results. 


THE EQUATIONS OF CONTINUITY 


In the following, for purposes of simplicity, only mixtures of gases containing three components 
will be considered. One finds for the continuity equation of the number densities, the result? 


ot pnd \D,, me(ne-+n;) Da; 
where 1+j#k. 

Dy2D 13D 32 a 

ond ap 
Zi= — pjXj5 (3) 


In the above equations, the X; are the components of external force per unit mass acting in the 
x direction (only variations in the x direction are considered). The diffusion coefficients, also, are 
defined in terms of the particle currents, and are not to be confused with the mass coefficients em- 
ployed in the paper previously cited. The coefficients appearing in (1) and (2) are thus, those which 
express the current moving past a stationary cross section of a constant pressure gas tube. 


Let us now consider a tube of length ‘“‘a” cm. A diaphragm divides the tube into two equal parts _ 


at a distance ‘“‘a”/2 cm from each end. Initially gas ‘‘3”’ occupies both sections at the same density, 
and gas ‘‘1” occupies the left section and gas ‘‘2” the right section. Gas ‘‘3”’ thus is the solvent, and 
the object of the investigation is to determine the behavior of this component with time, after the 
diaphragm is removed. In the experiment a constant temperature will be assumed, and also the 
only external forces considered will be those at the walls of the tube. 


1E. J. Hellund, Phys. Rev. 57, 319 (1940). 
2 Reference 1, Eqs. (46) and (47). 
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Since the pressure is a constant throughout the tube, one has 


On, One 
(4a) —+—+— =0, (4b) 
Ox Ox Ox 


which, inserted in (1), leads to the two simultaneous partial differential equations which describe 
the concentrations at any point of the tube, with the time. 


1 1 _ 9 fans 1 + | 6) 
dt Ax’ (m+m)Di dxLdx (me+m3)Ds2) nA J 


In (5), the coefficient of dn2/dx is thus the coefficient of osmotic diffusion and will be written, 


{ 
pa;3:2)=(“). (7). 
nF | 


Because of the character of the divisor, this coefficient is not necessarily constant for a constant 
total pressure and constant solvent density. The extent of the possible variation under the above ° 
conditions depends on the magnitudes of the various binary diffusion coefficients. Because of the 
positive character of all quantities appearing, and near equality of the D;; the variation, however, is 


not large. 
For times which are small one has approximately 


—dn;/dt=D(1; 3; 2)d°n2/dx?, (8) 
—On2/dt= —D(2; 1, 3)d2n2/dx?, (9) 
where 
1 Ne 1 
DQit,3)=—|{ + (10) 
Ds, nA 


In (9) the variability of the coefficient of dn2/dx has been ignored, and in (8) and (9) nae and 
d°n;/dx* have been set equal to zero. 
From (8) and (9) one obtains the result for m; that, 


n3(t) ex x<a/2; 
D(2;1, 3) 
where’ 
—4— cos (=) sin (=) exp | 1, | (12) 
wv 1 a 2 a 


The superscript, 0, is used to denote the initial value of the quantity for :=0. Eq. (11) tells little 
that was not immediately apparent from elementary physical considerations. The changes in density 
of ns take place most quickly where m2 changes most rapidly and remain constant for constant m2 
values. What is of greatest interest to the experimenter is not the variation of density throughout the 
tube, but the change in the total number of particles of a given kind in one section of the tube. 

Because of the increasing seriousness of the approximations made in (8) and (9) one cannot use (11) 


*Loeb, Kinetic Theory of Gases (McGraw-Hill, 1927), first edition, p. 233. 
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to obtain this information. One may, however, integrate Eqs. (5) and (6) directly as follows, 


ON; One Ons 
ot ax Ox 
ON: Ons One 

—— =—D(1; 2; 3) -—D(2; 1, 3), (14) 
ot ax Ox 


where large N denotes the total number of particles of a given kind in the tube between the limits | 


x=0, x=a/2. On the right side of (12) all quantities are evaluated at the point x=a/2, the values at 
x=0, not appearing because of the boundary conditions that 


On;/dx 0. (15) 


SOLUTION OF THE CONTINUITY EQUATIONS 


The relations expressed by (13) and (14) are exact and may be conveniently used to determine the 
functions Nand Nz. It is necessary however to transform, first, the partials of the densities, evaluated 
at x=a/2. 

If one assumes a Fourier expansion 


m= CH) c0s (=). (16) 


satisfying the boundary conditions dn;/dx=0 for x=0, or a, then one obtains, carrying only the 
first three terms of the series, 


ax = a i (17) 
One has then, if variables N;— Nj ]in«. = N;* are introduced, 
ON;* 
=—[D(1; 3; 2)N2*—D(3; 1, 2)N3*], (18a) 
ot a? 
=—[D(1; 2; 3)Ns*—D(2; 1, 3)Na*]. (18b) 
a 
One finds readily that the solutions, fitting the initial conditions are, 
exp| f xat| sinh (— f (19) 
t=0 2a? 0 2a? 0 
pt 
N2*= -»,| exp sinh (20) 
t=co 2a? 0 2a 0 
® t=0 


where 


nN, 1 1 ne 1 1 
(mit+ms)D1s (n)A\ (my +n2)Dis 


+ ) (22a) 


| 
Ww 
in 
(2 
in 
pe 
m 
ot 
te 
| pe 
p= in 
co 
sO 


(13) 


(14) 


limits 
ues at 


(15) 


le the 
lated 


(16) 
y the 


(17) 


18a) 


18b) 


(19) 


(20) 


(21) 


12a) 


OSMOTIC DIFFUSION IN GASES 741 


x={—+—+—} 22b 

. 
Ne { 1 1 

v= (22c) 
(n)A| | 


Remembering that these functions are all evaluated at x=a/2 and hence are not strongly depend- 
ent on the time, one effected the integrations by multiplying the integrands (evaluated for t= ~), 
by the elapsed time. For the binary diffusion case such an assumption is in fact exact, but in our 
case variations in time do occur. However even these variations are not appreciably effective, because 
of the character of the integrands, and the condition that n=n,-+-n2+n3=constant. 


DISCUSSION OF THE SOLUTIONS 


In view of the observations following (22), we shall use the following approximations for N(¢) 
and N 2(t) ’ 


N,(t)=Ns] inh -(@*—409)! 2 

5=tanh~! (25) 


where in the above equations %, VY, X, and Q are all evaluated at t= «. It may be easily shown that 
in the case of vanishing interaction between gases ‘‘3” and ‘‘1” and “3” and “2” that the solution 


(24) becomes 


N2= Ne |tmx(1—exp [ — }). (26) 


This compares with the exact series, derivable from (12), 


8 1 Or? 
1 -—(exp exp | ) (27) 
9 a? 


In practise only three terms of the summation 
in (27) are employed, so that initially this ex- 
pansion is unacceptable whereas (26) has been 
made to fit this boundary condition. On the 
other hand, for long elapsed times, the second 
term of (26) is in error by approximately 20 
percent when compared to the summation in 
(27). The source of this discrepancy is to be found 
in the approximation (17) for very short times 
after the diaphragm has been removed. One 
could improve the agreement by using (11) and 
(12) for short intervals and then joining to the 
solutions of (18a) and (18b). The gain in accuracy 


would not seem to compensate for the increased 
complication in the resulting solutions. 
Inasmuch as the initial high density gradients 
make doubtful the applicability of the diffusion 
equation for small elapsed times, it would seem 
more logical to determine Nz and Ns; experi- 
mentally at some later time, say, tr. Then the 
solutions for subsequent times could be deter- 
mined by evaluating the arbitrary constants for 
t=r. For moderate values of + error by this 
method becomes negligible. In connection with 
the above, it is perhaps of importance to point 
out that physically the Loschmidt experiment 
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Fic. 1. Plot of the CO pressure variation, in the section 
from x=0 to a/2, as a function of the elapsed time. 


must always introduce some uncertainty in re- 
gard to the initial conditions by virtue of the 
nature of the diaphragm itself and of its manner 
of removal. 

If then one solves (18a) and (18b) in this 
manner the resulting expression for N3 (and for 
N2) is materially altered. The equations which 
correspond to the group (23), (24), and (25), 
together with the results derived from them may 
be obtained quite easily and hence will not be 
given here. 

Returning to Eq. (23), one may evaluate the 
maximum pressures generated by gas “3."’ The 
result is, 


20P 2 (X—@’)/26’ 
\x+e’ 


where ® = (¢?—40y)!. The time at which this 
occurs is 


Ps, max =P 3 


t (29) 


The pressures P; are those for a uniform dis- 
tribution of m3 from x=0 to a/2. It may be no- 
ticed that Q is itself proportional to P;. By the 
conditions assumed initially, mz and m3 are equal, 
and thus the apparent lack of symmetry in (28) 
is spurious. The second term of (28) will vanish 
likewise if the cross section of gases “1” and ‘‘2” 
and their masses are equal. The osmotic pressures 
thus depend, as does thermal diffusion, on the 
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differences in the properties of two gases, Such 
phenomena are of infrequent occurrence and of 
great importance in the determination of force 
interactions between molecules. 


NUMERICAL RESULTS 


The variations of the osmotic pressure wit, 
time were calculated for a mixture of Hg, CH, 
and CO: which were numbered in the orde 
given to correspond to the previous th 
Thus CO: was taken to be the solvent os 
in the above. The number densities were gy 
equal to one another and of such a magnitug 
that the total pressure was 114 cm Hg. Th 
temperature used was 0° centigrade. For th 
above conditions, then, experimental values of 
the various binary diffusion coefficients are. 
D2=0.625; Dis=0.538; and D23=0.153. Fo 
these values one finds N3*/Ne2)max=0.179 which 
is equivalent to a pressure of 6.82 cm Hg 
Also one finds in this case that: = —0.0135, 
X=+0.521, V=D(1;2;3)=—0.129, 
3; 1)=+0.122, (#—42¥)'=0.251, and 
a? |max = 4.18. 


SUMMARY 


The continuity equation of a mixture of thre 
gases has been solved by approximate methods 
for the case where the distribution of one of the 
gases of the mixture is initially uniform. It is 
pointed out in connection with the solutions, 
with respect to the space-wise distribution, that 
the diffusion coefficient is not a constant and 
that the neglect of the composition dependence 
over considerable time intervals is not justified. 
First order solutions are obtained for the pressure 
surges of the “‘solvent’’ gas in the Loschmidt 
apparatus and the results (Fig. 1) plotted fora 
CO:2z, He, CH, mixture. In connection with the 
maximum pressures generated in the “solvent” 
gas of a ternary mixture it is possible to show 
from (28) that if gases ‘‘3,’’ ‘‘2,” “‘1”’ (by cyclic 
permutation of all three) are successively the 
solvent gases for the same partial pressures (ie, 
P,=P:=P;=}P), that the successive mat 
imum pressure differences occur in the ratios, 
D(2; 3; 1) : D(3;1; 2) : D(1; 2; 3) and that ther 
sum vanishes identically. 


De 
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Isotope Separation by Transient Pressure Diffusion 


E. J. Hettunp 
University of Washington, Seattle, Washington 
(Received February 9, 1940) 


The separation of molecules of nearly equal masses is investigated in connection with the 
Loschmidt experiment. The results are obtained by a perturbation treatment of the continuity 


equations for a ternary mixture. 


ARIOUS means have been employed in the 
separation of particles of nearly equal 
masses and effectively identical cross sections. 
Many of these depend on diffusion processes, 
which operate by virtue of the dependence of the 
velocity of migration on the mass of the particle 
in question. Most effective of these have been 
methods utilizing the diffusion through a porous 
solid. The advantage in the latter has been that 
the operation could be made continuous (e.g., 
the Hertz apparatus). Frequently, however, these 
methods operate discontinuously, requiring repe- 
tition after the attainment of an equilibrium 
state in which the separation is a maximum. 

Pressure diffusion, of gases alone, has been 
used in the ultracentrifuge wherein the high 
fields produced are balanced by opposite pressure 
gradients. In the equilibrium state no diffusion 
current exists and the separation depends in no 
way on the cross sections of the molecules of the 
mixture, no matter how many components 
there are. 

In the method to be described in the present 
paper the separation is dependent on the exist- 
ence of a third gas and one finds a dependence on 
the cross sections as well as the masses. In this 
respect it differs from most other methods of 
separation. The theory is based on the solution 
of the continuity equations of a ternary mixture, 
derived in a previous paper by the author.’ 

If one denotes the particle currents of mole- 
cules of type “”’ by J,', the continuity equations 
for a ternary mixture are? 


One 1 1 
J,{=—n |—-—} A 


On3(N—Nz 
Ox Bie Bae 


‘E. J. Hellund, Phys. Rev. 57, 319 (1940). 
*E. J. Hellund, “Osmotic diffusion in gases”’ (Phys. Rev. 
this issue), 


One(n—N2 Ne 
-— /na, (2) 
Ox Bis B32 


Bij = (ni+nj)Di;, (3) 
812812823 


In (1) and (2) the condition of constant total 
pressure has been used to eliminate the gradients 
of pressure of particles of type ‘‘1.”” Gases ‘3” 
and ‘‘2”’ in the following are to be treated as the 
isotope mixture which tends to separate under 
the influence of the diffusion of gas ‘‘1.”’ 

The initial conditions of the experiment are 
that the mixture ‘‘23” is contained in the “‘left’”’ 
half of a long tube of length ‘‘a” cm, and that 
gas “1” is contained in the right half. The dia- 
phragm which separates the two is removed and 
the resulting diffusion will partially separate ‘‘2"’ 
and ‘3,”’ the amount depending on the various 
coefficients of diffusion. 

If then 8:2 and 8,3 are assumed to differ by a 
small fraction of their values, one may consider 
that the solutions of (1) and (2) approximate 
closely the solutions obtained by assuming the 
mixture to diffuse as a simple gas. The diffusion 
coefficient of the mixture then would lie between 
the values for each isotope. The exact combina- 
tion to be used could be determined by requiring 
that the perturbation of (1) and (2) vanish if 
B23 were zero. It may be shown that this would 
give for the ‘‘zero-order’’ diffusion coefficient of 
the mixture the value, 


Let us define further, 
D* =823/n. 


(4) 
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Now the perturbations m2 and 6n3 must satisfy 
the same boundary conditions as mz and m3 and 
moreover vanish for =0, ©. In view of this the 
following approximation? may be made, 


cel? 
Ox ind | 7) 


Then one obtains from (1) and (2) two simul- 
taneous linear ordinary differential equations for 


a/2 a/2 
a= f “énedx and a=f bn3dx, 
0 0 


(8) 
de/dt (9) 


The coefficients a, 6 in (8) and (9) are depend- 
ent on the concentrations of gases ‘‘2” and “3” 
and the coefficients of diffusion D and D*. The 
C; and C, terms are likewise dependent on these 
variables and proportional in addition to the 
pressure gradient of gas “1” evaluated at x=a/2. 
The result of the diffusion, initially, is that the 
mixture ‘'2,3” diffusing into gas “‘1”’ has its con- 
centration changed from A, [A=12/(n2+n3) ], by 
the amount . 

(~) (10) 
D*+D\ 
where AD is the difference between the coeffi- 
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cients of diffusion of the mixture “3,1” ang the 
mixture “1,2.” 

In the time t=a*/r*D the expression (10) 
changes by the factor 8/(ex?—8), or approxi. 
mately 3. 

It is interesting to apply this result to the 
diffusion of He into air, which may be taken as 
Nz (80 percent), O2 (20 percent). The change jn 
concentration of the diffusing air, assuming a 
difference of one part in 200 for the diffusion of 


Ne and Os: into Hz would be of the order of 


1/2000. 

Thus the air mixture diffuses in this Case 
effectively as a simple gas. However, when the 
gas ‘‘1” is not light, the change may be much 
greater as, if the He should be replaced by some 
gas of large molecular weight and cross section, 
the change in the separation may be increased 
by a factor of 60 or 70. 

An alternative method of separation of jgo. 
topes employing the action of a third gas would 
be to place the isotope mixture in a tube closed 
at both ends by porous plugs which would be 
impervious to the isotope mixture, but permit 
the passage of some other gas which could there. 


fore be pumped through. The separation in this. 


case would increase with the increase in the 
diffusion current. 
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Note on Griineisen’s Constant for the Incompressible Metals 


J. C. SLATER 
Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 
(Received February 19, 1940) 


Griineisen’s constant y can be determined in terms of the 
thermal expansion, compressibility, and specific heat, as 
originally shown by Griineisen. It can also be found from 
the compressibility and change of compressibility with 
pressure. Theoretically the two methods should give ap- 
proximately the same result. This has proved to be the 
case experimentally, except for the more incompressible 
metals, where the discrepancies have been very great. 
Bridgman has now redetermined the change of com- 


Ik the frequencies of vibration of the normal 
modes of a crystal all vary as the inverse 7 
power of the volume, then it can be shown that 
there is a relation 


(thermal expansion) X (volume) 
(compressibility) X (Cy) 


(1) 


pressibility of iron with pressure, obtaining a much smaller 


' value than before, and as iron was used as a standard, this 


brings about a revision of other changes of compressibility 
with pressure, largest in proportion for the incompressible 
metals, whose volume change is about the same as for iron, 
It is shown that this revision is just enough to bring the 


two methods of finding y into agreement, for practically 


all the metals. 


This relation is due to Griineisen,' and y is often 
called Griineisen’s constant. By means of the 
relation (1), experimental values of y can be 
found, and for most simple solids they lie between 


1For discussion of this relation, see for instance J. C. 
Slater, Introduction to Chemical Physics (McGraw-Hill, 
1939), Chapter XIII, Section 4. 


~ 


Bri 
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1.5 and 2.5. In addition to (1), there is another 
quite different connection between y and experi- 
ment. This arises from the fact that the natural 
frequencies can be computed from the elastic 
constants, in particular from the _compressi- 
bility. The constant depends essentially on the 
change of frequency with volume, and hence on 
the change of compressibility with pressure. As 
a result of this, a value for y can be computed 
from Bridgman’s measurements of compressi- 
bility and its change with pressure. For many 
jonic crystals, and for the more compressible 
metals, the values computed in this way agree 
well with those found from Eq. (1), giving 
valuable verification of the theory, and a useful 
relation between experimental quantities, which 
could be used to find the thermal expansion from 
the change of compressibility with pressure, or 
vice versa. For the less compressible metals, 
however, the agreement is very poor, and to get 
better agreement, the change of compressibility 
with pressure should be considerably smaller (in 
the case of the least compressible metals several 
times smaller) than found by Bridgman.? 
Bridgman* has now redetermined the change 
of compressibility of iron with pressure, obtaining 
a value only about a third as large as the previ- 
ously measured value. This is important for all 
the other measurements, for iron is used as a 
standard, and the measurements on other ele- 
ments give really the difference between the 
change of volume of the element in question and 
iron. For the compressible metals, and the ionic 
crystals, the change of compressibility with 
pressure is so much greater than that of iron that 
a change in the value for iron makes only a small 
correction in the result, but a number of the 
metals of the transition groups have changes of 


compressibility with pressure of the same order - 


of magnitude as iron, and for these the new 
measurements make a very great difference. It is 
the purpose of this note to point out that 
Bridgman’s new measurements for iron entirely 
remove the discrepancy between the two 
methods of determining y for this element; and 
that if the old measurements of change of com- 

*For comments on this, see N. F. Mott and H. Jones, 
Properties of Metals and Alloys (Oxford, 1936), Chapter I, 
Section 4. omparisons of the two methods of determining 
7 are given in Slater, reference 1, Chapter XXIII, Section 
4, for the alkali halides; and Chapter XXVII, Section 2, 


for metals. 
*P. W. Bridgman, Phys. Rev. 57, 235 (1940). 
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pressibility with pressure for other incompres- 
sible metals are corrected by using the new 
results for iron, they likewise fall into line, with 
one or two exceptions. 

Bridgman’s‘ experimental results are stated 
in the form 


AV/Vo= —a,P+a2P%, (2) 


where AV the change in volume, V» the 
original volume, P the pressure, a; and az con- 
stants at a given temperature. For iron, at 30°C, 
he found in 1923 a;=5.87 X10~, ag=2.110-®, 
where the pressure is expressed in kilograms per 
square centimeter. In the new paper, the change 
of lengths rather than that of volume is given, 
but when the change of volume is computed from 
it (taking account of the fact that in the term in 
P* we cannot simply multiply by 3, but must 
cube the length to find the volume), we find 
a,=5.83X10-, a2=0.80X10-". That is, the 
value of a; is decreased by 0.04X10~, and that 
of a2 decreased by 1.30X10-". It is not hard to 
show that, to the accuracy with which the 
measurements are made, the measured values 
of a; and a2 for other metals are to be corrected 
by subtracting equal amounts from them. In this 
way, we have obtained corrected values from the 
observed values of reference 4. 

It is now easy to show’ that 7 can be found in 
terms of a; and a2 by the equation 


(3) 


In Table I we give values of a; and ae as 
originally found by Bridgman,‘ for most of the 
metals for which the most serious discrepancies 
existed, and the values of y computed from them 
by Eq. (3); next we give revised values of a; and 
2, as described in the preceding paragraph, with 


*P. W. Bridgman, Proc. Am. Acad. Arts and Sci. 58, 
165 (1923), for most of the metals mentioned in the 
present note. 

5 See Slater, reference 1, Chapter XIII, Section 1; and 
Chapter XIV, Section 4. In Chapter XIII, Eq. (1.1), (1.4), 
etc., neglect a, which represents the thermal change of 
volume. Then a), a2 of the present paper become identical 
with the quantities defined in that chapter. In Eq. (4.6), 
cmemnee IV, it is shown that vy=—4+P:/P,. In Eq. 
(1.10), Chapter XIII, neglecting ao, it is shown that 
P,=1/a;, P:=a;/a,*. From these the result of Eq. (3) of 
the present poper follows at once. This result does not 
with Eq. (42) of Mott and Jones (reference 2) which in our 
present notation is y=a;/a;*—}, or } unit greater than 
our value. We believe Mott and Jones to be wrong, for 
they have incorrectly defined the compressibility as 
—(1/Vo)(@V/aP)r, rather than —(1/V)(@V/aP)r, and 
have thus lost a term coming from the change of V with 
volume, leading to the difference between the formulas. 
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TABLE I. Original and revised values of ai, a2, and y. 
Data for a; and a, are from Bridgman, reference 4. Values 
of ¥ by Griineisen’s method [y(Gr)] are from the article by 

inetsen, ‘‘Zustand des festen Kérpers,” in Handbuch der 
Physik, Vol. X (Springer, 1926). 


a1X10? a2X10% a1X107 ~a2X10" 

METAL (OLD) (OLD) (OLD) (REV.) (REV.) (REV.) (Gr) 
Fe 5.87 2.1 54 583 080 1.68 1.60 
Co 5.39 21 66 535 0.8 2.1 1.87 
Ni 5.29 2.1 69 5.25 08 2.2 1.88 
Cu 7.19 2.6 44 7.15 13 1.9 1.96 
Pd 5.19 2.1 7.1 5.15 08 2.4 2.23 
+ 987 44 38 983 3.1 2.5 2.40 
2.93 1.5 17. 2.89 0.2 1.7 1.62 

Pt 3.60 18 13. 3.56 0.5 3.3 2.54 


a corresponding revised 7; finally we give values 
of y computed by Griineisen’s relation, Eq. (1). 
It is plain from the table that the revision has 
completely removed the discrepancies between 
the old computed values of y and Griineisen’s 
values of that quantity. The agreement between 
the revised y’s and Griineisen’s is within the 
experimental error of the a2’s, which are only 
known to a very rough approximation. We must 
remember in connection with this that an exact 
agreement is not required by the theory. It would 
be expected only if all natural frequencies of the 
solid varied in the same ratio with change of 
volume, and if the change of frequency with 
volume arose only from the change of compres- 
sibility, and not from a possible change in 
Poisson’s ratio. Neither of these conditions is 
likely to be exactly fulfilled, so that the two 
methods of calculating y may well differ by as 
much as the observed amounts. 

For the more compressible metals, the revised 
values of az lead to only a small change in y, so 
that the agreement is not essentially changed 
from the values given in Table XXVII-2.! 


There are a few incompressible metals, however, 


for which there are outstanding discrepancies. 
For manganese, Bridgman’ finds a,;=7.91 X10-’, 


*P. W. Bridgman, Proc. Am. Acad. Arts and Sci. 64, 
51 (1929). 


C. SLATER 


a2=5.3X10-"; corrected values would be g 
=7.87X10~, a2=4.0X10-, from which 
Griineisen’s value is 2.42. Since this discrepancy 
is so much larger than for other metals, there is 
good reason to suspect that Bridgman’s value of 
az is considerably too large. This is perhaps ex. 
plained by the fact that he states that jt js 
difficult to get manganese without small cracks, 
on account of a polymorphic transition; a 
cracked specimen would show an apparently high 
value of a:. For molybdenum and tantalum, 
there are discrepancies in the opposite direction, 
Bridgman finds a2 equal, respectively, to 12 
X<10-” and 0.25 X10-" for these metals. Maki 
our correction of subtracting 1.3X10-" woul 
give negative values of de. This seems most 
implausible on general grounds, so that the 
observations for these two metals are certainly 
to be questioned. Finally for gold he used two 
samples, one giving @2=2.1X10-", the other 
3.1X10-". These lead to y=1.8 and 4.8, re. 
spectively, against Griineisen’s value of 3.03, 
Evidently the order of magnitude is right, but 
the discrepancies between the two observations 
are so great as to throw a good deal of doubt on 
the accuracy of the measurements. A remeasure- 
ment of these doubtful metals, if not of all ‘the 
incompressible ones, would be very desirable. 
Since it appears as a result of the present note 
that Griineisen’s relation is really rather gener- 
ally satisfied, it would be justified in future 
measurements to question any values of change 
of compressibility with pressure which disagree 
with it violently. For convenience, we may note 
that for most materials y+ is in the neighbor- 
hood of 2.5. Thus, from Eq. (3), we may expect 
that approximately 


ag= 2.5a r. (4) 


This relation should be expected to give the 
quadratic term in the volume change, in terms 
of the linear term, with an error of usually not 
over fifteen or twenty percent. 
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LETTERS TO THE EDITOR 


Prompt publication of brief reports of important discoveries in physics may be secured by 
addressing them to this department. Closing dates for this department are, for the first issue of the 
month, the eighteenth of the preceding month, for the second issue, the third of the month. Because of / 
the late closing dates for the section no proof can be shown to authors. The Board of Editors does 

not hold itself responsible for the opinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


Cross Section of the Reaction N“(p,a)C” 


Of the four reactions in the “carbon cycle” which is 
assumed to supply the energy in ordinary stars,' the three 
p-v reactions have been investigated and their cross 
sections measured. The reaction C¥%(p,y)N™ has been 
studied by Roberts and Heydenburg,? C"(p,7)N™ by 
Dee, Curran and PetrZilka,? and by Curran (1) 


and Struthers.‘ where X is the wave-length of the incident protons, « the 
This letter is a report of the measurement of the last oss value of its “wave number” at the eurface of 


reaction of the cycle; vis. the nucleus, the range of the nuclear forces, the 
Gamow penetrability factor, and S a factor taking into 
reported work on the deuteron disintegration of nitrogen. 


The protons striking the aluminum foil separating the 
target gas from the cyclotron tank had an a al of e2 = 430, assuming a nuclear radius R=1.5X10-2X16) 
Bellamy ‘and. the stopping power of cm which corresponds to a barir height of 
erb, . 2.6 Mev. S can calculat rom the rules given by 
aluminum Bethe" and comes out about 1.25. Therefore the theoretical 
8 percent less than lor cross ion f 1) is 0.9X 10~*6 cm!, in very satisf. 
calibrating the foils. Therefore the foil used, of 0.865 cm 
air equivalent poy For a comparison with other astrophysical reactions, 
equivalent of f it may be most convenient to calculate the “effective 
He width” used in reference 1 Eq. (11). From the measured 
cm in P=1.0X107 
protons in the center of the target chamber was then Np 
0.36 Mev. The ionization eee oy ga The lifetime of a N"* nucleus at the center of the sun is 
re om steradians of the total solid angle about jy45 only about 20 years, provided no resonance effects 
the one taking place the range of the alphas was measured probably of the order 10° sateen Gn 0 tet cael tn 
and be be reference 1 (p. 446), which means that in each completed 
expected from e Q va carbon cycle only about one nucleus in a million is lost to 
mass pre Q value o Y the cycle. An account of the astrophysical consequences 
Burcham and Smith,® viz., 5. ev, range Of will be published in the Astrophysical Journal. 


would be expected. This is a sufficiently good check We wish to thank Professor H. C. Urey for the sample 
considering the uncertainty in the value for the proton of N" gas. : 


energy. A target of N" alone gave no alphas. This last 


N"5(p,a)C" reaction is calculated to be 1.3 cm? at 
0.36 Mev. 

Theoretically, we may safely assume that the p,a 
reaction is by far the most probable reaction which can be 
caused by the collision of a proton with N". Then its 
cross section should be® 


M. G. HoLtoway 


fact and the reasonable check of the range make it virtually De _— 3. A. Bec ' 
certain that the alphas came from the N'*(p,@)C" reaction. en Ce 

The yield was obtained by taking a number distance March 13, 1940. 
curve with a thyratron bias low enough to count ail alpha- 1H. A. Bethe, Phys. Rev. 55, 434 (1939). 


particles having 3 mm or more of their range in the ioniza- * R. B. Roberts and N. P. Heydenburg, Phys. Rev. 53, 374 (1938). 
tion chamber. The current of protons (~0.3ma) striking «S.C. Curran and J. E. Struthers, Nature 145, 224 (1940). 
the target was integrated by means of a circuit described * This reaction was reported, but no cross section measured by 


by Herb, Kerst and McKibben.* The total yield in 4 * M. G. Holloway and B. L. Moore, Phys. Rev. 56, 705 (1939). 
solid angle was 2.3X10~* alpha per proton at 0.36 Mev Phy Rev $2.78 (1037) ©. J. C. and ©. 


with a target containing 3.78X10' N atoms/cm® at go 8,0; Herb, D. W. Kerst, and J. L. McKibben, Phys. Rev. 51, 


691 (1937) 
0°C, 760 mm Hg and 0.046 cm thick (air equivalent at *H. A. Bethe, as continnae theory of the compound nucleus, to ’ 
0°C, 760 mm). From these figures the cross section of the *"Ry 4" Bethe, Rev. Mod. Phys. ° , 178 (1937). ' 
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Further Experiments on Fission of Separated 
Uranium Isotopes* 


Recently! we bombarded with slow neutrons small 
quantities of uranium isotopes which had been separated 
by means of a mass spectrometer. Measurement of the 
fission rates from the samples gave strong support to the 
view that it is U* (or possibly U™) rather than U8 which 
undergoes fission under slow neutron bombardment. 

This conclusion has now been verified with considerably 
larger samples of separated uranium isotopes. In addition, 
it has been possible to obtain tests of several samples of 
partially separated U** and U™. In the new work the 
uranium ion currents were increased by a factor of about 
ten times by increasing the electron current bombarding 
the UBr, vapor by a factor of ten. The U** ions were col- 
lected on an insulated platinum plate (115 mm) and the 
current to this plate measured with an electrometer tube 
amplifier. The U* and U™ ions were collected in the first 
set of samples on a single grounded platinum plate and in 
the second set on a pair of grounded plates. The dispersion 
of the apparatus is such that for infinitely great resolving 
power, ions differing in mass by one unit would fall along 
sharp parallel lines on the collector separated by 0.75 mm. 
From a measure of the variation of the current to the 
insulated collector with energy of ions, it was possible to 
determine the actual spread of the ion beam and hence 
the overlap of adjacent isotopes. It was concluded that 
for a region 0.75 mm wide, centered about the U™ position, 
one would collect approximately 25 percent as many U** 
ions as in an equally wide strip centered about the U** 
position. As the ion peaks are nearly symmetrical the 
converse also holds. That the resolution was good enough 
to separate U™ and U** was not apparent in the earlier 
work as the insulated collection plate was considerably 
wider (2.0 mm). 

Slow neutron tests —The separated isotopes were first 
tested by placing in an ionization chamber connected to a 
linear amplifier, and bombarding with neutrons from the 
Columbia cyclotron which had been slowed in paraffin. 


TABLE I. Slow neutron tests of separated U™, U™5 and U™, 


Us us 
Sample set No. 1: Mass U™* =3. ipg 
alphas/min. expected 1.1 0.05 1.1 
alphas/min. observed 1.0+0.1 0.25 +0.05 0.75 +0.1 
fissions/min. 0.1+0.1 4.4 +0.2 1.5 +0.1 
fission/alpha ratio 0.1 18. S. 
Sample Set No. 2: Mass U8 =4.4yug* 
alphas/min. expected 1.6 0.07 1.6 
alphas/min. observed 2.3+40.1 0.63 +0.1 2.2+0.2 
fissions/min. 0.5+0.2 10.3 +0.5 7.8+0.5 
fission/alpha ratio 0.22 1 3.5 


* Considerable difficulty in obtaining large ion currents was encoun- 
tered in this run. As a result, a great deal of testing and adjustment was 
necessary. Unfortunately, no record of the time or currents was kept 
while adjustments were being made. The mass of the U™* target was 
thus larger than 4.4ug. This would account for a good share of the dis- 

between the and observed alpha-counts. The resolu- 


crepancy 
tion of the apparatus was definitely poorer in this run than the number 
one run. 
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TABLE II. Fast neutron tests of separated U8, 


U8 No. 2 U™8 No, 3 
alphas/min. 2.3+40.1 
fissions/min. without Cd 0.5+0.1 a3 
fissions/min. with Cd 0.5+0.1 1.4402 
— 


For the first set of samples, where the U™ and Ys were 
deposited on the same plate with some overlapping, the 
alpha-particle and fission counts were made by : 
each of the sample regions in turn with aluminum foil 
0.008 cm thick. 

Table I shows (a) the number of alpha-particles/mip, to 
be expected on the basis of the masses of the uranium ig, 
topes deposited as calculated from the mass spectromete 
ion currents assuming complete separation, (b) the ob. 
served number of alpha-particles/min., (c) the number of 
fissions/min., and (d) the fission/alpha-ratio, for the yagi. 
ous isotope samples. The background has been subtracted, 

The number of alpha-particles/min. observed js jp 
agreement with the expected number within the limits ¢ 
error. The amount of overlapping of U™ and U** is likewig 
close to that calculated from the resolution of the mas 
spectrometer. 

The number of fissions/min. observed and the fission 
alpha ratios show conclusively that U** is responsible for 
at least seventy-five percent of the slow neutron fission 
observed in unseparated U. Since U™ is present only to{ 
part in,17,000, it cannot be excluded that U™ does undergo 
some slow neutron fission, but the total contribution ¢ 
U™, if any, is small. 

Fast neutron tests —Further tests were made of Ut 
see if the small numbers of fissions observed from it in the 
slow neutron tests, were due to fast neutrons. Measure. 
ments were made without paraffin, and with and withow 
Cd surrounding the fission chamber, on the U™* samp 
from set No. 2 and also on a much larger sample, No, 3 
(18.9 micrograms), subsequently collected under improved 
conditions. The data are recorded in Table II. 

These results indicate that U** undergoes fast neutron 
fission only. Furthermore, the number of fissions/micn 
gram of U** observed under these neutron intensity co 
ditions, is sufficient to account for practically all the fas 
neutron fission observed from unseparated U. 

The aid of the Research Corporation to Columbia 
University is gratefully acknowledged. The collection ¢ 
the isotopes was supported by a grant from the University 


of Minnesota Graduate School. 
A. O. 
Department of Physics, 
University of Minnesota, 
Minneapolis, Minnesota, 
E. T. Boors 
J. R. Dunning 
A. V. Grosse 
of Physics, 
olumbia University, 
New York, New York, 
April 3, 1940. 


* Publication assisted by the Ernest Kempton Adams Fund ie 
Physical Research of Columbia University. 
1 Nier, Booth, Dunning and Grosse, Phys. Rev. 57, 546 (1940). 
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LETTERS TO 


Shower Production by Mesotrons in Different Materials! 


The apparatus used has already been described* and 
comprises a number of Geiger counter trays, each forming 
a sensitive area 20 by 20 cm, mounted in a vertical column 
and arranged in relation to the associated electrical circuits 
in such manner that the apparatus is only allowed to 
record when a ray passes through all of the trays. Below 
the top tray is placed a block of lead 18 cm thick, or its 
mass equivalent in other materials, so that all records by 
the apparatus are initiated by a ray which is capable of 
passing through at least this thickness of material. When a 
ray passes through all of the trays, any shower rays ac- 
companying it are recorded separately by the individual 
counters which they excite, the record being a photograph 
of spots of light from the mirrors of a number of electro- 
ne for each counter. The whole apparatus is 


aed below a 30-foot column of water contained in a 


large water tank. 

With the above apparatus the authors measured the 
numbers of one-ray showers, two-ray showers, three-ray 
showers, etc., produced in lead last August, and the present 
experiments extend the measurements to lead, tin, iron, and 

um. 

In dealing with the various materials, it was necessary 
to have a thickness sufficiently great to insure that shower 
production through mesotrons had attained equilibrium. 
Thus, it was necessary to have the thickness comparable 
with or greater than the range of the electrons. In order 
that one might be free from complications involved in 
variation of mesotron intensity by different absorptions 
in the different elements, it was arranged that, in the case 
of each element, the total equivalent thickness as regards 
mass absorption should be the same, the balance being 
adjusted for convenience in each case by the choice of a 
suitable thickness of lead above the material under 
examination. 

Figure 1 shows the results. Horizontally are plotted the 
atomic numbers of the elements concerned. The ordinates 
represent the numbers of showers, of the kind cited, asso- 
ciated with the passage of 1000 mesotrons through the 


apparatus. 
The two-electron, and perhaps the three-electron, 
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showers showed frequency of occurrence which is inde- 
pendent of the atomic number within the limits of accuracy 
of the experiments; but the one-electron showers show a 
marked dependence upon atomic number, a dependence 
sufficiently strong to dominate the situation even if we 
should plot a curve representing the total number of shower” 
rays associated with 1000 mesotrons. These results are not 
in harmony with the elementary conclusions from Bhabha’s 
theory, which predicts results practically independent of 
atomic number. As shown in our former paper,’ however, 
the data for the actual numbers of showers agree, as 
regards order of magnitude, with the theoretical predic- 
tions. 

It is a curious fact that the regularity of the graphs for 
the actual data is greater than one might expect on the 
basis of the statistical uncertainty indicated in the usual 
way by the vertical. This is a matter which may have a 
rather profound significance in relation to processes 
involved. 

W. F. G. Swann 


W. E. RAMsEY 


Bartol Research Foundation, 
The Franklin Institute, 
Swarthmore 


March 26, 1940. 


1 Presented at the New York Meeting of the American Physical 
Society, February 23-24, 1940 (Phys. Rev. 57, 554(A) (1940)). 

2W. F. G. Swann and W. E. sey, Phys. Rev. 56, 378 (1939); 
W. F. G. Swann, Rev. Mod. Phys. 11, 242 (1939). 


Fission of the Separated Isotopes of Uranium* 


Samples of the separated isotopes of uranium have 
been prepared with a mass spectrograph. UCI, was 
vaporized in a furnace, and the electron source for ion 
formation was placed inside the furnace. The upper 
entrance slit of the analyzer was also mounted on the 
furnace to avoid condensation of UC. With a beam 
current of U***+ of 6X10-* amp., about 1.8 microgram of 
U*™* was collected on a Pt strip in three hours. The corre- 
sponding amounts of U**+** were collected on an adjacent 
Pt strip. A third Pt strip, mounted on the opposite side of 
the U** collector, showed no measurable a-emission. The 
a-count for the U** sample was 0.6+-0.1 per*minute, and 
for the U****™ sample 0.6+0.1 per minute. 

The fission of the separated isotopes was tested by 
bombarding with slow neutrons from the Columbia 
cyclotron. The U™*#5 sample gave 3.7+0.4 fissions/ 
minute, and the U** sample gave 0.1+0.1 fission/minute. 
This shows definitely that U** is not responsible for slow 
neutron fission in uranium, and confirms the results 


K. H. Kincpon 
H. C. Pottock 


Comper. 


E. T. Bootn 


J. R. Dunninc 
Columbia University, 
New York, New York, 
March 20, 1940. 
* Publication assisted by the Ernest Kempton Adams Fund for 


ysical Research of Columbia University. 
1 Nier, Booth, Dunning and Grosse, Phys. Rev. 57, 546 (1940). 
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750 LETTERS TO 


Collisions of Alpha-Particles with Chlorine Nuclei 


The cloud tracks of about 700,000 thorium (C+C’) 
alpha-particles have been photographed with a stereo- 
scopic camera. A large-diameter cloud chamber filled with 
methyl-chloride and helium was used. In addition to the 
range and velocity measurements obtained for carbon recoil 
nuclei,’ forty-one selected forks ascribed to collisions of 
alpha-particles with chlorine nuclei have been analyzed, 
the results of which are presented here in Fig. 1. The 
analysis was carried out with the same procedure as that 
used for the carbon collisions by one of us (G. A. W.), the 
value 36 for the atomic weight of chlorine being used 


throughout. 
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Fic. 1. Range-velocity curve for chlorine recoil atoms in methyl-chloride 
and helium. 


The range-velocity curve for recoil fluorine nuclei as 
obtained by Feather? is included for comparison (broken 
curve). No outstanding similarity is observed in the shape 
of the fluorine and chlorine curves for the region of recoil 
ranges considered, although such might have been ex- 
pected from Anthony’s* observations. The range-velocity 
curve for chlorine obtained from that for fluorine by the 
empirical relation given by Blackett and Lees‘ gives poor 
agreement with the present results. 

HANSEN 


Sloane Ph boratory, 
Yale etm GERALD A. WRENSHALL 
New Haven, 


March 26, 


1G. A. Wrenshall, to appear soon in the Physical Review. 
Feather, Proc. Roy. Soc. 194-200 (1933). 
aR. L. Anthony, Phys. [2] 50, 726-732 (1936). 


4P. M. S. Blackett and D. S. , Proc. Roy. Soc. Al34, 658-671 


A New Radioactive Isotope of Masurium, ,;Ma’'* 


An accurate study of the 24-min. period produced in 
molybdenum by bombardment with slow neutrons made 
it clear that this period should be analyzed into two. 
This point was first noticed in a very short exposure, in 
which the decay curve showed a distinct convex curvature. 
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Fic. 1. The decay curves obtained for the chemically separateq 
masurium and molybdenum fractions. 


A very quick chemical separation of the molybdenum 
sample into the masurium fraction and the molybdenym 
fraction after a slow neutron bombardment, enabled us to 
find a new masurium isotope, Ma’ which emits B-rays 
of about 1.14 10* ev in energy (K-U upper limit) witha 
half-life of 9+1 min., while the mother isotope Mola 
decays with a half-life of 19+1 min. emitting B-rays of 
1.78 X 10® ev in energy. The decay curves for the chemically 
separated molybdenum and masurium fractions are given 
in Fig. 1. 

Another masurium isotope Ma” has been found by 
Seaborg and Segré,' the half-life of which is 6 hours by 
it emits very soft 8-rays.** For the investigation of some 
special chemical nature of masurium, it would be much 
more convenient to use Ma’ rather than Ma™ because 
it is produced very strongly by bombardment with sly 
neutrons and also it emits energetic 8-rays. 

We are indebted to Dr. Y. Nishina and Professor K 
Kimura for their valuable help and discussions. The 
assistance given by Mr. Matsuura, Mr. Saito and Mr. Mori 
is highly appreciated. We wish to acknowledge the financial 
support given by the Japan Society for the Promotion of 
Scientific Research, Oji Paper Manufacturing Company, 
Mitsui Ho-Onkwai Foundation, Tokyo Electric Light 
Company and Japan Wireless Telegraph Company. 
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S. Kojima 
G. Miyamoto 
M. IKAWA 
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‘okyo Imperial University, 
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* Reported at the meeting held on December 13, 1939 at the Institute 
of Physical and Chemical Research. 
1G. T. Seaborg and E. Segré, Phys. Rev. 55, 808 (1939). 
2 By bombardment with protons, ner periods, 0.5 min., 31 mia, 
53 min. .~ 110 hr. have been found, but not much has been 
me Du Bridge, Barnes, Buck and Strain, Phys. Rev. 53, 447 (1938); 
herr, Henderson, White, Delsasso and Ridenour, Phys. Rev. 53, 4 
(1938); Ewing, Perry and McCreary, Phys. Ly 59, 1136 (1939). 
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On the Dynamics of Complex Fission 

We wish to report preliminary results of an investigation 
of the large amplitude distortions of a nucleus undergoing 
fission. The nucleus is assumed to have a uniform charge 
and mass density and a constant volume for all distortions. 
The radius vector r is expanded in a series of Legendre 

ynomials r=R(i+2a,P,), and the values of ao and a; 
are adjusted to insure constancy of the volume and center 
of gravity. An estimate of the quantum excitation energies 
shows that only the first five harmonics need be considered. 
The potential energy of distortion, including all terms so 
far evaluated, is given by 


AE —x)a2* + (0.7143 —0.4082x)a;? 
+(1—0.3704x)ae+ (1.2727 —0.3306x)a;? 
—(0.03810+0.07619x)a.°— (0.07619 +-0.2503x)aza;* 
— (0.1143 +0.3429x)as%a4— (0.05195 +0.2226x)a3*a4 
— (0.1732 (0.2171 —0.2563x)a2* 

— (1.3615 — 1.0934x)as*as*+ (0.09350 +0.03206x)a2*}, 


where 2x = (32%e*/5.R) /(42R?O) and O is the surface tension. 

The following conclusions may be drawn from this 
formula. (1) There are large coupling terms between odd 
and even harmonics. These provide a possible explanation 
of the marked asymmetry of fission.’ Bohr and Wheeler* 
have neglected these coupling terms. (2) The coupling 
between az and a, is such as to give an equatorial bulge 
to the drop when it is on the road to fission.* The drop is 
pinched in near both poles; this would lead to ternary 
fission were it not for the excitation of the third harmonic 
which causes one of the preformed droplets to come off 
first, leaving a residual fragment which has not a sufficiently 
high excitation energy for its low charge number to divide 
further. However, it seems possible that such a ternary 
fission may occasionally be observed, in view of the fact 
that the saddle point for this mode of fission corresponds 
to a;=0, and moreover, appears to lie at very low energies. 
Thus, whether fission occurs principally through the coup- 
ling of a, and a; or through the coupling of a; and ay, in 
either case the fragments are of unequal sizes. According 
to the experiments, it is very improbable that the nucleus 
should divide into nearly equal fragments." * (3) The 
coupling between a: and a; leads to an energy surface 
which has been investigated for a wide range in x. We have 
assumed a; to be of the order of a2! and taken all terms in 
AE up through the fifth order in a2. For x=0.75 we find 
a saddle point at a.=0.66, a;=0.30 with an activation 
energy E;=13 Mev (assuming the surface energy of 
uranium to be 500 Mev); for x=0.85 there is a saddle 
point with a,=0.65, as=0.25 and Mev. For such 
large amplitudes, the terms considered can give only a 
rough description of the fission process; moreover, theo- 
retical uncertainties in the values of the nuclear radius and 
surface tension make any exact estimate impossible. These 
results are nevertheless of the right order of magnitude to 
explain the observed fission of uranium. (4) The energy 
surface for the a2—a, coupling has been investigated for 
terms up through the fourth order in a; assuming that a, 
is of the order of a;*. For x=0.75 there is a saddle point 
at @2=0.47, a4=0.057 with E;=4.4 Mev and for x=0.85 
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we find one at a;=0.31, a4=0.029 with Ey;=1.1 Mev. 
(5) With all the mutual coupling terms for a2, a3 and a, 
included, there are two roots of the extremalizing equations 
which are of interest. One of these corresponds to a; =0 and 
reproduces the results already stated in (4). The other 
leads to a reduction of the amplitudes and activation 
energies evaluated in (3). We have not yet included a 
sufficient number of terms to give a significant estimate of 
this effect, but it appears to be of the order of several Mev. 
In view of this fact and especially because of the modifica- 
tions that may arise with the inclusion of higher order 
terms, we cannot conclude that the a.—a, mode of fission 
requires less activation energy than the a;—a3 mode. (6) An 
estimate of the particle sizes to be expected from the a;—a; 
mode of fission can be made by using the relative ampli- 
tudes of ai, a2 and a; at the saddle point. The two fission 
products have approximately 0.4 and 0.6 the mass of the 
original nucleus, this result being insensitive to the value 
of x. For the a2—a, mode of fission the particle sizes are 
roughly in the ratio of 1:3 for binary fission and 1: 2:1 
for ternary fission. The first of these is consistent with the 
experimental energy peaks." 

Finally we want to emphasize that certain conclusions 
of this note are subject to possible modification due to the 
inclusion of higher order terms which are probably not 
small. A more detailed discussion of the dynamics of fission 
will be published shortly.’ 
R. D. PRESENT 


J. K. Knipe 
Purdue University, 
Lafayette, Indiana, 
April 1, 1940. 


1 Booth, Dunning and Slack, +4 ~ Rev. 55, 981 (1939); M. H. 
Kanner and H. H. MS , ote. v. 57, 372 (1940). 

2?N. Bohr and J. A. Wheeler, Phys. Rev. 56, 426 (1939). 

3 In contradiction to the statement of Bohr and eeler, reference 2, 
p. 433, that a4 contributes a a about the equatorial belt such 
as to lead continuously to a dumb-bell shaped figure. We do not check 
the higher order coefficients of their Eq. (22) and we note further that 
Eq. (23) is inconsistent in sign and magnitude with Eq. (22). 

*L. A. Turner, Rev. Mod. Phys. 12, 1 (1940), cf. table on p. 23. 

* Related biophysical problems are being considered. ‘ 


Radio Isotopes of Chromium 


We have investigated the induced radioactivities of 
chromium by the aid of our cyclotron. The results are 
summarized as follows: 

(I) Deuteron bombardment. The exposure of chromium 
metal to about 100 microampere hours of 3-Mev deuterons 
yields chemically identified chromium isotopes of half- 
lives 1.6 hours and about 14 days. 

(II) Slow neutron bombardment. In this case activity 
induced in chromium was very weak. The following 
periods were established : 2.8 hr., 14 hr., 1.7 hr. The purity 
of the sample was chemically tested and Al (1.1 percent), 
Si (0.17 percent), Fe (0.7 percent), Mn(>0.02 percent) 
were shown as impurities. Therefore the 2.8-hour activity 
seems to be due to Si* (2.5 hours) or Mn** (2.5 hours), 
and the 14-hour activity is probably due to Na™ produced 
in aluminum impurity by the reaction Al-n-a-Na, as a 
considerable number of fast neutrons were always among 
slow neutrons. The period 1.7 hours is very near to the 
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period 1.6 hours observed in the case of deuteron bom- 
bardment, and therefore can probably be ascribed to an 
isotope of chromium, but unfortunately no chemical iden- 
tification was done because of very weak intensity. 

(III) Fast neutron bombardment. Periods of 3-4 hours, 
14 hours and about 12 days were observed. In this case the 
activities were very weak. The 3—4-hour activity is very 
likely due to V® (3.7 hours) according to the reaction 
Cr®-n-p. The 14-hour activity is probably Na™ as in the 
case of slow neutron bombardment. The period 12 days is 
probably due to the same activity (14 days), as was ob- 
served in the case of deuteron bombardment. 

Now chromium has four stable isotopes Cr®*, Crs, 
Cr, According to the above results we can attribute the 
1.6-hour activity to Cr®, and the longer period (about 14 
days) activity to Cr®. 

A more detailed account of this work will be given later 


in the Scientific Papers of the Institute of Physical and | 


Chemical Research. 
In conclusion we wish to express our best thanks to 


Dr. Y. Nishina for his kind interest. 
TosHio AMAKI 


TAKEO IIMORI 
Asao SUGIMOTO 


Nuclear Research 
Institute of Physical and Chemical Research, 


Tokyo, Japan, 
January 10, 1940. 


Spectroscopic Evidence for the B, Molecule 


Modern valence theories lead one to expect that a 
molecule formed from two boron atoms should be stable.! 
The stability of B:H¢ also suggests that a stable bond 
between two boron atoms is possible. But no evidence for 
the existence of the B: molecule has yet been reported. The 
observation of a Bz spectrum, apart from showing that B: 
exists, also would be of interest because it would make 
possible a determination of the spin and the statistics of 
the boron nuclei. 

We have therefore tried to observe a Bz spectrum in a 
discharge through helium of about 10 mm pressure to 
which a trace of BCl; was added. Under these conditions 
we have been able to find between 3300 and 3170A a 
system of new bands having a simple fine structure (P and 
R branches only). The following considerations prove 
definitely that the emitter of these bands is the Bz molecule 
whose existence is thus established for the first time. 

(1) A clear intensity alternation is observed in the two 
bands that are least overlapped by others, the even- 
numbered lines being the strong ones. This shows definitely 
that the carrier of the spectrum is a homonuclear molecule. 

(2) Each band is accompanied by a weaker isotopic 
band, whose intensity and displacement from the main 


THE EDITOR 


band is accountable only on the assumption that the main 
bands are due to B"B", the weaker bands to BURw The 
B"B" bands, as is to be expected, do not show any jn. 
tensity alternation. There are indications also of bands dye 
to BB, 

Only two sequences of the new system, Ap=( and 
Av= +1, have thus far been observed. A preliminary vibra. 
tional analysis yields the following formula for the band 
heads: 


v= 30546. 1+ (929.3v—2.75v*) — (1035.20” —9,58y"), 


A preliminary rotational analysis yields Bo” = 1,205 
Bo'=1.155 cm=! and consequently for the internucles, 
distances ro” = 1.595 X 10-8 cm, ro’ = 1.628 X 10-8 cm, 

A rough estimate leads to a value of the intensity ratio 
of the strong and weak lines between 2:1 and 15 | 
which would mean that the spin of the B™ nucleus js very 
probably J = §. However, more accurate intensity measure. 
ments are necessary to fix definitely this value. By such 
measurements we also hope to determine the spin of the 
B’® nucleus. 

The fact that the bands have only one P and one R 
branch shows that they represent a © —= transition, Since 
the lines with even values of J are the strong ones it follows, 
if Fermi statistics is assumed for the B" nucleus, that the 
lower state is a 2, state and the upper state 2.~. From the 
electron configuration one would expect the ground state 
of to be and a transition to bea 
one. It seems very likely that the new bands are this trap. 
sition, the triplet splitting being too small for detection, |t 
is significant to find that the assumption of Fermi statistics 
for the B™ nuclei leads to reasonable results since it js the 
statistics one would expect for this nucleus on the basis of 
modern nuclear theory. 

When nitrogen was added to the discharge in He and 
BCI; two new band systems were found, one around 
3600A, the other around 3225A. The former can be 
assigned definitely, the latter probably to the BN molecule 
whose spectrum has not been previously observed. When 
hydrogen was added two new electronic transitions of the 
BH molecule were found at 3415 and at 3100A. 

A more detailed report of our experiments will be sub- 
mitted to the Canadian Journal of Research. 

This investigation has been supported by a grant from 
the Penrose fund of the American Philosophical Society and 
by the National Research Council of Canada. To both 
organizations we would like to express our thanks. 


A. E. 
G. HERZBERG 
~ of Physics, 
niversity of Saskatchewan, 
toon, Saskatchewan, Canada, 
March 25, 1940. 


1See, for example, G. Herzberg, Molecular Spectra and Moleale 
I (Prentice-Hall, New York, 1939), p. 392. 


